Table of Contents Author Guidelines Submit a Manuscript
Mathematical Problems in Engineering
Volume 2017, Article ID 3498350, 7 pages
Research Article

Optimal Skip Entry Trajectory for Lunar Return Vehicle with Trim-Flaps

1College of Electric and Information Engineering, Taiyuan University of Science and Technology, Taiyuan 312400, China
2College of Electric and Information Engineering, Jiangsu University of Technology, Changzhou 213001, China

Correspondence should be addressed to Yin Wang; moc.361@714wpx

Received 1 December 2016; Accepted 16 January 2017; Published 6 February 2017

Academic Editor: Stefan Balint

Copyright © 2017 Yin Wang and Keming Yao. 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.


A medium lift-to-drag ratio lunar return vehicle with trim-flaps is presented in this paper. The trajectory optimization design under heat-rate constrain for skip entry lunar return vehicle is analyzed. The optimization problem with a first-order state constraint is introduced. The trajectory applying the Pontryagin maximum principle under the performance of minimum heat is optimized, and the optimal expression of lift coefficient is derived. The simulation studies show that this research method can decrease the heat-rate effectively.

1. Introduction

When the lunar return vehicle reenters the atmosphere, its speed is up to 11 km/s, and its kinetic energy is about twice that of the near-earth orbiter. Low lift-drag ratio lunar return vehicle, such as the Apollo spacecraft, using the way of direct reentry has the large reentry overload and the high rate of heat flow in its return process, so it is not conducive to the safety return of return vehicle. And its poor mobility makes it difficult to achieve precise landing.

Appropriately increasing the return capsule lift-drag ratio and using the way of skip reentry can improve the safety performance of the return vehicle. It also can significantly improve the maneuvering ability of the lunar return vehicle and offer more choices for spacecraft landing site [13]. The vehicle with trim-flaps is a kind of medium lift-drag ratio vehicle, which has a simpler structure than that of lift type vehicle. Since the 1990s, scholars like Shui et al. [47] have done some studies in conceptual design and return optimal trajectory for vehicles with trim-flaps.

Skip entry is the way that return vehicle flies out of the atmosphere boosted by the lift after the vehicle reentry with a smaller reentry angle and then does some ballistic flight outside the atmosphere; then it reenters the atmosphere again [7, 8]. Istratie et al. have done a series of studies in the initial reentry optimal trajectory design of the near-earth orbit skip reentry [811].

In view of such fact that a lot of research has been done for reentry trajectory optimization design of near-earth orbit returning spacecraft, only initial reentry of vehicle skip reentry is discussed in this paper. A return vehicle with trim-flap is used as the lunar return vehicle model; the optimization problem with a first-order state constraint is introduced; the trajectory applying the Pontryagin maximum principle under the performance of minimum heat is optimized and the optimal expression of lift coefficient is derived in this paper.

2. Reentry Model of Sample Return Vehicle

The physical model of return vehicle with trim-flaps is shown in Figure 1.

Figure 1: Return vehicle with trim-flaps.

The maximum lift-drag ratio of vehicle is about 0.7. At the end of the vehicle there is the aerodynamic wing, and it is used to control the attitude of the vehicle in the reentry process. The lift-drag ratio curve of the vehicle is shown in Figure 2.

Figure 2: Curves of lift-to-drag ratio.

In the initial stage of vehicle design, the motion of the center of mass is only needed to be considered in order to understand the aerospace vehicle flight path and flight performance. In this paper, the earth is assumed to be round; only the motion equation of the longitudinal plane is considered. If the control system is in the ideal work, the effect of the attitude control system force can ignored, and the return vehicle is in unpowered reentry, so the equations of motion for vehicle can be described as follows [710]:where is vehicle speed, is vehicle flight path angle, is vehicle radial distance, is vehicle height, is vehicle air density, is vehicle acceleration of gravity, is vehicle quality, is vehicle reference area, is vehicle drag coefficient, and is vehicle lift coefficient.

Andwhere is sea level atmospheric density and is atmospheric scale height coefficient. For a certain aerospace vehicle, and are constants determined.

The input variable of this control system is vehicle lift coefficient . When the motion parameters of the vehicle reentry point are known and the control variable is determined, the motion equation has a unique solution.

3. The Description of Optimal Reentry Problem

In the skip reentry process, there are some effects on the vehicle by the aerodynamic force and the force of gravity, and the attitude is adjusted to change the flight path. After some time of flight, the vehicle escapes the atmosphere by lift, and then it enters the new Kepler orbit. After doing some ballistic movement, it will reenter the atmosphere again.

The attitude of the vehicle in the air flight is adjusted to change its lift-drag ratio. So the optimal flight trajectory can be found, and it can make the vehicle reach the minimum total heat absorption and meets the heat flow constraint.

During the flight, the heat flux equation of vehicle can be given as follows [9, 10]:Heat flow constraint is , and it can also be described as follows:where , is the maximum heat flow value, and it is also a positive real number.

Because does not explicitly include the control parameter and explicitly includes the control parameter on the first derivative of the time, the heat flow constraint is a first-order state variable inequality constraint.

With derivation of this number , we can obtain is defined asDuring the flight, the lift coefficient constraint iswhere is the maximum lift coefficient, and it is also a positive real number.

The optimization objective function is the relative value of the total heat absorption. It can be described as follows:The initial state of the vehicle is The terminal constraint of the vehicle is

The optimal control problem under the heat flow constraints can be described as follows: constraints (1), (4), and (7) should be met during the flight; constraint (9) should be met at the initial moment; constraint (10) should be met at the terminal time. Under the condition that all of the above constraints are met, the minimum problem of the performance index is to be researched.

4. Optimal Trajectory Design

According to Pontryagin maximum principle, Hamiltonian can be obtained as follows [1013]:where , , are adjoint variables.

The adjoint equations are described as follows:where is Kuhn-Tucker multiplier.

According to the literature [12], the optimal curve can be divided into the free section , the constraint section , and the free section . is the time to enter the restriction moment. is the time to exit the restriction moment.

In the free section, , , according to , we haveIn the restricted section, , must meet [14, 15]According to , we have Thus, the adjoint equations can be written as follows:At the moment , we havewhere is a real number.

From , we haveAccording to (13), (15), and (18), we have At the moment , we have

From the above analysis, the optimal reentry problem under the heat flow restriction is the two-point boundary value problem which must meet the connection requirements of the constraints section ends, when the initial state and the terminal state are known. Adjacent extreme method can be used to solve this problem. is selected as the initial moment for calculation.

At the moment , , , , , , and are unknown. According to , vehicle height at this time can be described as follows:According to , we have The value of can be calculated by adjacent extreme method. can be retroactively integrated from time to time to calculate the value . According to (8), the value of can be calculated. Substituting into (17), the values of , , and can also be calculated. Then, reverse operations of integral are done from , and the end time is determined by the initial state.

5. Simulation and Analysis

Calculations are done for two cases in this paper. One is in the case of restricted lunar return vehicle heat flow; the other is just the opposite. Vehicle parameters are given as follows: = 600 kg, = 0.754 m2,,,.

The initial values are given as follows: = 11 km/s,, = 110 km.

A large number of numerical simulations show that the aerospace vehicle has minimum value in the initial reentry. Heat flow constraint value can not take arbitrarily small value. In this paper, heat flow constraint value is given as follows: = 6.5 × 104 W/cm2.

For no power return vehicle, its terminal velocity should be less than 7800 m/s. In order to make the safe reentry, the secondary reentry angle can not be too large. Thus, the end state is set by = 110 km, = 7750 m/s,.

Through the simulation of optimization model, the actual terminal state is set by = 110000.2 m, = 7753.4 m/s,, = 141.2 s, = 3.88 s.

Curves of vehicle height, velocity, flight path angle, lift coefficient, lift-to-drag ratio, heat flux, pressure, and overload in the atmospheric flight changing with time in the two cases are given by Figures 310. The solid line is the heat constraint optimization curve, and the dashed line is no heat constraint optimization curve in the figures.

Figure 3: Time histories of the altitude.
Figure 4: Time histories of the velocity.
Figure 5: Time histories of the flight path angle.
Figure 6: Time histories of the lift coefficient.
Figure 7: Time histories of the lift-drag.
Figure 8: Time histories of the heat-rate.
Figure 9: Time histories of the dynamic pressure.
Figure 10: Time histories of the overload.

During the initial stage with heat constraint, lift coefficient can maintain maximum. This can decrease the speed of the vehicle faster than that in unconstrained case. So it can meet the heat flux constraints during the flight. The maximum dynamic pressure value is smaller than that in unconstrained case, and the maximum dynamic pressure value is greater than that in unconstrained case. The lift coefficient in keeping the maximum value will become smaller after a period of time. This will make the lift-drag ratio become larger. As the reentry vehicle is at low altitude, where the air density is larger, increasing the lift-drag ratio can slow down the speed of the vehicle descent rate. So it can meet the terminal velocity requirement. The flight time in the atmosphere of the return vehicle in heat flow constraint conditions is slightly longer than that in unconstrained conditions.

6. Conclusion

A medium lift-to-drag ratio lunar return vehicle with trim-flaps is presented as the programming of the lunar return vehicle in this paper. By the numerical simulation above, the following conclusions can be obtained:(1)In heat flow constraint conditions, the maximum dynamic pressure value of return vehicle is greater than that in unconstrained case. For lunar return vehicle with trim-flaps, the heat flux constraints can not only reduce the quality of thermal protection system of the return vehicle, but also reduce the loads on the control wing.(2)In heat flow constraint conditions, the maximum overload of the return vehicle is 12.6 g, which is 10.6 g larger than that in unconstrained case. That is acceptable for the unmanned vehicle. But the maximum overload will need to be constrained for the manned vehicle.

Competing Interests

The authors declare that they have no competing interests.


This study is supported by the National Natural Science Found of China (Grant no. 91016017), Natural Science Foundation of Shanxi (Grant no. 2015021089), and Doctoral Research Found of Taiyuan University of Science and Technology (Grant no. 20132020).


  1. C. Bai, J. Guo, X. Xu, and G. He, “Robust Skip earth entry guidance for a low L/D spacecraft,” in Proceedings of the 12th IEEE International Conference on Mechatronics and Automation (ICMA '15), pp. 1653–1658, Beijing, China, August 2015. View at Publisher · View at Google Scholar · View at Scopus
  2. P. Lu, C. Brunner, S. Stachowiak, G. F. Mendeck, M. Tigges, and C. Cerimele, “Verification of a fully numerical entry guidance algorithm,” in Proceedings of the AIAA Guidance, Navigation, and Control Conference, pp. 1969–2017, 2016.
  3. F. B. Wang and C. H. Dong, “Reentry trajectory optimization design for lunar return through coevolutionary algorithm,” Journal of Beijing University of Aeronautics and Astronautics, vol. 40, no. 5, pp. 629–634, 2014. View at Publisher · View at Google Scholar · View at Scopus
  4. Z. Shui, J. Zhou, and Z. Ge, “Lateral entry guidance for lunar return vehicles,” Modern Applied Science, vol. 6, no. 3, pp. 100–106, 2012. View at Publisher · View at Google Scholar · View at Scopus
  5. K. Mohan and M. Patterson, “Optimal trajectory and control generation for landing of multiple aircraft in the presence of obstacles,” in Proceedings of the AIAA Guidance, Navigation, and Control Conference, AIAA 2012-4826, Minneapolis, Minn, USA, August 2012.
  6. J. R. French, S. Alan Stern, M. Vozoff, T. McCallum, and C. Deiterich, “Architecture for lunar return using existing assets,” Journal of Spacecraft and Rockets, vol. 50, no. 4, pp. 838–847, 2013. View at Publisher · View at Google Scholar · View at Scopus
  7. J. Burkhardt, F. Zimmermann, and U. M. Schöttle, “Operational use of guided reentry capsules—system design solutions and mission safety considerations,” Aerospace Science and Technology, vol. 8, no. 7, pp. 635–644, 2004. View at Publisher · View at Google Scholar · View at Scopus
  8. V. Istratie, “Three-dimensional optimal skip entry with terminal maximum velocity,” in Proceedings of the 22nd Atmospheric Flight Mechanics Conference, pp. 19–26, AIAA, New Orleans, La, USA, August 1997. View at Scopus
  9. V. Istratie, “Optimal skip entry with terminal maximum velocity and heat constraint,” in Proceedings of the 7th AIAA/ASME Joint Thermophysics and Heat Transfer Conference, Albuquerque, NM, USA, June 1998.
  10. D. Ross and J. Russell, “Mars Ascent Vehicle (MAV): designing for high heritage and low risk,” in Proceedings of the 2012 IEEE Aerospace Conference, pp. 1–6, Big Sky, Mont, USA, March 2012. View at Publisher · View at Google Scholar
  11. V. Istratie, “Optimal skip entry into atmosphere with minimum heat and constrains,” in AIAA Atmosphere Flight Mechanics Conference, Denver, CO, Aug. 14–17, 2000, Collection of Technical Papers (A00-39676 10-08), pp. 288–294, American Institute of Aeronautics and Astronautis, Reston, Va, USA, 2000. View at Google Scholar
  12. V. Istratie, “Optimal skip entry into atmosphere with minimum heat,” in Proceedings of the 2003 AIAA Atmospheric Flight Mechanics Conference and Exhibit, pp. 11–14, Austin, Tex, USA, August 2003.
  13. H. Seywald, “Variational solutions for the heat-rate-limited aeroassisted orbital transfer problem,” Journal of Guidance, Control, and Dynamics, vol. 19, no. 3, pp. 686–692, 1996. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at Scopus
  14. H. Zhang, J. Zhou, and D. Wu, “The optimal aeroassisted orbital transfer heat flow constraints,” Shanghai Spacecraft, no. 4, pp. 6–11, 1999. View at Google Scholar
  15. A. Cassioli, D. Izzo, D. Di Lorenzo, and M. Locatelli, “Global optimization approaches for optimal trajectory planning,” Modeling and Optimization in Space Engineering, vol. 73, pp. 111–140, 2013. View at Google Scholar