High-resolution trajectory refinement via hodograph ballooning

Solutions to optimal control problems obtained via direct methods are usually of low resolution. As the desired resolution becomes higher, the number of parameters in the resulting parameter optimization problem increases nonlinearly. A concatenated approach to trajectory optimization previously introduced by the authors to obtain high resolution optimal solutions with low storage requirements is applicable only to a specific class of problems, such as minimum time and minimum energy optimal control problems. The hodograph ballooning approach introduced in this paper generalizes the concatenated approach to the most general optimal control problem. Introduction The development of powerful nonlinear programming software [1, 2] has made direct approaches to trajectory optimization the method of choice for most practical applications [3-5]. The discretization of an optimal control problem to a finite dimensional parameter optimization problem avoids many of the troublesome pathologies that are typical for the optimization over infinite dimensional function spaces [6-8]. Furthermore, direct approaches usually do not. require guessing the switching structure, and the dynamical system is devoid of the nonintuitive costate variables. However, to achieve good resolution it is usually necessary to refine some type of discretization mesh, and, ultimately, to use a large number of design parameters. For typical gradient search methods, this increase in the dimension of the optimization problem has nearly always a negative effect on the convergence rate and the robustness. Additionally, CPU time and RAM requirements for each iteration increase with the second power of the number of parameters (third power for non-sparse solvers) so that direct optimization approaches can become quite expensive. In [9], a method was introduced that achieves excellent resolution without solving high-dimensional nonlinear programming problems. This was achieved by lo'Senior Engineer, Member AIAA. t Supervising Engineer, Member AIAA. 117 Research Drive, Hampton, VA 23666. Copyright ©1996 by the American Institute of Aeronautics and Astronautics, Inc. All rights reserved cally optimizing the trajectory on an iterative sequence of short, partly overlapping subarcs. However, this approach is not applicable to the most general problem formulation. Loosely sneaking, it was required in [9] that the cost state does not appear explicitly in the right-hand side of the equations of motion. The present paper introduces a method through which the approach introduced in [9] becomes applicable to the most general problem formulation. Application of this method, however, requires a priori knowledge of the optimal cost. In practice it is observed that the optimal cost can be obtained with very high precision with reasonably low order discretizations. Once a good guess for the optimal cost is obtained presented in this paper can be applied to refine the trajectory, or, if the optimal cost value is not known a priori with sufficient precision, it can be applied iteratively with the prescribed cost value adjusted in an outer loop. Problem Formulation Let us consider the following simple optimal control problem stated in Mayer form: ^ 4 ( x ( t f ) , t f ) (1 )