Table of Contents Author Guidelines Submit a Manuscript
Mathematical Problems in Engineering
Volume 2015 (2015), Article ID 635395, 14 pages
Research Article

Trajectory Optimization for a Cruising Unmanned Aerial Vehicle Attacking a Target at Back Slope While Subjected to a Wind Gradient

1School of Mechatronical Engineering, Beijing Institute of Technology, Beijing 100081, China
2China Academy of Launch Vehicle Technology, Beijing 100076, China

Received 5 January 2015; Accepted 11 June 2015

Academic Editor: Andrzej Swierniak

Copyright © 2015 Tieying Jiang 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.


The trajectory of a tubular launched cruising unmanned aerial vehicle is optimized using the modified direct collocation method for attacking a target at back slope under a wind gradient. A mathematical model of the cruising unmanned aerial vehicle is established based on its operational and motion features under a wind gradient to optimize the trajectory. The motion characteristics of  “altitude adjustment” and “suicide attack” are taken into full account under the combat circumstance of back slope time key targets. By introducing a discrete time function, the trajectory optimization is converted into a nonlinear programming problem and the SNPOT software is applied to solve for the optimal trajectory of the missile under different wind loads. The simulation results show that, for optimized trajectories, the average attack time decreased by up to 29.1% and the energy consumption is reduced by up to 25.9% under specified wind gradient conditions. , , and have an influence on the flight trajectories of cruising unmanned aerial vehicle. This verifies that the application of modified direct collocation method is reasonable and feasible in an effort to achieve more efficient missile trajectories.

1. Introduction

Enemy armored vehicles, missile launchers, and support vehicles always use the favorable terrain to mask their location, such as a back slope. In modern warfare, however, it is not practical to cause effective damage in such situations, using conventional ammunition following its original trajectory. These types of targets can be effectively attacked through the best trajectory and attack strategy of a cruising unmanned aerial vehicle (CUAV) [1, 2]. Taking the tubular launched CUAV as an example, it will unfold and change to a conventional aircraft after launched by the tube for storage, transportation, and launch. After a short attitude adjustment, it will avoid mountains and obstacles rapidly and can attack time sensitive key targets at a back slope. A large number of studies on operational use and trajectory optimization of CUAVs have been conducted by scholars around the world. Reference [3] provided the optimal strategy for guiding a CUAV in reconnaissance and search during military operational procedures. In [4], the optimal ballistic trajectory problem of CUAV at attack stage with the adoption of genetic algorithm was solved. In [5], the operational efficiency of CUAV for attacking time sensitive key targets was studied. For trajectory planning, [6] determined the optimal path of small and medium-sized UAV attacking targets in urban warfare by application of TSP method. In [7], the optimal climbing trajectory for UAV through the study of UAV climbing performance was given. Reference [8] studied the shortest trajectory of small UAV spying on multiple targets through refueling. Reference [9] established a trajectory generation method of terrain following and threat avoidance for UAVs at low altitude. References [10, 11] presented the optimal strategy of multiple UAVs under cooperative search and attack. Considering the influence of wind field on UAV flight performance and aiming at studying the effects of wind field on trajectory, references [1216] studied the optimal trajectory of a small UAV when taking periodic motions in dynamic wind gradient field. In [17] the influences of wind power and solar energy on flight endurance of small UAV were studied. Reference [18] investigated the use of fuzzy logic in the flight trajectory of a UAV under the condition of known wind velocities. Reference [19] studied the influence of wind velocities between city buildings on the motion of a UAV. Reference [20] studied the optimal trajectory of a UAV collecting wind power through dynamic soaring. Considering that a CUAV combines both the characteristics of ammunition and a UAV, the features of unpowered gliding and posture adjustment as well as suicide attacks during the whole operation process should be fully included in the research of trajectory optimization.

In the present work, the motion features of portable tubular launched CUAV were first analyzed, and the kinematic and dynamic equations of CUAV dynamics at different stages under the background of attacking targets on a back slope were also analyzed according to ENU coordinate system. The model was discretized with a discrete time function according to constraint conditions such as initial position, terminal position, velocity, attitude, flight performance, and terrain and the solution was optimized by using software SNOPT. The optimal trajectory for attacking targets at back slope under different wind gradients and the methods in improving attack performance of CUAV by using wind gradient were studied.

2. Problem Description

Operating procedures of a portable tubular launched CUAV can be divided into three stages. The first stage is to provide enough initial kinetic energy with the aid of external input. The second is the form transformation stage, in which the CUAV enters into safety attitude with unpowered flight. The third stage is the mission flight.

After further analyses of the process, a diagram which takes time as the horizontal axis and CUAV features as the vertical axis was established. Based on the stage division, initial launching stage, unfolding stage, and mission stage were marked on the time axis, while the features of the CUAV were identified from motion features and form. To facilitate the systematic description of the process, we assumed that every feature is delineated instantaneously and the final analysis diagram is shown in Figure 1. According to the stage division, during the working process, the correspondence of time and CUAV form can be found. The CUAV keeps folding from the launch point to the unfolding point; after that the form transforms to conventional aircraft.

Figure 1: Operation process analysis of tubular launched CUAV.

Form change of tubular launched CUAV will directly affect its motion features. Interior ballistic motion is performed in changing from the tube from the launch stage to the unfolding stage. When it unfolds, the CUAV first enters into a safety attitude in unpowered flight, and then both the control system and the dynamic system begin to work with the power supply and the missile adjusts its attitude according to the route defined by the trajectory planning. Finally it arrives at the target area and attacks enemy target.

Combat background here is as follows: one of the combat groups considered is at the foot of the mountain. A high altitude UAV detects that certain enemy time sensitive key target is masking its location somewhere at a back slope and its location becomes clear. The UAV sends the specific coordinates to the combat group. The mission is for the combat group to attack the target at the back slope with a tubular launched CUAV, which is equipped with a one-time use, low-cost ammunition. The objective function of trajectory optimization for CUAV is set asin which is the flight time after CUAV unfolds.

3. Mathematical Model

3.1. Model Assumptions

(1)Magnitude and direction of the wind velocity are only related to attitude and are invariable over time.(2)Quality of CUAV is a constant and energy capacity of battery is enough to support the CUAV arriving in target area.(3)Launching angle , initial launching velocity and form transformation time of CUAV are fixed values as they are systematic design parameters.(4)CUAV instantly unfolds when it is ejected from launch tube. The transformation process of “loader” is too fast to cause any influence on attitude.

3.2. Kinematic and Dynamical Models

Particle kinematic and dynamical models of CUAV were established from unfolding according to the division of stages.

() The initialization of the form transformation stage is at (0, ), in which point the CUAV transitions to an unpowered flight stage. Particle kinematics and dynamics equations for this stage can be written asin which is the altitude, are (east, north) position, is the airspeed, is the heading angle measured clockwise from the north, is the acceleration of gravity, and and are lift and drag when cruising mission is unpowered flying. Definewhere and are the lift coefficient and the drag coefficient when cruising mission is unpowered flying, is air density, and is wing area.

(, , ) and (, , ) are wind speed and wind acceleration along the axis, respectively. Consider

() The mission flight time of CUAV is (, ). Its kinematics equations can be conveyed by (2) and its dynamics equations can be written as follows:in which is the thrust from airscrew and and are lift and drag between and . We havewhere and are lift coefficient and the drag coefficient between and . can be determined fromin which is the parasitic drag coefficient and and are constants defining a drag polar. Other variables are as listed above.

is ground velocity of CUAV and is defined asThe state variables of CUAVs in (2) and (6) are and the control variables are . According to the model assumptions, the position and posture of the CUAV at can be described as follows:where are the state variables.

3.3. Formal Modeling in Battlefield Environment

The battlefield environment includes factors such as the terrain, climate, and hydrology. According to the requirements of mission operation, the operational mountain discrete model and airspace gradient wind model were established.

3.3.1. Discrete Mountain Model

The three-dimensional terrain of mountain is shown in Figure 2. A cube area ABCD-EFGH was first built, which completely covered the mountain area, as shown in Figure 3. The area was divided into and sections along the -and -axes, respectively. These sections and three-dimensional figures cross each other and form point coordinates, so that the three-dimensional terrain is discretized into a grid. Space coordinates for each grid represent position and height of the point.

Figure 2: Three-dimensional mountain model.
Figure 3: Discretization of model.
3.3.2. Wind Velocity Model in Operation Airspace

Actual wind velocity changes according to a linear gradient, wind index or logarithm [21]. The wind speed is given bywhere is the maximum value of wind, is the maximum value of flight altitude, and is the type of wind gradient.

If , the profile of wind velocity changes as a straight line; if , it changes according to similar index; if , the profile changes according to the similar logarithm, as shown in Figure 4.

Figure 4: Wind profile changes.

To simplify (11), here we setThen the wind gradient model can be simplified asAssuming that is the wind direction angle between the wind direction and the -axis direction, then

3.4. Constraint Condition

() Initial position constraint conditions are where , , and denote the unfolding point coordinates when launching the CUAV; and refer to the velocity and launching angle of CUAV ejecting from launch tube, respectively.

() Terminal position constraint conditions arewhere denotes the position coordinates for back slope target; , , and are the allowable deviation in the position from the terminal attack point; and are velocity and attitude constraints to terminal attack of CUAV.

() Process constraint conditions areEquation (17) is the flight performance and mountain obstacle constraints, in which , , , , , , , , and represent position, velocity, flight path angle, propeller power, flight path drift angle, change rate of flight path angle, and change rate of flight path drift angle when the CUAV is at time , respectively. Also, , , , and are the limitation on flight performance. and are the limitation on the performance of maneuverability. In the following table, min and max, respectively, refer to the corresponding allowable minimum and maximum parameters; is the maximum power of battery; is the efficiency of battery power converting to propeller thrust power, which is specified as 0.6 [16]; represents the distance between two points, and and are coordinates of discrete points and safe flight distance, respectively.

() Control variable constraint conditions are

4. Modified Direct Collocation Method

Direct collocation method is one of main methods employed in trajectory optimization research. It generally divides the whole system course into equal segments and two endpoints of each segment which describes the changes of state variables over time with polynomial. It assumes that the change of control variables is linear [21] and the tubular launched CUAVs have the motion features of attitude adjustment and suicide attack. If we use the conventional average discrete method, the discrete model will be weakened or even the two motion features will be lost. Therefore, the time discrete function was established here for the CUAV and Simpson method was used on polynomial. The basic steps are listed as follows.

() Time course was discretized; namely,The corresponding state variables and control variables to time points can be, respectively, recorded as and . Discrete function of time is as follows:where and are discrete point distribution factor; is time discrete symmetry factor.

() Each state variable was expressed with cubic Hermit polynomial on its subinterval:, and , in which boundary conditions can be obtained as follows:and thenTo solve it, we can obtainAccording to above equations, the polynomial expressing the corresponding state can be solved based on the given state variables at the node of each group. The collocation point was defined as the midpoint of subinterval, , and nowwhere is (2) and (6) evaluated at .

() In order to improve the cubic polynomial fit to the change of state variables, the derivative at the collocation point on the subinterval should be equal to calculated from polynomial:The mathematical model of the CUAV was first discretized by adding a discrete time function and then the nonlinear programming problem was solved by using SNOPT in the Matlab software package with the established constraint conditions taken into account.

5. Simulation Calculation and Results Analysis

5.1. Simulation Parameter Setting

Initial position constraint isterminal position constraint isprocess constraint iscontrol variable constraint isthe time discrete function was set asOther parameters of CUAV and environmental parameters are shown in Table 1.

Table 1: Other parameters of loiter munition and environmental parameters.
5.2. Result and Analysis

The simulation was performed on a 2.6GHz PC and algorithm is real-time as computation time is 21s under Matlab environment. As can be observed from Figure 5, the generated trajectory is smooth, so that the CUAV can attack target avoiding mountain obstacles. From Figures 6~7, it can be seen that the optimized trajectory has obvious features of attitude adjustment and suicide attack after adding the discrete time function. When the CUAV unfolds, an attitude adjustment stage is clearly showed and the attitude is adjusted to maximum velocity 50m/s. When CUAV adjusts to cruising state, it climbs steadily and closes to the target around mountain obstacle. At the end of attack stage there is a clear dive-hike process and the flight speed is adjusted to terminal constraints velocity 30m/s. The trajectory without time discrete function loses the attitude adjustment information in the initial stage and the information reflecting attitude and control in suicide attack stage is less. The results indicate that discrete time function is effective and its related discretized model contains more motion features of CUAVs.

Figure 5: Operational flight trajectory of the CUAV ().
Figure 6: Influence of discrete time function on the trajectory.
Figure 7: Changing curve of state variable and control variable ().

The optimal trajectories of attack under wind gradient and without wind are quite different. As shown in Figures 8~9, the CUAV can avoid mountain obstacles and climb with smaller flight path drift angle at the beginning of attack trajectory. It climbs to maximum altitude 2000m at about 300s and then increases the flight path drift angle. At this point ground velocity of CUAV obviously increases under the influence of wind speed, up to nearly 75m/s as shown in Figures 10~11. The CUAV can avoid obstacles and keep at 2000m height. When it gets close to the target, it dives and attacks it.

Figure 8: CUAV flight trajectories (, ).
Figure 9: Changing curve of state variable and control variable variables (, ).
Figure 10: Inertial speed versus time.
Figure 11: Wind speed versus time.

In order to analyze the influence of the gradient in wind speed on attack time, here influence function is obtained as follows:where is the attack time for CUAV without wind gradient and is the attack time under wind gradient (Figure 16). Similarly the influence function of wind to energy consumption can be defined.

As shown in Table 2, when and , attack time of the CUAV attacking target can be reduced by 13.3% on average and energy consumption reduced by 7.72%.

Table 2: Influence of wind gradient on attack time and energy consumption.

The flight trajectories at different of wind direction angles and maximum values of wind are summarized in order to further analyze the influence of battlefield environment on the CUAV flight trajectory.

Figure 12 presents the flight trajectories of CUAV for different wind direction angles when and . Figure 13 shows the influence of wind direction angle on the attack time and energy consumption, and when , the attack time and energy consumption of CUAV can decrease 24.1% and 19.6%, respectively. The changes of state variable and control variable variables for different wind direction angles are shown in Figure 14.

Figure 12: CUAV flight trajectories of different wind direction angles.
Figure 13: Influence of on the attack time and energy consumption.
Figure 14: Changing curve of state variable and control variable variables at different .

Figure 15 presents the flight trajectories of CUAV for different maximum values of wind when and . Obviously, there have obvious difference between the trajectories with and without wind gradient. But the changes of have less influence on the flight trajectory when . The attack time and energy consumption of CUAV can decrease 29.1% and 25.9%, respectively. The changes of state variable and control variable variables for different maximum values of wind are shown in Figure 17.

Figure 15: CUAV flight trajectories for different maximum values of wind.
Figure 16: Influence of on attack time and energy consumption.
Figure 17: Changing curve of state variable and control variable variables for different maximum values of wind.

6. Conclusion

In the present work, the operational processes of a tubular launched CUAV were analyzed, and the mathematical model for describing the dynamics was established for a missile attacking time sensitive key targets at a back slope under wind loading. In view of the motion features of attitude adjustment and suicide attack for tubular launched CUAVs, a discrete time function was added in the direct collocation method. In this way, the problem was converted into a nonlinear programming problem and then the CUAV trajectory was optimized by using the software SNOPT. Simulation results indicate that this method can effectively keep the motion features of attitude adjustment and suicide attack. The simulation results show that, for optimized trajectories, the average attack time decreased by up to 29.1% and the energy consumption is reduced by up to 25.9% under specified wind gradient conditions, which has an important significance in enriching and perfecting the operational use of CUAV.

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.


  1. X.-L. Ji and G.-L. He, “Aerodynamic characteristics of gun-launched loitering munitions and its shape design,” Transaction of Beijing Institute of Technology, vol. 28, no. 11, pp. 953–961, 2008 (Chinese). View at Google Scholar · View at Scopus
  2. B. Li, J. Li, G. He, and D. Li, “Research on cooperative combat for integrated reconnaissance-attack-BDA of group LAVs,” Mathematical Problems in Engineering, vol. 2014, Article ID 123142, 6 pages, 2014. View at Publisher · View at Google Scholar · View at MathSciNet · View at Scopus
  3. M.-F. Guo, N.-J. Fan, and Z.-H. Yuan, “Battlefield operational strategy of loitering munition,” Acta Armamentarii, vol. 27, no. 5, pp. 944–947, 2006 (Chinese). View at Google Scholar · View at Scopus
  4. Y. Wang, D. Li, and Q. Shen, “Path Planning for the attack stage of a loitering unit,” Transactions of Beijing Institute of Technology, vol. 28, no. 1, pp. 8–10, 2008 (Chinese). View at Google Scholar
  5. X.-L. Ji, D.-L. Li, and G.-L. He, “Combat efficiency analyses for loitering missile attacking time-critical target,” Transaction of Beijing Institute of Technology, vol. 30, no. 1, pp. 46–49, 2010 (Chinese). View at Google Scholar · View at Scopus
  6. T. Shima, S. Rasmussen, and D. Gross, “Assigning micro UAVs to task tours in an urban terrain,” IEEE Transactions on Control Systems Technology, vol. 15, no. 4, pp. 601–612, 2007. View at Publisher · View at Google Scholar · View at Scopus
  7. P. Zhang, X. Wang, and X. Chen, “Climb trajectory optimization of UAV based on improved particle swarm optimization,” Computer Simulation, vol. 29, no. 4, pp. 92–94, 2012 (Chinese). View at Google Scholar
  8. K. Sundar and S. Rathinam, “Algorithms for routing an unmanned aerial vehicle in the presence of refueling depots,” IEEE Transactions on Automation Science and Engineering, vol. 11, no. 1, pp. 287–294, 2014. View at Publisher · View at Google Scholar · View at Scopus
  9. R. Zardashti, A. A. Nikkhah, and M. J. Yazdanpanah, “Constrained optimal terrain following/threat avoidance trajectory planning using network flow,” Aeronautical Journal, vol. 118, no. 1203, pp. 523–539, 2014. View at Google Scholar · View at Scopus
  10. H.-F. Guo, D.-L. Ding, W.-C. Wu, and Y.-L. Liu, “Long-range penetration and cooperative search decision-making of multiple UAVs,” Acta Armamentarii, vol. 35, no. 2, pp. 248–255, 2014 (Chinese). View at Publisher · View at Google Scholar · View at Scopus
  11. C.-L. Tang, C.-Q. Huang, H.-W. Du, H.-Q. Huang, D.-L. Ding, and C. Luo, “Study of trajectory planning for UCAV formation cooperative attack,” Acta Armamentarii, vol. 35, no. 4, pp. 523–530, 2014 (Chinese). View at Publisher · View at Google Scholar · View at Scopus
  12. Y. J. Zhao, “Optimal patterns of glider dynamic soaring,” Optimal Control Applications and Methods, vol. 25, no. 2, pp. 67–89, 2004. View at Publisher · View at Google Scholar · View at MathSciNet · View at Scopus
  13. Y. J. Zhao, “Taking advantage of wind energy in UAV operations,” in Proceedings of the InfoTech at Aerospace: Advancing Contemporary Aerospace Technologies and Their Integration, pp. 26–29, American Institute of Aeronautics and Astronautics, Arlington,Va, USA, 2005.
  14. Y. J. Zhao, “Extracting energy from downdraft to enhance endurance of uninhabited aerial vehicles,” Journal of Guidance, Control, and Dynamics, vol. 32, no. 4, pp. 1124–1133, 2009. View at Publisher · View at Google Scholar · View at Scopus
  15. H. Y. Singhania, Practical strategies of wind energy utilization for uninhabited aerial vehicles in loiter flights [Ph.D. thesis], The University of Minnesota, Minnesota, Minn, USA, 2008.
  16. W. Guo, Y. J. Zhao, and B. Capozzi, “Optimal unmanned aerial vehicle flights for seeability and endurance in winds,” Journal of Aircraft, vol. 48, no. 1, pp. 305–314, 2011. View at Publisher · View at Google Scholar · View at Scopus
  17. G. C. Bower, T. C. Flanzer, and I. M. Kroo, “Conceptual design of a small UAV for continuous flight over the ocean,” in Proceedings of the 11th AIAA Aviation Technology, Integration, and Operations Conference (ATIO '11), pp. 1–17, American Institute of Aeronautics and Astronautics, Virginia Beach, Va, USA, 2011.
  18. G. P. Kladisa, J. T. Economoua, K. Knowlesa, J. Lauberb, and T.-M. Guerra, “Energy conservation based fuzzy tracking for unmanned aerial vehicle missions under a priori known wind information,” Engineering Applications of Artificial Intelligence, vol. 24, no. 2, pp. 278–294, 2011. View at Publisher · View at Google Scholar · View at Scopus
  19. C. White, E. W. Lim, S. Watkins, A. Mohamed, and M. Thompson, “A feasibility study of micro air vehicles soaring tall buildings,” Journal of Wind Engineering and Industrial Aerodynamics, vol. 103, no. 1, pp. 41–49, 2012. View at Publisher · View at Google Scholar · View at Scopus
  20. V. Bonnin and C. C. Toomer, “Energy-harvesting mechanisms for UAV flight by dynamic soaring,” in Proceedings of the AIAA Atmospheric Flight Mechanics Conference (AFM '13), pp. 732–745, American Institute of Aeronautics and Astronautics, Boston, Mass, USA, 2013.
  21. C. R. Hargraves and S. W. Paris, “Direct trajectory optimization using nonlinear programming and collocation,” Journal of Guidance, Control, and Dynamics, vol. 10, no. 4, pp. 338–342, 1987. View at Google Scholar · View at Scopus