Next Article in Journal
Optimal Circle-to-Ellipse Orbit Transfer for Sun-Facing E-Sail
Next Article in Special Issue
On the Size of the Safety Area around the Launch Trajectory of a Rocket
Previous Article in Journal
Selection Rules for Resonant Longitudinal Injector-Coupling in Experimental Rocket Combustors
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Receding Horizon Trajectory Generation of Stratospheric Airship in Low-Altitude Return Phase

The School of Aeronautics and Astronautics, Shanghai Jiao Tong University, Shanghai 200240, China
*
Author to whom correspondence should be addressed.
Aerospace 2022, 9(11), 670; https://doi.org/10.3390/aerospace9110670
Submission received: 30 June 2022 / Revised: 6 October 2022 / Accepted: 26 October 2022 / Published: 29 October 2022
(This article belongs to the Special Issue Aircraft Trajectory Design and Optimization (Volume II))

Abstract

:
The contribution of this paper is the proposal of a new receding horizon trajectory generation method for stratospheric airships’ return phase. Since the energy consumption, wind field and path constraints are restrictions during the return phase of airships at low altitude, it is crucial to develop novel trajectory optimization methods to ensure that the airship returns to the landing site. In this article, optimization objects and conditions of the return trajectory were constructed, considering the energy consumption and wind field. Then, a modified interior point method (MIP) was used to transform the inequality constraint, which is a simple and adaptable method used to improve the solving efficiency on the basis of a modified multiple shooting method (MMS). In addition, an adaptive gradient descent regulator was improved to reduce the influence on the optimization result due to different selections of the initial search point, and the convergence was made faster and more stable. Finally, considering the performance and path constraints of the airship, the effectiveness of the scheme was verified by numerical experiments under wind effects and a comparison of different methods.

1. Introduction

Stratospheric airships can extend the duration of the flight time (over a month) at near-space altitude (18–22 km) based on buoyancy and propulsion systems. Based on the advantages of a long endurance, large carrying capacity and wide work coverage, airships as a stratospheric real-time monitor or space transmission platform have attracted wide attention [1,2,3,4,5,6], having important applications in the fields of earth observation, emergency rescue and civil communication, etc. [7,8,9,10,11].
Stratospheric airships rely on the lift generated by lighter-than-air gases such as helium to maintain being stationary in the air, and altitude changes can be realized through ballonet valves opening and closing. The near-space region is approximately a quasi-zero wind area, having a small wind force. However the wind is very strong at low altitude. Thus, the applicable propulsion system is required to be different in the two regions. Based on this, the wind field, prescribed flight range, low-altitude propulsion system and forbidden region become constraints for the process of airships landing. A main challenge in the airship landing process is the planning and optimization of the airship’s descent trajectory.
The airship descent process is mainly divided into two steps: early-stage descent and low-altitude return, respectively, shown in Figure 1. The first step is to release the lighter-than-air gas in ballonets, reducing the buoyancy of the airship, and to drop rapidly through the troposphere. The second step is to abandon the ballast, decreasing the descent speed of the airship at an altitude of 3 km. Then, the low-altitude propulsion system is activated to bypass the forbidden region and reach the landing site within the designated flight range. In this process, the goal of trajectory optimization is the low-altitude propulsion system energy consumption and flight time, and the limiting factors are path constraints and the wind field. However, very little research has been conducted in this area so far [12,13,14,15].
Zhang and Li [16] studied an energy optimization problem based on the position potential method to improve the duration time of airships, and conducted a case study to validate this strategy. A three-dimensional optimal trajectory design was addressed for a stratospheric airship platform, considering a real jet stream and flight constraints [17]. However, path constraints such as forbidden zones were not addressed. Zhu and Li [18] aimed to improve the solar energy system of a near-space airship by optimizing the airship yaw angle. The recovery trajectory optimization of the stratospheric airship was discussed to find the economical recovery flight trajectory to the station-keeping site [19]. Lanteigne and Blouin [20] transformed the optimal trajectory problem into an optimal control problem and solved it using a pseudo-spectral method, providing a reference for the numerical optimization of trajectory. Although these methods [21,22,23,24,25] can also obtain airships flight trajectories, there seem to be some limitations.
  • As described above, the influence of the wind field is mainly considered in the stratosphere, and trajectory generation in the low-altitude region is not covered.
  • In previous studies, the limiting factor is mainly considered as the wind field, while the path constraint is not mentioned.
  • The airship is an aircraft with large inertia, and its state cannot be changed rapidly. The trajectory generated by many studies is not smooth; for example, the speed and attitude change rapidly in a short time, which will challenge the propeller and structure of the airship.
This paper takes the low-altitude return phase as the research background, expresses the limiting factors, and proposes a modified simple-structure and efficient numerical optimization method. The main contribution of this paper can be summarized as follows.
  • The vertical wind shear model and constant wind model at different altitudes were established according to the wind field data results. The landing site, wind field, forbidden region and flight range were transformed into a terminal constraint, penalty function and path constraint, respectively.
  • The trajectory generation of a stratospheric airship was solved by converting the boundary value problem into parameter optimization according to the modified multiple shooting method, and transforming the inequality constraint into a penalty function by using the modified interior point method.
  • An adaptive gradient descent regulator was used to reduce the influence on the optimization result due to different selections of the initial search point, and the convergence was made faster and more stable.
The section structure of this paper is as follows. The dynamic model and kinematic model of stratosphere airships are established in Section 2. The wind field model is also denoted. Then, numerical solution methods are proposed in Section 3. In Section 4, two numerical optimization results are carried out, which are the influence of wind effects on trajectory generation and the comparison between different methods.

2. Dynamic Model and Problem Formulation

The six-degrees-of-freedom (6-DOF) dynamics and kinematics equations of a stratospheric airship are introduced in this section [26]. The reference frames of an airship, including the earth reference frame (ERF), the body reference frame (BRF) and the path parallel frame (PPF), are defined in Figure 2, based on Zhou’s work [27]. P , Ω denote the relative position and the attitude of airships, respectively, and V , W v , w are the velocity and the angular velocity in the body reference frame of airships, respectively.
Define
ζ = P Ω , η = V W
P Ω V W T = x g ϕ u p y g θ v q z g ψ w r

2.1. Kinematics and Dynamics Equations of the Stratospheric Airship

Consider the kinematics of airships in the following form:
P ˙ = R I B V
Ω ˙ = Φ W
where R I B , Φ are the direction cosine transformation matrix and the angular transformation matrix from BRF to ERF, respectively. The velocity and the angular velocity can be converted from BRF to ERF. Both of them can be considered in the following Equations (5) and (6).
R I B = cos ψ cos θ cos ψ sin θ sin ϕ sin ψ cos ϕ cos ψ sin θ cos ϕ + sin ψ sin ϕ sin ψ cos θ sin ψ sin θ sin ϕ + cos ψ cos ϕ sin ψ sin θ cos ϕ cos ψ sin ϕ sin θ cos θ sin ϕ cos θ cos ϕ
Φ = 1 sin ϕ tan θ cos ϕ tan θ 0 cos ϕ sin ϕ 0 sin ϕ sec θ cos ϕ sec θ
The dynamic model can be denoted as
M · V W = F GB + F I + F A + F T
where M represents the inertia matrix, denoted by Equation (8), where m is the total mass of the airship, m i i , i = 1 6 represents the additional mass, I x , I y , I z , I x y , I x z and I y z stand for the inertia and x G , y G , z G represents the center of airship gravity.
M = m + m 11 0 0 0 m z G m y G 0 m + m 22 0 m z G 0 m 26 + m x G 0 0 m + m 33 m y G m 35 m x G 0 0 m z G m y G I x + m 44 I x y I x z m z G 0 m 53 m x G I x y I y + m 55 I y z m y G m 62 + m x G 0 I x z I y z I z + m 66
where F GB is the resultant force of gravity and buoyancy, F A , F I and F T are the aerodynamic force, inertial force, Coriolis force and the propeller force, respectively. The detailed expressions of these matrics mentioned above are given in [28].
Remark 1.
From the system input as the propeller force F T , it is clear that the stratospheric airship is under-actuated and lateral motion cannot be controlled directly.
Assumption 1.
It is known from [29] that the rolling angle ϕ of the airship is uncontrollable, which is set to zero in three-dimensional trajectory optimization.

2.2. Wind Field Models

Based on Azinheira’s work [30] and Sarma’s work [31], it is known that the horizontal wind movement is considered to be the steady wind with the disturbance of gust, which can affect the accuracy of the trajectory during the airship return phase. The amplitude of the wind varying with altitude is called the vertical wind shear, affecting the flight safety.
The wind profile generated by HWM14 [32] was used in the position of 85 E , 42 N , which is the outdoor experiment site. Horizontal wind velocities in different seasons are shown in Figure 3. The horizontal wind in summer and autumn is obviously weaker than that in winter and spring, which is also the reason why it is better to select the outdoor experiment in the summer. At an altitude of approximately 19.5 km, there will be an area with nearly zero wind velocity, which is called the quasi-zero wind layer, and is the frequently selected stationary point.
The summer wind field was selected as the numerical experiment environment as shown in Figure 4. Equation (9) denotes the fitting relationship between wind velocity and altitude.
w E ( h ¯ ) = i = 0 7 p ( 8 i ) e h ¯ i , w N ( h ¯ ) = i = 0 7 p ( 8 i ) e h ¯ i
where w E and w N are the east wind speed and north wind speed, respectively. In order to keep the data at the 0 mean value, h ¯ = h μ d / σ d is used as the independent variable for fitting, and μ d = 15,000 m and σ d = 8660.4 m. p i e and p i n ( i = 1 8 ) are the fitting coefficients, respectively, and the values are summarized as follows:
p e = [ 0.971 , 9.714 , 35.940 , 53.170 , 4.215 , 56.070 , 20.270 , 19.610 ] T , p n = [ 0.193 , 2.959 , 9.812 , 4.780 , 6.499 , 2.615 , 1.835 , 0.470 ] T
The acceleration of wind is derived in Equation (11).
w ˙ = R s ˙ w i + R s w ˙ i = R s w ˙ i Ω · R s w i = R s w ˙ i Ω · w
where w i is the wind speed in ERF, R s is the transform matrix from ERF to BRF and Ω is the angular velocity of the airship in BRF. The force generated by the wind is expressed by Equation (12).
F w = M Bw w ˙ + Ω · M Bw w
where M Bw = m Bw 0 3 × 3 0 3 × 3 J B w = m B I 3 + M v 0 3 × 3 0 3 × 3 J B + J v , m B is the mass of the airship and J B is the moment of buoyancy generated by the airship. M v and J v are the added mass and added moment of the airship, respectively. The effect of vertical wind shear is represented by Equation (13) with altitude changing.
F w ( h ) = M B w R s w ˙ i ( h ) M B w Ω · w ( h ) + Ω · M B w w ( h ) = M B w R s w ˙ i ( h ) M B w Ω · R s 1 w i ( h ) + Ω · M B w R s 1 w i ( h )
Assumption 2.
Based on the steady horizontal wind, the angular and vertical velocity of the wind are not considered.

3. Optimal Numerical Solution

Under the description of the above models, the receding horizon trajectory optimization model was established, and an objective function was designed by using the terminal cost and process cost to characterize the time and energy consumption of the trajectory. In the case of that, the path constraint, differential constraint and boundary constraints are aimed at expressing the forbidden region and the dynamics of the airship, as well as the environment model and target point, respectively. In this section, the trajectory optimization problem is transformed into a nonlinear programming problem (NLP) by using a parametric method to generate a smooth and effective trajectory. A flowchart of this process is shown in Figure 5.

3.1. Conversion of Bolza Problem

The airship receding horizon trajectory optimization problem is solved by converting it to an initial value problem on each optimization period, which has the Bolza form as Equation (14).
min ζ ( · ) , η ( · ) J = t 0 t f P ( t , ζ ( t ) , η ( t ) , u t ; s ) d t + R t f , ζ t f , η t f s . t . ζ t 0 η t 0 T x ^ 0 = 0 l t , ζ ˙ ( t ) , η ˙ ( t ) , ζ ( t ) , η ( t ) , u t ; s = 0 , t t 0 , t f g ζ ( t ) , η ( t ) ; s 0 , t t 0 , t f h ζ t f , η t f ; s = 0
where ζ η T and x ^ 0 denote the state variables and the initial state value, respectively, and s is the combined parameter variable. P and R are the process energy and terminal time consuming object, respectively. l, g and h represent the Lipschitz continuous function, path constraint and terminal boundary constraint of the airship, respectively.
u i k ( t ) = D i F Ti k ( t ) s i n μ i k ( t ) F Ti k ( t ) c o s μ i k ( t )
The action horizon is divided into I k : = t k 1 , t k , ( k = 1 , 2 , , N ) . u i k ( t ) , F Ti k ( t ) and μ i k ( t ) are each control component, each thrust force and the force vector angles. D i is the input transformation in Equation (15), approximated as
u i k ( t ) ord = 1 M + 1 γ i , ord k ( t ) ord λ i , ord k
where γ i , ord k ( t ) ord , λ i , ord k , o r d are the basis function, linear combination coefficient and the order number of the basis function. Cubic spline approximation is selected as Equation (17).
u i k ( t ) = u i ˙ k ( t k 1 ) ( t k t ) 3 6 τ + u i ˙ k ( t k ) ( t t k 1 ) 3 6 τ + ( t k t ) τ ( u i k ( t k 1 ) u i ˙ k ( t k 1 ) τ 2 6 ) + ( t t k 1 ) τ ( u i k ( t k ) u i ˙ k ( t k ) τ 2 6 ) , k = 1 , 2 , , N
where τ = t k t k 1 is the interval time. Thus, the control component is parameterized as σ i k = u i ˙ k ( t k 1 ) u i ˙ k ( t k ) u i k ( t k 1 ) u i k ( t k ) T .
Remark 2.
Due to the the airship complex dynamic characteristics and non-directivity of the objective function, it is unrealistic to solve the problem by using the analytic method [33,34], namely, the method of finding extremum by a functional. Refs. [35,36,37] used piecewise constants or polynomials to discretize control variables and was not suitable for the smoothing input on the optimized object. In this paper, the start and end states of each time interval are substituted into a cubic spline approximation formula. Robust and gentle trajectories can be generated more quickly using this method.
The MMS method approximately transforms the boundary value problem (BVP) into an NLP problem by means of discretization as Equation (18).
x ˙ k ( t ) = l σ k , x k ( t ) , t , t t k , t k + 1 , x k t k = χ k
where x ( t ) = ζ ( t ) ˙ η ( t ) ˙ ζ ( t ) η ( t ) is the state variables, and χ k = x k ( t 0 ) , ( k = 1 , 2 , , N ) is the initial value of the state vector in each sub-interval. In this way, the airship receding horizon trajectory optimization problem is transformed into a finite dimensional NLP problem to determine the control parameters σ k and auxiliary parameters χ k as Equation (19).
min χ k , σ k J = k = 1 N t k 1 t k P t , x k ( t ) , σ k d t + R x t f s . t . x ˙ ( t ) = k = 1 N l t , x k ( t ) , σ k x k t k = χ k x k t k + 1 χ k + 1 = 0 g t , x k ( t ) , σ k 0 h t , x k ( t ) , σ k = 0 u L σ k u U t k 1 t t k , k = 1 , 2 , , N
where u U and u L are the upper and lower boundary values of control parameters, respectively.
In order to ensure the continuity of the state curve, the MMS method adds the matching condition as the equality constraint.
Remark 3.
The main direct methods are the orthogonal collocation method (OC) and multiple shooting method (MS). Because conjugate variables are difficult to express explicitly, it is difficult to deal with the complex structure and repeatable deployment of airships. The key toward parameterization is that the new cubic spline is used to replace the Runge–Kutta segment approximating method to generate a smooth and efficient trajectory.

3.2. Parameter Nonlinear Optimization

The MIP method was adopted to deal with the large-scale nonlinear optimization problem transformed by the MMS method. The general form of the optimization problem is as follows.
min x R n J = f ( x ) / / s . t . g i ( x ) 0 , i = 1 , , n h j ( x ) = 0 , j = 1 , , m
Inequality constraints are transformed into a penalty object as Equation (21).
J = f ( x ) + i = 1 n B g i ( x ) s . t . h j ( x ) = 0
with g i ( x ) 0 , B g i ( x ) = 0 ; g i ( x ) > 0 , B g i ( x ) = . B g i ( x ) can be approximately replaced by a continuous function 1 μ ln g i ( x ) .
J = f ( x ) + 1 μ i = 1 n ln g i ( x ) s . t . h j ( x ) = 0
as μ tends to infinity, 1 μ ln g i ( x ) tends to B g i ( x ) . The corresponding Lagrangian function is Equation (23).
L ( x , u , ν ) = f ( x ) + i = 1 n u i g i ( x ) + j = 1 m ν j h j ( x ) .
Theorem 1.
Let the optimal value of the original problem be Y * , and the optimal value of the dual problem be R * . If Y * = R * , then there exists strong duality between the original problem and dual problem [38].
According to the strong duality of the primal problem and equivalent problem,
L ( x , u , ν ) = f ( x ) + i = 1 n u i g i ( x ) + j = 1 m ν j h j ( x ) = f ( x ) + i = 1 n 1 μ g i ( x ) g i ( x ) + j = 1 m ς j μ h j ( x ) = 0 .
Theorem 2.
If f ( x ) , g i ( x ) , h j ( x ) are differentiable, let x * be the optimal solution to the original optimization problem, and ( u i * , ν j * ) be the optimal solution of the dual problem. The strong duality is satisfied. Then, f x * + i = 1 n u i * g i x * + i = 1 m ν j * h j x * = 0 , where
g i x * 0 , i = 1 , , n h j x * = 0 , j = 1 , , m u i * 0 , i = 1 , , n u i * g i x * = 0 , i = 1 , , n .
The optimal solution of the original optimization problem is p * ; then,
p * L x * ( t ) , u * ( t ) , ν * ( t ) = f x * ( t ) + i = 1 n u i * ( t ) g i x * ( t ) + j = 1 m ν j * ( t ) h j x * ( t ) = f x * ( t ) + i = 1 n 1 μ g i x * ( t ) g i x * ( t ) + j = 1 m ν j * ( t ) h j x * ( t ) = f x * ( t ) n μ .
The new μ parameter is obtained from
μ i + 1 = min n ϵ tol , max κ μ μ i , μ i θ μ
Thus, when μ , f x * ( t ) p * is the optimal solution of the original problem.
The MIP method transforms the inequality constraint problem into an unconstrained optimization problem. Finally, the descending direction of the iteration point is obtained by the Newton iteration method, etc.
Remark 4.
The basic idea of the MIP method is to construct an interior penalty function so that the iteration point is always able to move in the feasible region constituted by various constraints [39,40]. However, it is a key point to choose the appropriate interior point as the initial point in the feasible region. In this paper, the iteration value of the barrier parameters μ can be adaptive-obtained.

4. Optimal Flight Trajectory Result

In this section, the summer wind field is considered as the environmental factor, because the horizontal wind speed in summer is relatively lower. The main parameters of airships [41] are able to be seen in Table 1. There are two numerical optimization solution analyses, which are, respectively, the analysis of the influence of the wind field on the optimization and the comparison of results between different methods.
Receding horizon optimization was used to divide the whole planning horizon into several moments at equal intervals, and the information of the future period and real-time detection was used to carry out online planning. Cost function optimization was carried out within each prediction cycle, and optimization results were input within the input cycle. Relevant parameters of the optimal numerical experiment are listed in the Table 2.

4.1. Optimization Effect under Wind Field

After ballonet valves are opened and the lighter-than-air gas is released, the low-altitude return starts to be executed. The flight path is set to return to the predetermined landing site from the stationary position, and the airship cannot pass through the forbidden regions with areas 15.90 km 2 , 7.07 km 2 and 12.57 km 2 during the flight (cylindrical areas shown in Figure 6). In the process of trajectory optimization, the upper limit of the linear velocity of the airship | u | is specified as 15 m / s and the upper limit of the descent speed | w | is specified as 3 m / s . The initial position is (0.48,0.6,3) km and the terminal position is (22,22,0) km.
Two cases were compared to reflect the influence of wind on the optimized trajectory. In the first case, wind effects were considered. In the second case, wind effects were not added. The three-dimensional optimization schematic comparison of two cases is expressed as Figure 6 and Figure 7.
The position and attitude changes during the airship’s return phase in the two cases are shown in Figure 8. Because the east wind direction is y axis forward, the north wind direction is x axis negative. Therefore, at position y, due to the influence of wind, the speed will be increased to reach the target landing site in advance, whereas, at position x, the speed will be reduced to delay the arrival of the target landing site. If there is no wind, the optimized trajectory is through the first and second forbidden region to reach the landing site to reduce mileage. However, with wind effects, such a trajectory is not feasible for airship dynamics. Thus, the airship will initially start at a yaw angle ψ of 62 deg to gradually reduce the yaw angle ψ to circumnavigate the forbidden region and eventually arrive at the landing site.
The comparison of horizontal velocity and the yaw angle velocity | r | during the airship’s return phase is shown in Figure 9 and Figure 10. The capacity of the airship angular velocity is of great significance for airship structural stability and stable operation, and represents the attitude change rate. The upper limit of the yaw angle velocity of the airship | ψ | is specified as 0.05 deg / s . With wind effects, the yaw angle velocity r will decrease and then increase reversely. When t = 2500 s , r will decrease again. Without the influence of the wind field, r will decrease the entire time until t = 3000 s , and r will increase in reverse.
The comparison of energy and propulsion system power is shown in Figure 11. As can be seen from the variation in kinetic energy and potential energy during the airship’s return phase, wind will increase the amplitude of the kinetic energy change, but the power accumulation of the propulsion system does not increase significantly because the airship is in a downwind state. In general, wind will affect the attitude change and planning strategy of airships, and the optimized trajectory needs to avoid the drastic fluctuation in attitude to ensure the flight safety of the airship.

4.2. Comparison between Results of Different Methods

The stratospheric airship is a large inertia aircraft with a complex structure, so its state cannot be changed quickly. Therefore, it is very important to select a numerical solution method that can deal with a complex dynamics system and generate an efficient and smooth trajectory. The trajectory optimization of complex systems was firstly transformed into an NLP problem by a parameterized method, in which, the direct methods of the MS method and OC method were mainly used.
The OC method does not need numerical integration and has a high solving efficiency. A small number of collocation points can achieve a high conversion efficiency. It is an extremely effective method with sequential quadratic programming (SQP) for aircraft trajectory optimization. However, its conjugate variable cannot be approximated explicitly, and a poor estimation of it has a great impact on mesh refinement and the sensitivity of optimization performance indexes. In Figure 12 and Figure 13, a comparison of optimization results of the MMS method with the MIP method, the MS method with the interior-point method (IP) and the OC method with the SQP method is shown.
Under the same receding horizon parameters, airship parameters, terminal conditions and path conditions, certain differences between the optimized trajectories appear. In the t = 2000 s , the yaw angle ψ under the OC method changed greatly, which brought a great challenge to the steering moment actuator of the airship. The reason for this is that it is difficult to estimate the conjugate variables reasonably.
The improvements to the modified method can be seen more clearly in Figure 14 and Figure 15. The modified method generates the best trajectory and reduces the mileage at the location. The speed of the airship can also be controlled to be smaller, with a reduction of 2 m / s for u and v. In addition, the yaw angle velocity r decreases by 0.02 deg / s . Figure 16 shows the variation curves of the airship’s energy and propulsion system power for different methods. The smoothness and stability of the energy change of the MMS+MIP method can be seen from the changing trend of potential energy and kinetic energy. The power integral of the airship’s propulsion system represents the main energy consumption during the return phase, and large power variations are unacceptable for the airship’s propulsion system. The power integration area of the MMS+MIP method is also significantly smaller than that of other methods. The advantages and disadvantages of the MS method and OC method are shown in Table 3.
In general, by modifying the connection link of the MS method and improving the convergence speed of the interior point adaptively, this method is beneficial to the smooth trajectory generation of an airship with large inertia.

5. Conclusions

This paper is concerned with the receding horizon trajectory optimization method of a stratospheric airship in a low-altitude return phase. The limit range of state variables and control variables was set according to the wind field, path constraints and energy consumption. The boundary value problem was transformed into a multi initial value problem using the MMS method, and inequality constraints were transformed into penalty functions using the MIP method. Essentially, the idea is to ensure that the airship is able to reach a prescribed geometric terminal target state with minimum dynamic model inputs smoothly. In addition, an adaptive gradient descent regulator was used to make the optimization solution fast and robust.
In the future, the energy acquisition of solar energy will be added as an optimization condition, and a trajectory tracker will be designed for the generated trajectory. The real-time sensing model of the wind field is also expected to make the trajectory generation more accurate.

Author Contributions

Conceptualization, Y.J. and J.T.; methodology, Y.W.; software, Y.J.; validation, Y.J., Y.W., J.T. and P.Z.; formal analysis, P.Z.; investigation, Y.J. and J.T.; resources, D.D.; data curation, Y.W.; writing—original draft preparation, Y.J. and J.T.; writing—review and editing, Y.W. and P.Z.; visualization, Y.J.; supervision, D.D.; project administration, P.Z.; funding acquisition, P.Z. and J.T. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded in part by the National Science Foundation of China (No. 62073216 and No. 51906141) and Shanghai Sailing Program under Grant (No. 19YF1421500).

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Tang, J.; Pu, S.; Yu, P.; Xie, W.; Li, Y.; Hu, B. Research on Trajectory Prediction of a High-Altitude Zero-Pressure Balloon System to Assist Rapid Recovery. Aerospace 2022, 9, 622. [Google Scholar] [CrossRef]
  2. Wang, Y.; Zhou, W.; Luo, J.; Yan, H.; Pu, H.; Peng, Y. Reliable intelligent path following control for a robotic airship against sensor faults. IEEE/ASME Trans. Mechatron. 2019, 24, 2572–2582. [Google Scholar] [CrossRef]
  3. Zhou, P.; Wang, Y.; Duan, D. Adaptive fault-tolerant directional control of an autonomous airship against actuator fault. In Proceedings of the 2016 35th Chinese Control Conference (CCC), IEEE, Chengdu, China, 27–29 July 2016; pp. 10538–10542. [Google Scholar]
  4. Tang, J.; Xie, W.; Wang, X.; Chen, C. Simulation and Analysis of Fluid–Solid–Thermal Unidirectional Coupling of Near-Space Airship. Aerospace 2022, 9, 439. [Google Scholar] [CrossRef]
  5. Wu, Y.; Wang, Q.; Duan, D.; Xie, W.; Wei, Y. Neuroadaptive output-feedback trajectory tracking control for a stratospheric airship with prescribed performance. Aeronaut. J. 2020, 124, 1568–1591. [Google Scholar] [CrossRef]
  6. Tang, J.; Duan, D.; Xie, W. Shape Exploration and Multidisciplinary Optimization Method of Semirigid Nearing Space Airships. J. Aircr. 2022, 59, 946–963. [Google Scholar] [CrossRef]
  7. Tang, J.; Xie, W.; Wang, X.; Chen, Y.; Wu, J. Study of the Mechanical Properties of Near-Space Airship Envelope Material Based on an Optimization Method. Aerospace 2022, 9, 655. [Google Scholar] [CrossRef]
  8. Teng, Y.; Song, Y.; Tong, S.; Zhang, M. Acquisition Performance of Laser Communication System Based on Airship Platform. Acta Opt. Sin. 2018, 38, 0606005. [Google Scholar] [CrossRef]
  9. Yu, Z.; Zhang, Y.; Jiang, B.; Su, C.Y.; Fu, J.; Jin, Y.; Chai, T. Distributed fractional-order intelligent adaptive fault-tolerant formation-containment control of two-layer networked unmanned airships for safe observation of a smart city. IEEE Trans. Cybern. 2021, 52, 9132–9144. [Google Scholar] [CrossRef]
  10. Zhao, Z.; Dong, X.; Feng, J.; Liang, X.; Hu, C. A Simulating Method of Airship-Borne Polarimetric Weather Radar for Typhoon Observation. In Proceedings of the IGARSS 2020—2020 IEEE International Geoscience and Remote Sensing Symposium, IEEE, Waikoloa, HI, USA, 19–24 July 2020; pp. 5306–5309. [Google Scholar]
  11. Tang, J.; Wang, X.; Duan, D.; Xie, W. Optimisation and analysis of efficiency for contra-rotating propellers for high-altitude airships. Aeronaut. J. 2019, 123, 706–726. [Google Scholar] [CrossRef]
  12. Chu, X.; Lin, X.; Lan, W. Energy-optimal rectilinear trajectory for stratospheric airships in constant wind field. In Proceedings of the 2013 IEEE International Conference on Cyber Technology in Automation, Control and Intelligent Systems, IEEE, Nanjing, China, 26–29 May 2013; pp. 264–269. [Google Scholar]
  13. Guo, X.; Zhu, M. Ascent trajectory optimization for stratospheric airship with thermal effects. Adv. Space Res. 2013, 52, 1097–1110. [Google Scholar] [CrossRef]
  14. Zhu, B.J.; Yang, X.X.; Deng, X.L.; Ma, Z.Y.; Hou, Z.X.; Jia, G.W. Trajectory optimization and control of stratospheric airship in cruising. Proc. Inst. Mech. Eng. Part I J. Syst. Control. Eng. 2019, 233, 1329–1339. [Google Scholar] [CrossRef]
  15. Mueller, J.; Zhao, Y.; Garrard, W. Sensitivity and solar power analysis of optimal trajectories for autonomous airships. In Proceedings of the AIAA Guidance, Navigation, and Control Conference, Chicago, IL, USA, 10–13 August 2009; p. 6014. [Google Scholar]
  16. Zhang, L.; Li, J.; Jiang, Y.; Du, H.; Zhu, W.; Lv, M. Stratospheric airship endurance strategy analysis based on energy optimization. Aerosp. Sci. Technol. 2020, 100, 105794. [Google Scholar] [CrossRef]
  17. Lee, S.; Bang, H. Three-dimensional ascent trajectory optimization for stratospheric airship platforms in the jet stream. J. Guid. Control. Dyn. 2007, 30, 1341–1351. [Google Scholar] [CrossRef]
  18. Zhu, W.; Li, J.; Xu, Y. Optimum attitude planning of near-space solar powered airship. Aerosp. Sci. Technol. 2019, 84, 291–305. [Google Scholar] [CrossRef]
  19. Wang, J.; Meng, X.; Li, C. Recovery trajectory optimization of the solar-powered stratospheric airship for the station-keeping mission. Acta Astronaut. 2021, 178, 159–177. [Google Scholar] [CrossRef]
  20. Blouin, C.; Lanteigne, E.; Gueaieb, W. Trajectory optimization of a small airship in a moving fluid. Trans. Can. Soc. Mech. Eng. 2016, 40, 191–200. [Google Scholar] [CrossRef] [Green Version]
  21. Das, T.; Mukherjee, R.; Cameron, J. Optimal trajectory planning for hot-air balloons in linear wind fields. J. Guid. Control. Dyn. 2003, 26, 416–424. [Google Scholar] [CrossRef] [Green Version]
  22. Ceruti, A.; Marzocca, P. Heuristic optimization of Bezier curves based trajectories for unconventional airships docking. Aircr. Eng. Aerosp. Technol. 2017, 89, 76–86. [Google Scholar] [CrossRef]
  23. Zhang, J.; Yang, X.; Deng, X.; Lin, H. Trajectory control method of stratospheric airships based on model predictive control in wind field. Proc. Inst. Mech. Eng. Part G J. Aerosp. Eng. 2019, 233, 418–425. [Google Scholar] [CrossRef]
  24. Bestaoui, Y.; Kahale, E. Time optimal 3D trajectories for a lighter than air robot with second order constraints with a piecewise constant acceleration. J. Aerosp. Inf. Syst. 2013, 10, 155–171. [Google Scholar] [CrossRef]
  25. Mueller, J.B.; Zhao, Y.J.; Garrard, W.L. Optimal ascent trajectories for stratospheric airships using wind energy. J. Guid. Control. Dyn. 2009, 32, 1232–1245. [Google Scholar] [CrossRef]
  26. Liu, S.; Sang, Y.; Jin, H. Robust model predictive control for stratospheric airships using LPV design. Control. Eng. Pract. 2018, 81, 231–243. [Google Scholar] [CrossRef]
  27. Zhou, W.X.; Xiao, C.; Zhou, P.F.; Duan, D.P. Spatial path following control of an autonomous underactuated airship. Int. J. Control. Autom. Syst. 2019, 17, 1726–1737. [Google Scholar] [CrossRef]
  28. Chen, L.; Zhou, G.; Yan, X.; Duan, D. Composite control of stratospheric airships with moving masses. J. Aircr. 2012, 49, 794–801. [Google Scholar] [CrossRef]
  29. Zheng, Z.; Zhu, M.; Shi, D.; Wu, Z. Hovering control for a stratospheric airship in unknown wind. In Proceedings of the AIAA Guidance, Navigation, and Control Conference, National Harbor, MA, USA, 13–17 January 2014; p. 0973. [Google Scholar]
  30. Azinheira, J.R.; de Paiva, E.C.; Bueno, S.S. Influence of wind speed on airship dynamics. J. Guid. Control. Dyn. 2002, 25, 1116–1124. [Google Scholar] [CrossRef]
  31. Sarma, A.; Hochstetler, R.; Wait, T. Optimization of airship routes for weather. In Proceedings of the 7th AIAA ATIO Conf, 2nd CEIAT Int’l Conf on Innov and Integr in Aero Sciences, 17th LTA Systems Tech Conf; followed by 2nd TEOS Forum, Belfast, Ireland, 18–20 September 2007; p. 7880. [Google Scholar]
  32. Gul, B.; Ameen, M.A.; Verhulst, T.G. Correlation between foEs and zonal winds over Rome, Okinawa and Townsville using Horizontal Wind Model (HWM14) during solar cycle 22. Adv. Space Res. 2021, 68, 4658–4664. [Google Scholar] [CrossRef]
  33. Athans, M.; Canon, M. On the fuel-optimal singular control of nonlinear second-order systems. IEEE Trans. Autom. Control 1964, 9, 360–370. [Google Scholar] [CrossRef]
  34. Bongiorno, J. Minimum-energy control of a second-order nonlinear system. IEEE Trans. Autom. Control 1967, 12, 249–255. [Google Scholar] [CrossRef]
  35. Diedam, H.; Sager, S. Global optimal control with the direct multiple shooting method. Optim. Control. Appl. Methods 2018, 39, 449–470. [Google Scholar] [CrossRef]
  36. Greco, C.; Di Carlo, M.; Vasile, M.; Epenoy, R. Direct multiple shooting transcription with polynomial algebra for optimal control problems under uncertainty. Acta Astronaut. 2020, 170, 224–234. [Google Scholar] [CrossRef]
  37. Chen, Y.; Scarabottolo, N.; Bruschetta, M.; Beghi, A. Efficient move blocking strategy for multiple shooting-based non-linear model predictive control. IET Control Theory Appl. 2020, 14, 343–351. [Google Scholar] [CrossRef] [Green Version]
  38. Wächter, A.; Biegler, L.T. On the implementation of an interior-point filter line-search algorithm for large-scale nonlinear programming. Math. Program. 2006, 106, 25–57. [Google Scholar] [CrossRef]
  39. Guo, C.; Li, D.; Zhang, G.; Ding, X.; Curtmola, R.; Borcea, C. Dynamic interior point method for vehicular traffic optimization. IEEE Trans. Veh. Technol. 2020, 69, 4855–4868. [Google Scholar] [CrossRef]
  40. Jiang, H.; Kathuria, T.; Lee, Y.T.; Padmanabhan, S.; Song, Z. A faster interior point method for semidefinite programming. In Proceedings of the 2020 IEEE 61st Annual Symposium on Foundations of Computer Science (FOCS), IEEE, Durham, NC, USA, 16–19 November 2020; pp. 910–918. [Google Scholar]
  41. Wang, X.L.; Fu, G.Y.; Duan, D.P.; Shan, X.X. Experimental investigations on aerodynamic characteristics of the ZHIYUAN-1 airship. J. Aircr. 2010, 47, 1463–1468. [Google Scholar] [CrossRef]
Figure 1. Flow chart of the airship descent process.
Figure 1. Flow chart of the airship descent process.
Aerospace 09 00670 g001
Figure 2. Reference frames of a stratospheric airship.
Figure 2. Reference frames of a stratospheric airship.
Aerospace 09 00670 g002
Figure 3. (left) North wind velocity. (right) East wind velocity.
Figure 3. (left) North wind velocity. (right) East wind velocity.
Aerospace 09 00670 g003
Figure 4. North and east component of the wind velocity in summer.
Figure 4. North and east component of the wind velocity in summer.
Aerospace 09 00670 g004
Figure 5. Flow chart of receding horizon trajectory generation during the airship return phase.
Figure 5. Flow chart of receding horizon trajectory generation during the airship return phase.
Aerospace 09 00670 g005
Figure 6. Three-dimensional schematic diagram of optimized trajectory.
Figure 6. Three-dimensional schematic diagram of optimized trajectory.
Aerospace 09 00670 g006
Figure 7. Top view of optimized trajectory.
Figure 7. Top view of optimized trajectory.
Aerospace 09 00670 g007
Figure 8. Position, attitude and change rate of the airship.
Figure 8. Position, attitude and change rate of the airship.
Aerospace 09 00670 g008
Figure 9. Horizontal speed and yaw angle velocity of the airship.
Figure 9. Horizontal speed and yaw angle velocity of the airship.
Aerospace 09 00670 g009
Figure 10. Control forces and torques of the airship.
Figure 10. Control forces and torques of the airship.
Aerospace 09 00670 g010
Figure 11. Energy and power changes during the airship’s return phase.
Figure 11. Energy and power changes during the airship’s return phase.
Aerospace 09 00670 g011
Figure 12. Three-dimensional schematic diagram of different methods.
Figure 12. Three-dimensional schematic diagram of different methods.
Aerospace 09 00670 g012
Figure 13. Top view of optimized trajectories between different methods.
Figure 13. Top view of optimized trajectories between different methods.
Aerospace 09 00670 g013
Figure 14. Position, attitude and change rate of the airship between different methods.
Figure 14. Position, attitude and change rate of the airship between different methods.
Aerospace 09 00670 g014
Figure 15. Horizontal speed and yaw angle velocity of the airship between different methods.
Figure 15. Horizontal speed and yaw angle velocity of the airship between different methods.
Aerospace 09 00670 g015
Figure 16. Energy and power changes between different methods.
Figure 16. Energy and power changes between different methods.
Aerospace 09 00670 g016
Table 1. Airship model parameters.
Table 1. Airship model parameters.
Airship ParametersValue
Length, m 25.0
Maximum diameter, m 7.576
Fineness ratio of the hull 3.3
Volume of the hull, m 3 750.0
Surface area, m 2 480.388
Location of maximum diameter, m 9.840
Moment center, m 2 12.001
Reference area, m 2 82.544
Reference length, m 25.0
Volume Reynolds number1.8–9.3 × 10 6
Table 2. Receding horizon parameters.
Table 2. Receding horizon parameters.
Receding Horizon ParametersValue
Sampling interval, s 0.5
Predict cycle600
Input interval, s 0.2
Action cycle500
Table 3. Advantages and disadvantages of MS method and OC method.
Table 3. Advantages and disadvantages of MS method and OC method.
MethodAdvantageDisadvantage
MS methodThe structure is simple and easy to implement, and the result is smooth and stableThe convergence domain is narrow and the time grid is evenly divided
OC methodHigh solving efficiency and high conversion efficiencyConjugate variable cannot be approximated explicitly; high sensitivity of mesh refinement and optimization index
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Jing, Y.; Wu, Y.; Tang, J.; Zhou, P.; Duan, D. Receding Horizon Trajectory Generation of Stratospheric Airship in Low-Altitude Return Phase. Aerospace 2022, 9, 670. https://doi.org/10.3390/aerospace9110670

AMA Style

Jing Y, Wu Y, Tang J, Zhou P, Duan D. Receding Horizon Trajectory Generation of Stratospheric Airship in Low-Altitude Return Phase. Aerospace. 2022; 9(11):670. https://doi.org/10.3390/aerospace9110670

Chicago/Turabian Style

Jing, Yuhao, Yang Wu, Jiwei Tang, Pingfang Zhou, and Dengping Duan. 2022. "Receding Horizon Trajectory Generation of Stratospheric Airship in Low-Altitude Return Phase" Aerospace 9, no. 11: 670. https://doi.org/10.3390/aerospace9110670

APA Style

Jing, Y., Wu, Y., Tang, J., Zhou, P., & Duan, D. (2022). Receding Horizon Trajectory Generation of Stratospheric Airship in Low-Altitude Return Phase. Aerospace, 9(11), 670. https://doi.org/10.3390/aerospace9110670

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop