Perception, Navigation, and Control for Unmanned Aerial Vehicles: Theory and Applications
View this Special IssueResearch Article  Open Access
Yunjie Wang, Chen Jiang, Shuai Deng, Haowen Wang, "Helicopter Autorotation Trajectory Planning Method Using Functional TensorTrainBased Dynamic Programming Algorithms", International Journal of Aerospace Engineering, vol. 2020, Article ID 1372052, 13 pages, 2020. https://doi.org/10.1155/2020/1372052
Helicopter Autorotation Trajectory Planning Method Using Functional TensorTrainBased Dynamic Programming Algorithms
Abstract
Helicopter autorotation trajectory planning problems have been dealt within computationally expensive optimal control algorithms. This paper presents an efficient helicopter autorotation trajectory planning method, using functional tensortrain (FT) based dynamic programming (DP) algorithms. The autorotation trajectory planning method is shown realtime feasible, which involves general helicopter autorotation dynamics at the same time. To validate the dynamic feasibility of the trajectories, a trajectorytracking controller using active disturbance rejection control (ADRC) is designed to ensure a helicopter model tracks the trajectories. Finally, a helicopter autorotation simulation with a sixdegreeoffreedom highfidelity multibodybased helicopter model is demonstrated for validation.
1. Introduction
Autorotation is a primary measure for helicopters, whether for manned or unmanned helicopters, to land safely after engine power failure. In autorotation, the main rotor is driven by the upward airflow through the rotor, making the flight similar to a gliding fixedwing aircraft. To achieve a safe autorotation landing, unique control strategies are needed [1]. Generally speaking, the collective pitch should be carefully handled to maintain a sufficient and steady rotor speed before reaching the ground.
Studies on means of achieving autorotation have been focused on solving optimal control problems [2–4]. Traditionally, simplified helicopter dynamic models are used, including 2D pointmass models [4–6], threedegreeoffreedom rigidbody models [7, 8], and loworder sixdegreeoffreedom rigidbody models [3]. With the dynamic model formulated, optimal autorotation problems are solved by numerical methods. Gradient algorithms, such as sequential gradientrestoration algorithms, are used in [5, 7, 9]. Other methods, such as direct methods, discretize the problem first and turn it into a nonlinear programming problem. Then, the nonlinear programming problem can be solved using algorithms such as sequential quadratic programming (SQP), as demonstrated in [6, 10, 11]. In general, these algorithms are offline, computationally expensive for current computing power, which is for now unrealistic for online usages. Furthermore, there is a lack of validations either through highfidelity simulations or experiments in the above studies.
Recently, Taamallah [12] proposed the first realtime feasible, modelbased trajectory planning method and designed a modelbased trajectorytracking controller to ensure the helicopter tracks the trajectories. In specific, they employed optimal planning based on differential flatness, assuming a helicopter as a rigid body. The trajectory planning problem is solved, regarding all the rotor forces and moments as plant inputs of the rigid body. Taamallah’s work provides novel directions for online trajectory planning of autorotation problems. However, as a consequence, forces and moments are simplified and decoupled with helicopter dynamics in such methods. It is necessary to mention that, in autorotation procedures, the main rotor’s ability to generate forces is highly restrained by the states of helicopters, especially by the main rotor’s rotating speed and the inflow state of the rotor [13]. The rotor speed is the crucial factor that determines whether a trajectory leads to a safe autorotation landing. Thus, the dynamics of autorotation is of great significance for trajectory planning and should be involved explicitly.
In this paper, we present a realtime feasible autorotation trajectory planning method using functional tensortrain (FT) based dynamic programming (DP) algorithms, which ensures that general autorotation dynamics is satisfied along the trajectory. For validation of the dynamic feasibility of the trajectories, we also present a trajectorytracking controller based on active disturbance rejection control (ADRC) to make a helicopter model track the trajectories. The validation of the trajectories using the controller is then implemented on a sixdegreeoffreedom, highfidelity, multibodybased helicopter simulation model.
The trajectory planning method using functional tensortrain (FT) based dynamic programming (DP) algorithms will be presented in Section 2. The ADRCbased trajectorytracking controller is described in Section 3. In Section 4, we describe the sixdegreeoffreedom, highfidelity, multibody helicopter simulation model using the Tsinghua Rotorcraft Utility Simulation Tool (TRUST). In Section 5, autorotation trajectories are demonstrated with various initial conditions, and simulation results for validation are demonstrated and discussed. Finally, conclusions are presented in Section 6.
2. Autorotation Trajectory Planning Using FTBased DP Algorithms
As mentioned in Section 1, traditional trajectory planning methods for autorotation have been dealt within computationally expensive offline algorithms. In 2017, the first realtime feasible trajectory planning method for autorotation was demonstrated in [12], which is based on differential flatness of the rigidbody dynamics. Such methods simplify forces as direct inputs, thus leaving the helicopter autorotation dynamics not included during trajectory planning procedures. In this section, we introduce a realtime feasible trajectory planning method, using functional tensortrain (FT) [14] based dynamic programming (DP) algorithms, which guarantees a strict satisfaction of helicopter dynamics along the trajectory.
2.1. FTBased DP Algorithms
Functional tensortrainbased (FTbased) dynamic programming (DP) algorithms are newly proposed algorithms for solving highdimensional stochastic optimal control (SOC) problems, which are mainly discountedcost infinitehorizon Markov decision process (MDP) problems. Here, we give a brief review of the FTbased DP algorithms, and details can be found in [15, 16].
Consider a system described by stochastic differential equations (SDE) as follows: where is the state vector, is the control input, is a vector of independent unit Wiener processes, and and are generally nonlinear functions. The cost functional is defined as where is a discount factor, is the indicator function of the state boundary, and and are stage cost function and terminal cost function. The SOC problem is to find a control within a specified set on the time interval , such that the cost is minimized.
Next, using the MCA method [17], continuous SOC problems are discretized into the discrete Markov decision processes (MDPs). The discretized problem is turned into searching for a value function that satisfies the following recursive equation: where is the optimal discretized value function and is the transition probability function.
Then, a discountedcost infinitehorizon MDP is formulated and can be solved by the FTbased DP algorithms, which are FTbased value iteration algorithm, FTbased policy iteration algorithm, and FTbased multilevel algorithm, respectively.
For the traditional discretestate Markov decision processes, computational requirements grow exponentially with dimensionality. For example, if a MDP has 8 dimensions, and each dimension has a discretization of 10 points, such a problem involves a search space of points. In order to mitigate such curse of dimensionality, FTbased DP algorithms use lowrankfunctions, namely functional tensortrain, to represent value functions. The basic idea of functiontrain (FT) is to make a continuous analogue [14] of the tensortrain decomposition [18]. To be specific, it is a continuous version of tensortrain crossapproximation (TTcrossapproximation) [19], with the formulation as follows: where is a dimensional multivariable function. is a set of univariate functions, which are also called cores: where are the FT ranks evaluated by a continuous version of TTrounding [14] and are univariate hat functions.
With methods described above, the exponentially growing computational complexity of for a typical dynamic programming problem is compressed to a polynomially growing complexity of where is the operations during value evaluation, is the step of operations within each Bellman equation step.
To summarize, this method formulates a SOC problem as a dynamic programming problem and implements a continuous tensor decomposition method to compress such a problem, resulting in significant improvements of computing speed and storage savings. Such method has been shown to be able to work in realtime [15].
2.2. Formulating Autorotation Trajectory Planning Problems within FTBased DP Algorithms
In this part, we describe a general method for solving autorotation trajectory planning problems in the framework of FTbased DP algorithms. First of all, we need to describe the dynamics of a helicopter in a SOC form, as equation (1) shows.
2.2.1. Helicopter Dynamic Model Formulation
The helicopter dynamic function in the SOC equations is given by a nonlinear threedegreeoffreedom rigidbody helicopter autorotation dynamic model. Note that in this paper, two different helicopter dynamic models are described. The first one is the nonlinear threedegreeoffreedom rigidbody helicopter dynamic model, which is used in this section as the system dynamics in equation (1). The second one is the sixdegreeoffreedom, highfidelity, multibodybased helicopter simulation model, which will be used as a validation model in simulation and will be described in Section 5. The threedegreeoffreedom dynamic model is chosen here for computation considerations, and such model is capable of predicting steady collective pitch manipulations and rotor power consumption [13].
States variables for the nonlinear threedegreeoffreedom rigidbody helicopter dynamic model are chosen as which include height , horizontal velocity , rate of descent , pitch angle , pitch rate , and rotating speed of the main rotor . The control variables are collective pitch and longitudinal cyclic pitch . The dynamic model equations are formulated as follows [20]: where and are rotor forces resolved in the body axes. Note that all the axes are defined in accordance with [13]. and are pitch moment and rotor torque, respectively. and are the drag of the fuselage and the lift of the horizontal tail resolved in the body axes. and are the angle of attack of the fuselage and the flightpath angle. and are vertical and horizontal distance of rotor hub to the center of gravity. is the horizontal distance of the horizon tail from the center of gravity. Details of forces generated by the fuselage and the horizontal tail can be seen in [20].
Next, the rotor performance during autorotation is given by the following procedures. The induced velocity is calculated by [9]: where is the empirical correction factor of nonuniform inflow and is the induced velocity at hover. is the induced velocity parameter which has included the influence of the vortexring state, given by [9]: where and are normalized velocities resolved in body axes:
The ground effect factor is expressed as [13] where and is the rotor height when the aircraft is on the ground. Note that we ignore the influence of horizontal speed on the ground effect, and such assumption is reasonable for a typical autorotation procedure.
The rotor thrust coefficient is expressed as [13] where is the blade liftcurve slope, is the rotor solidity ratio, is the blade twist, and is the blade tiploss factor.
The rotor drag force coefficient is made up of the rotor profile drag and the induced drag . The rotor pitch moment coefficient is . and are expressed as follows: where is the blade’s mean profile drag coefficient, is the flap offset, is the number of blades, and is the blade’s first moment about the flap hinge. is the coning angle of rotor disk, and is the first harmonic cyclic flap, given by where is the normalized inflow velocity, expressed as
The power coefficient required by the main rotor is equal to the rotor torque coefficient , which is made up of the rotor induced power, the rotor profile power, the rotor parasite power, and the rotor climb power. A general expression is given by [13] where is the fuselage equivalent parasite drag area.
2.2.2. Formulation of SOC Problem for Autorotation Trajectory Planning
Following the helicopter autorotation dynamic model described above, and for numerical considerations [5], the stochastic differential equations are obtained by normalizing equations of (8). The normalized states and inputs are as follows:
In addition, the Jacobian matrix of control inputs is calculated by numerical difference method for any state and input . The diffusion function is set to for each term.
Next, the cost functions are designed. The stage cost is expressed by
The terminal cost is given by setting an absorbing region [15], which encourages a safe landing, as shown below: where and denote weighting factors of stage cost and terminal cost, respectively. We add cost functions to control inputs in order to avoid fierce manipulations. is the average horizontal speed, estimated by empirical flight test results and numerical simulation results. A safe autorotation landing is defined by
The design of cost functions is shown to be important for such autorotation problems, and parameters are needed to be adjusted by numerical experiments.
2.2.3. Trajectory Generation Using FTBased Algorithms
After the problem is formulated in the form of a SOC problem in Section 2.2.2, solutions are then obtained using FTbased algorithms.
Although FTbased algorithms are proposed for stochastic optimal control problems, we slightly alter such algorithms to generate a trajectory, rather than to obtain control inputs directly. The reason is that even though the nonlinear dynamics of helicopter described in Section 2.2 gives nice results in terms of predicting steady collective pitch manipulations and rotor power consumption, such model is not able to make good predictions of the helicopter dynamic response [21]. Therefore, the control inputs generated by the SOC controller are not adopted, and we make use of the trajectory instead. As mentioned before, the results of rotor power and collective pitch along the trajectory can be regarded reasonable.
We employ the FTbased oneway multigrid algorithm proposed in [15, 22]. For each discretization level, steps of FTbased policy iterations are applied. The trajectory of normalized states is obtained by an integration procedure of given initial conditions to the controller:
In specific, a fourthorder order RungeKutta method is implemented. Thus, preliminary trajectories are obtained by making unnormalized by equation (18). Because the trajectories are generated by integrations based on the dynamic equations of the helicopter, they satisfy the specified helicopter dynamics by nature.
However, although such trajectories are generated by solving dynamic equations, the trajectories still need to be smoothed. There are two main reasons. The first reason is that for realtime realizability, the integration time step (for example, 0.05 s in our study) is not small enough for a controller to generate sufficient differential information. The second reason is that although we have cost functions on control inputs, the control values still result in discontinuous variations, making the trajectory not smooth enough for our controller. Thus, we apply a fast interpolation method, i.e., cubic Hermite interpolation, to obtain smooth trajectories.
3. Autorotation TrajectoryTracking Controller Based on ADRC
In order to show the practicability of trajectories generated by FTbased DP methods, we demonstrate a trajectorytracking controller based on active disturbance rejection control (ADRC) and use the controller to make a highfidelity helicopter model track the trajectories. Tracking of position and velocity is vital for a successful autorotation landing, and the dynamics during autorotation may be more complicated than the formulations described in Section 2.2. Besides, there also exist unexpected disturbances during real flights. Thus, we use active disturbance rejection control (ADRC) methods to design the tracking controller. Active disturbance rejection control (ADRC) is proposed by Han [23], and the ADRC controller is capable of estimating inner modelling errors or outer disturbances, thus making compensation accordingly.
In specific, we implement the ADRCbased trajectorytracking controller described in [24], and such controller has been successfully validated through flight tests [24, 25]. The structure of the controller is shown in Figure 1: where and are position vectors and velocity vectors of the trajectory and represents the reference velocities that needed to be tracked by the inner loop controller.
The order of the ADRC controller differs with different channels. For the autorotation application, we use a 3rdorder ADRC controller for forward, vertical, lateral channels, and a PI controller for yaw channel. In general, an ADRC controller is made up of a tracking differentiator (TD), an extended state observer (ESO), and a nonlinear state error feedback (NLSEF). Disturbances or modelling errors of the system are observed by the ESO, and then NLSEF is utilized to restrain them. A typical architecture of an ADRC controller is shown in Figure 2.
Where v is the input signal, b_{0} is the input gain factor, is the output of the ESO. For brevity, only the 3rdorder ADRC controller is presented here. The TD of a 3rdorder is formulated as [25] where is an optimal control function defined in [23] and r and are controller parameters. The 3rdorder ESO is formulated as where h is the integration step of the system, is the observer gain, and function is defined in [23]. The 3rdorder NLSEF is formulated as
Note that because in this paper we focus on the longitudinal performance of the controller, the yaw controller can be regarded as a yaw stabilizer.
The input signal of the ADRC is processed by the following procedure: where , , , and are positions and velocities along the trajectory, is the vertical factor, and is the input signal of the ADRC controller.
The inner loop controller is a PI controller with damping feedbacks, and such controller is proven to be working well with ADRC in trajectorytracking applications [25].
Following the above process, a trajectory controller is obtained. It is needed to mention that when the helicopter is near the ground, the controller is reset to make the helicopter descent slowly. Such techniques to ensure a stable landing is also reported by [12].
4. The Tsinghua Rotorcraft Utility Simulation Tools (TRUST)
This section briefly introduces the Tsinghua Rotorcraft Utility Simulation Tools (TRUST), with which we model the S58 helicopter and validate the trajectories using the ADRC controller. Tsinghua Rotorcraft Utility Simulation Tools (TRUST) stems from [26], which is based on the framework of multibody dynamics. The TRUST simulator is able to simulate many kinds of rotorcrafts, including regular helicopters, tandem helicopters, and compound helicopters, with articulated or flexible rotor hubs. Regular helicopter models consist of main rotor dynamics, tail rotor dynamics, rigidbody fuselage dynamics, horizontal tails, and vertical tails. Such models have proven to be of good agreement with various test data [26].
4.1. Model Description
Considering both fidelity and simulation feasibility, we adopt the following modelling methods.
4.1.1. Main Rotor
Main rotor is modelled as an articulated rotor with rigid blades of certain twist. Each blade is attached to the hub through pitch, lag, and flap hinges. The inflow is modelled as the threestate PittPeters dynamic inflow [27]. Blade lift and drag forces are calculated using the blade element theory, and liftcurve slope is obtained from a nonlinear tabulate of angle of attack and relative airflow speed.
In addition, because the main rotor may traverse from powered lift state into vortexring state (VRS) during autorotation, a modification [28] of vortexring state is applied to the threestate PittPeters model.
Due to the multibody dynamic framework, each blade’s motion is solved individually by the Newton/Euler equations. Thus, flap angles of blades are not assumed to be small, nor the blades motions are treated as a consequence of a periodical disk, which involves inevitable accuracy loss for blade motion responses [29].
4.1.2. Tail Rotor
Tail rotor is modelled similar to the main rotor, except that the rotor hub is rigid and no VRS correction is applied.
4.1.3. Fuselage
Fuselage is regarded as a rigid body, with linear aerodynamic lift and drag coefficients.
4.1.4. Horizontal and Vertical Tails
The aerodynamic forces of horizontal and vertical tails are calculated by flat plate models.
For the purpose of validating autorotation trajectory planning and tracking, the nonlinear helicopter model described above is capable of predicting dynamic responses of the helicopter [28, 30].
4.2. Dynamics of TRUST Model
As mentioned in Section 2, there are two different helicopter dynamic models used in our study. The nonlinear threedegreeoffreedom rigidbody helicopter dynamics is capable of calculating steady collective pitch manipulations, rotor forces, and power, but not sufficient for dynamic response predictions. Due to the fact that the rotor forces in such models are essentially zeroorder processes, dynamic responses of the rotor cannot be well described.
In the TRUST helicopter model, a more sophisticated dynamic inflow model is involved, and rotor blade motions are accurately described in the multibody framework. Blade lift and drag forces are from nonlinear tabulates of experiment data. Furthermore, helicopter dynamics is well known for coupling effects between different channels [31]. Hence, for validation purposes, a highorder, sixdegreeoffreedom, helicopter dynamic model should be considered.
5. Numerical Experiments
In the first part, we show the trajectory planning results of various autorotation initial conditions based on the Sikorsky S58. In the second part, to validate the dynamic feasibility of the trajectory planning results, we implement the ADRC controller and validate the trajectory planner by making a 6DOF nonlinear highfidelity S58 model track the trajectory.
The Sikorsky S58 is a singleengined helicopter equipped with a 4bladed articulated rotor with a radius of 8.535 m. The takeoff weight used in our research is 4500 kg. Detailed parameters can be found in [32, 33], and main parameters among which can be found in Table 1.

5.1. Autorotation Trajectory Planning Results
The parameters of the cost functions in equations (19) and (20) are set as follows:
Note that costs for control inputs are added to avoid extreme manipulations. The empirical gliding speed is set as . The definition of a successful autorotation landing is defined as
In addition, we also implement a start cost function before the value iterations. A start cost function is used to initialize the global value function before the iterations begin. The start cost function can be expressed as follows: where is set to 4.0 m in this case. We use the oneway method of points for each dimension and set the max rank approximation of the core functions to . Using the oneway method, the initial value function of is obtained from the solutions of . Convergence is plotted in Figure 3.
The left panel shows when , the value function converges to a certain value around 2000 quickly. The relative error of the value function between different iterations is plotted in the right panel. It is shown that for the discretization of , the relative error is small for a value of about , which is near the converging threshold of .
Various autorotation cases using the trajectory planner are tested, with different initial height and horizontal speed:
Case 1. .
Case 2. .
Case 3. .
Case 4. .
Trajectories are obtained following procedures described in Section 2. We apply a Hermite interpolation of 20 points sampled from the preliminary trajectories. Note that such number of points can be regarded as sufficient, considering there are 10 points in [11] and 16 points in [12]. Results are shown in Figure 4.
As Figure 4 shows, the helicopter enters a steady descent quickly, and begins to decrease at about 20 meters above the ground. Such behaviours agree generally with [9]. For high initial speed cases, the forward speed is decreased immediately after the engine fails, in order to get the helicopter prepared for the final landing. Such behaviours are in accordance with [5]. All the trajectories terminate within the landing constraints defined in (28). Thus, the trajectory planning method using FTbased DP algorithms is able to generate successful autorotation trajectories that satisfy the helicopter autorotation dynamics described in Section 2. Furthermore, comparisons are made with results obtained from Sequential Quadratic Programming (SQP) using the SNOPT software package [34]. In specific, besides the landing specifications assumed in (28), path constraints are required, mainly as . A typical result is demonstrated using initial conditions of Case 1 and is shown in Figure 5.
Both trajectories in Figure 5 lead to safe autorotation landing, considering landing constraints are all satisfied. Besides, the two trajectories show similarities in terms of the rate of descent. The main difference is that the trajectory using the FTbased DP method tends to maintain the rotating speed by entering the steady autorotation right after the engine failure, while the trajectory obtained by SQP method appears to be more agile since the rotor speed does not reach the lower limit of the path constraint. In this respect, the autorotation trajectory using the FTbased DP method is more conservative in terms of keeping a higher rotor speed.
Both the above trajectory planning results are obtained using one core of a 3.20 GHz Inter i78700 CPU. Time costs for each trajectory generation using the FTbased DP method (shown as Algorithm 1) and SQP method (shown as Algorithm 2) are listed in Table 2. The average time cost of the FTbased DP method is 0.60 s, which indicates that the method is realtime feasible.

Before entering autorotation, the time delay of disengaging the clutch is a key factor that determines whether the subsequent maneuvers lead to a safe landing. A long time delay makes autorotation a tough task because of a low rotor speed. Here, the influence of time delay is considered by setting lowered rotating speeds. Initial conditions are given by , where for each trajectory.
From Figure 6, we can see that for lowered rotor speeds caused by different time delays, the trajectory planner is still capable of generating trajectories without modifying any cost function parameter. The rotating speed is recovered along the trajectory by certain maneuvers.
As for realtime applications, because the trajectory planner is able to generate a global trajectory by any specified initial condition, the computing speed shown in Table 2 can be considered realtime applicable with such hardware. The trajectory computing time cost in our study is about 0.6 seconds, which is similar to the time of a fast human pilot reaction. However, as mentioned above, the time delay is a key factor of the autorotation procedure, which should be shortened with best efforts. Therefore, although the trajectory planner is shown to be able to deal with situations of different time delays, an update frequency of at least 1.5 Hz is recommended for realtime applications, considering the time cost of other parts (such as detection of engine failure). Thus, we recommend running such algorithms with CPU of at least 3.1 GHz or higher.
5.2. TrajectoryTracking Simulation with the TRUST 6DOF Helicopter Dynamic Model
In this part, for validation of the dynamic feasibility of the trajectory, we demonstrate a sixdegreeoffreedom flight simulation, with the trajectory planner described in Section 2 and the trajectorytracking controller described in Section 3. For brevity, control parameters can be found in the Appendix. As mentioned before, the simulation is based on a sixdegreeoffreedom nonlinear multibodybased helicopter model of S58, with a VRS correction of the PittPeters dynamic inflow and an accurate dynamic description of the blade’s flap motions.
We demonstrate the simulation results of initial conditions from Case 1. The initial conditions of states are given by
We assume the engine failure starts at in the forward flight, and the simulations are terminated once the height descends to zero.
Simulation results are shown in Figure 7. Remember that is positive downward in the body axes and is positive in the noseup direction. The red lines denote the reference values that are needed to be tracked. The black lines denote the 6DOF nonlinear helicopter model responses. The blue line indicates that the helicopter is close to the ground, and the controller is reset to make a slow landing. The value of the reset height is chosen as 1.0 m by simulation experiments.
As the results shown in Figure 7, the velocities and attitude angles of the helicopter are , , , , and when landing, which satisfy the safe autorotation landing specifications defined above. The height history and rotating speed of the main rotor are shown in lower right and the upper right panel of Figure 7, which indicate that the rotating speed is still beyond 75%, although the ground effect is neglected in the sixdegreeoffreedom nonlinear helicopter model. Notice that the performance of horizontal speed is not as good as the vertical channel. However, we do not set horizontal distance as any kind of target variable in this study. We also find that increasing horizontal controlling gains affects the vertical channel instead. Apparently, a more sophisticated modelbased controller should improve the tracking performance, but our purpose of validating the trajectories is achieved. Simulation results show that following the trajectories generated by the trajectory planner, the helicopter is successfully guided to a safe landing. Thus the trajectory is both realtime feasible and dynamic feasible.
6. Conclusions
This paper demonstrates a trajectory planning method for autorotation, which is realtime feasible and guarantees a strict satisfaction of specified helicopter dynamics along the trajectory. A trajectorytracking controller is demonstrated to ensure the helicopter fly along the trajectories. Successful autorotation trajectories with various initial conditions are shown and discussed, and the time cost of trajectory generation procedures is around 0.60 s. A comprehensive validation is made based on a sixdegreeoffreedom highfidelity nonlinear S58 model using the TRUST simulator. Simulation results show that although the horizontal velocities are not tracked as well as the vertical channel, the trajectories generated by the trajectory planner and the trajectory tracker are capable of guiding a helicopter into a successful autorotation. Thus, the trajectory planner is realtime feasible, and the generated trajectories are dynamic feasible. Further potential improvements for the trajectory planning method could be selecting the derivatives of the control as control inputs of the threedegreeoffreedom helicopter dynamics, extending the twodimensional helicopter dynamics to a threedimensional helicopter dynamic model, and automatic tuning of the weight functions.
Appendix
A. TrajectoryTracking Controller Parameters
The controller parameters described in Section 3 are listed below. Parameters of the NLSEF are shown in Table 3.

Parameters of the ESO and TD are listed in Table 4.

Parameters of the inner loop are listed in Table 5.

Data Availability
The data used to support the findings of this study are available from the corresponding author upon request.
Conflicts of Interest
The authors declare that they have no conflicts of interest.
Acknowledgments
This work has been supported by the Zhuhai Longhua Helicopter Science and Technology Co., Ltd.
References
 B. L. Aponso, E. N. Bachelder, and D. Lee, “Automated autorotation for unmanned rotorcraft recovery,” in AHS International Specialists’ Meeting  Unmanned Rotorcraft: Design, Control and Testing, Proceedings, pp. 21–33, 2005. View at: Google Scholar
 S. Taamallah, “A qualitative introduction to the vortexringstate, autorotation, and optimal autorotation,” in 36th European Rotorcraft Forum, ERF 2010, pp. 464–492, Paris, France, 2010. View at: Google Scholar
 P. Bibik and J. Narkiewicz, “Helicopter optimal control after power failure using comprehensive dynamic model,” Journal of Guidance, Control, and Dynamics, vol. 35, no. 4, pp. 1354–1362, 2012. View at: Publisher Site  Google Scholar
 M. W. Floros, “Descent analysis for rotorcraft survivability with power loss,” in American Helicopter Society 65th Annual Forum, Grapevine, TX, USA, May 2009. View at: Google Scholar
 A. Y. Lee, A. E. Bryson Jr., and W. S. Hindson, “Optimal landing of a helicopter in autorotation,” Journal of Guidance, Control, and Dynamics, vol. 11, no. 1, pp. 7–12, 1988. View at: Publisher Site  Google Scholar
 E. N. Bachelder and B. L. Aponso, “An autorotation flight director for helicopter training,” Annual Forum ProceedingsAmerican Helicopter Society, vol. 59, no. 2, American Helicopter Society, Inc, 2003. View at: Google Scholar
 Y. Okuno and K. Kawachi, “Optimal control of helicopters following power failure,” Journal of Guidance, Control, and Dynamics, vol. 17, no. 1, pp. 181–186, 1994. View at: Publisher Site  Google Scholar
 C. L. Bottasso, A. Croce, D. Leonello, and L. Riviello, “Rotorcraft trajectory optimization with realizability considerations,” Journal of Aerospace Engineering, vol. 18, no. 3, pp. 146–155, 2005. View at: Publisher Site  Google Scholar
 W. Johnson, “Helicopter optimal descent and landing after power loss,” Tech. Rep., NASA TM 73244, Ames Research Center, National Aeronautics and Space Administration, Moffett Field, CA, USA, 1977. View at: Google Scholar
 E. B. Carlson, “An analytical methodology for category a performance prediction and extrapolation,” Annual Forum ProceedingsAmerican Helicopter Society, vol. 57, no. 2, American Helicopter Society, Inc, 2001. View at: Google Scholar
 A. A. Jhemi, E. B. Carlson, Y. J. Zhao, and R. T. N. Chen, “Optimization of rotorcraft flight following engine failure,” Journal of the American Helicopter Society, vol. 49, no. 2, pp. 117–126, 2004. View at: Publisher Site  Google Scholar
 S. Taamallah, X. Bombois, and P. M. J. Van den Hof, “Trajectory planning and trajectory tracking for a smallscale helicopter in autorotation,” Control Engineering Practice, vol. 58, pp. 88–106, 2017. View at: Publisher Site  Google Scholar
 W. Johnson, Helicopter Theory, Courier Corporation, 2012.
 A. A. Gorodetsky, S. Karaman, and Y. M. Marzouk, “Functiontrain: a continuous analogue of the tensortrain decomposition,” 2015, http://arxiv.org/abs/1510.09088. View at: Google Scholar
 A. Gorodetsky, S. Karaman, and Y. Marzouk, “Highdimensional stochastic optimal control using continuous tensor decompositions,” The International Journal of Robotics Research, vol. 37, no. 23, pp. 340–377, 2018. View at: Publisher Site  Google Scholar
 A. A. Gorodetsky, Continuous LowRank Tensor Decompositions, with Applications to Stochastic Optimal Control and Data Assimilation, Diss. Massachusetts Institute of Technology, 2017.
 H. Kushner and P. G. Dupuis, Numerical Methods for Stochastic Control Problems in Continuous Time, vol. 24, Springer Science & Business Media, 2013.
 I. V. Oseledets, “Tensortrain decomposition,” SIAM Journal on Scientific Computing, vol. 33, no. 5, pp. 2295–2317, 2011. View at: Publisher Site  Google Scholar
 I. Oseledets and E. Tyrtyshnikov, “TTcross approximation for multidimensional arrays,” Linear Algebra and its Applications, vol. 432, no. 1, pp. 70–88, 2010. View at: Publisher Site  Google Scholar
 Y. Okuno, K. Kawachi, A. Azuma, and S. Saito, “Analytical prediction of heightvelocity diagram of a helicopter using optimal control theory,” Journal of Guidance, Control, and Dynamics, vol. 14, no. 2, pp. 453–459, 1991. View at: Publisher Site  Google Scholar
 U. T. P. Arnold, J. D. Keller, H. C. Curtiss, and G. Reichert, “The effect of inflow models on the predicted response of helicopters,” Journal of the American Helicopter Society, vol. 43, no. 1, pp. 25–36, 1998. View at: Publisher Site  Google Scholar
 C. H. E. E.S. Chow and J. N. Tsitsiklis, “An optimal oneway multigrid algorithm for discretetime stochastic control,” IEEE Transactions on Automatic Control, vol. 36, no. 8, pp. 898–914, 1991. View at: Publisher Site  Google Scholar
 J. Han, “From PID to active disturbance rejection control,” IEEE transactions on Industrial Electronics, vol. 56, no. 3, pp. 900–906, 2009. View at: Publisher Site  Google Scholar
 C. Jiang, H. W. Wang, J. K. Li, and L. J. Li, “Trajectorytracking hybrid controller based on ADRC and adaptive control for unmanned helicopters,” Chinese Journal of Engineering, vol. 11, 2017. View at: Google Scholar
 C. Jiang, A Model ReferenceActive Disturbance Rejection Hybrid Trajectory Tracking Controller for Unmanned Helicopter, Diss. Tsinghua University, 2017.
 S. Deng, C. Jiang, Y. Wang, and H. Wang, “Helicopter flight dynamics simulation with continuestime unsteady vortex latticefree wake and multibody dynamics,” The Proceedings of the 2018 AsiaPacific International Symposium on Aerospace Technology (APISAT 2018). APISAT 2018. Lecture Notes in Electrical Engineering, vol 459, X. Zhang, Ed., Springer, Singapore, 2018. View at: Publisher Site  Google Scholar
 D. A. Peters and N. HaQuang, “Dynamic inflow for practical applications,” Journal of the American Helicopter Society, vol. 33, no. 4, 1998. View at: Google Scholar
 D. A. Peters and C. He, “Technical note: modification of massflow parameter to allow smooth transition between helicopter and windmill states,” Journal of the American Helicopter Society, vol. 51, no. 3, pp. 275–278, 2006. View at: Publisher Site  Google Scholar
 J. R. Majhi and R. Ganguli, “Modeling helicopter rotor blade flapping motion considering nonlinear aerodynamics,” Computer Modeling in Engineering and Sciences, vol. 27, no. 1/2, pp. 2526, 2008. View at: Google Scholar
 G. H. Gaonkar and D. A. Peters, “Effectiveness of current dynamicinflow models in hover and forward flight,” Journal of the American Helicopter Society, vol. 31, no. 2, pp. 47–57, 1986. View at: Publisher Site  Google Scholar
 G. D. Padfield, Helicopter Flight Dynamics, John Wiley & Sons, Chichester, UK, 2008.
 W. J. Hanley, G. DeVore, and S. Martin, An Evaluation of the Height Velocity Diagram of a Heavyweight, High Rotor Inertia, Single Engine Helicopter. No. FAAADS84, Federal Aviation Administration Washington Dc Aircraft Development Service, 1966.
 J. Scheiman, A tabulation of helicopter rotorblade differential pressures, stresses, and motions as measured in flight. Vol. 952, National Aeronautics and Space Administration, 1964.
 P. E. Gill, W. Murray, and M. A. Saunders, User’s Guide for SNOPT Version 7: Software for LargeScale Nonlinear Programming, Report NA 052, Department of Mathematics, University of California, San Diego, CA, USA, 2006.
Copyright
Copyright © 2020 Yunjie Wang et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.