Quasi-steady flight to quasi-steady flight transition for abort landing in a windshear: Trajectory optimization and guidance

This paper is concerned with the optimal transition and the near-optimum guidance of an aircraft from quasi-steady flight to quasi-steady flight in a windshear. The abort landing problem is considered with reference to flight in a vertical plane. In addition to the horizontal shear, the presence of a downdraft is considered.It is assumed that a transition from descending flight to ascending flight is desired; that the initial state corresponds to quasi-steady flight with absolute path inclination of –3.0 deg; and that the final path inclination corresponds to quasi-steady steepest climb. Also, it is assumed that, as soon as the shear is detected, the power setting is increased at a constant time rate until maximum power setting is reached; afterward, the power setting is held constant. Hence, the only control is the angle of attack. Inequality constraints are imposed on both the angle of attack and its time derivative.First, trajectory optimization is considered. The optimal transition problem is formulated as a Chebyshev problem of optimal control: the performance index being minimized is the peak value of the modulus of the difference between the instantaneous altitude and a reference value, assumed constant. By suitable transformations, the Chebyshev problem is converted into a Bolza problem. Then, the Bolza problem is solved employing the dual sequential gradient-restoration algorithm (DSGRA) for optimal control problems.Two types of optimal trajectories are studied, depending on the conditions desired at the final point. Type 1 is concerned with gamma recovery (recovery of the value of the relative path inclination corresponding to quasi-steady steepest climb). Type 2 is concerned with quasi-steady flight recovery (recovery of the values of the relative path inclination, the relative velocity, and the relative angle of attack corresponding to quasi-steady steepest climb). Both the Type 1 trajectory and the Type 2 trajectory include three branches: descending flight, nearly horizontal flight, and ascending flight. Also, for both the Type 1 trajectory and the Type 2 trajectory, descending flight takes place in the shear portion of the trajectory; horizontal flight takes place partly in the shear portion and partly in the aftershear portion of the trajectory; and ascending flight takes place in the aftershear portion of the trajectory. While the Type 1 trajectory and the Type 2 trajectory are nearly the same in the shear portion, they diverge to a considerable degree in the aftershear portion of the trajectory.Next, trajectory guidance is considered. Two guidance schemes are developed so as to achieve near-optimum transition from quasi-steady descending flight to quasi-steady ascending flight: acceleration guidance (based on the relative acceleration) and gamma guidance (based on the absolute path inclination).The guidance schemes for quasi-steady flight recovery in abort landing include two parts in sequence: shear guidance and aftershear guidance. The shear guidance is based on the result that the shear portion of the trajectory depends only mildly on the boundary conditions. Therefore, any of the guidance schemes already developed for Type 1 trajectories can be employed for Type 2 trajectories (descent guidance followed by recovery guidance). The aftershear guidance is based on the result that the aftershear portion of the trajectory depends strongly on the boundary conditions; therefore, the guidance schemes developed for Type 1 trajectories cannot be employed for Type 2 trajectories. For Type 2 trajectories, the aftershear guidance includes level flight guidance followed by ascent guidance. The level flight guidance is designed to achieve almost complete velocity recovery; the ascent guidance is designed to achieve the desired final quasi-steady state.The numerical results show that the guidance schemes for quasi-steady flight recovery yield a transition from quasi-steady flight to quasi-steady flight which is close to that of the optimal trajectory, allows the aircraft to achieve the final quasi-steady state, and has good stability properties.This research was supported by NASA Langley Research Center, Grant No. NAG-1-516, by Boeing Commercial Airplane Company, and by Air Line Pilots Association.The authors are indebted to Dr. R. L. Bowles (NASA-LRC) and Dr. G. R. Hennig (BCAC) for helpful discussions.     

Optimal penetration landing trajectories in the presence of windshear

This paper is concerned with optimal flight trajectories in the presence of windshear. The penetration landing problem is considered with reference to flight in a vertical plane, governed by either one control (the angle of attack, if the power setting is predetermined) or two controls (the angle of attack and the power setting). Inequality constraints are imposed on the angle of attack, the power setting, and their time derivatives.The performance index being minimized measures the deviation of the flight trajectory from a nominal trajectory. In turn, the nominal trajectory includes two parts: the approach part, in which the slope is constant; and the flare part, in which the slope is a linear function of the horizontal distance. In the optimization process, the time is free; the absolute path inclination at touchdown is specified; the touchdown velocity is subject to upper and lower bounds; and the touchdown distance is subject to upper and lower bounds.Three power setting schemes are investigated: (S1) maximum power setting; (S2) constant power setting; and (S3) control power setting. In Scheme (S1), it is assumed that, immediately after the windshear onset, the power setting is increased at a constant time rate until maximum power setting is reached; afterward, the power setting is held constant; in this scheme, the only control is the angle of attack. In Scheme (S2), it is assumed that the power setting is held at a constant value, equal to the prewindshear value; in this scheme, the only control is the angle of attack. In Scheme (S3), the power setting is regarded as a control, just as the angle of attack.Under the above conditions, the optimal control problem is solved by means of the primal sequential gradient-restoration algorithm (PSGRA). Numerical results are obtained for several combinations of windshear intensities and initial altitudes. The main conclusions are given below with reference to strong-to-severe windshears.In Scheme (S1), the touchdown requirements can be satisfied for relatively low initial altitudes, while they cannot be satisfied for relatively high initial altitudes; the major inconvenient is excess of velocity at touchdown. In Scheme (S2), the touchdown requirements cannot be satisfied, regardless of the initial altitude; the major inconvenient is defect of horizontal distance at touchdown.In Scheme (S3), the touchdown requirements can be satisfied, and the optimal trajectories exhibit the following characteristics: (i) the angle of attack has an initial decrease, which is followed by a gradual, sustained increase; the largest value of the angle of attack is attained near the end of the shear; in the aftershear region, the angle of attack decreases gradually; (ii) initially, the power setting increases rapidly until maximum power setting is reached; then, maximum power setting is maintained in the shear region; in the aftershear region, the power setting decreases gradually; (iii) the relative velocity decreases in the shear region and increases in the aftershear region; the point of minimum velocity occurs at the end of the shear; and (iv) depending on the windshear intensity and the initial altitude, the deviations of the flight trajectory from the nominal trajectory can be considerable in the shear region; however, these deviations become small in the aftershear region, and the optimal flight trajectory recovers the nominal trajectory.A comparison is shown between the optimal trajectories of Scheme (S3) and the trajectories arising from alternative guidance schemes, such as fixed controls (fixed angle of attack, coupled with fixed power setting) and autoland (angle of attack controlled via path inclination signals, coupled with power setting controlled via velocity signals). The superiority of the optimal trajectories of Scheme (S3) is shown in terms of the ability to meet the path inclination, velocity, and distance requirements at touchdown. Therefore, it is felt that guidance schemes based on the properties of the optimal trajectories of Scheme (S3) should prove to be superior to alternative guidance schemes, such as the fixed control guidance scheme and the autoland guidance scheme.Portions of this paper were presented at the AIAA 26th Aerospace Sciences Meeting, Reno, Nevada, January 11–14, 1988 (Paper No. AIAA-88-0580).This research was supported by NASA-Langley Research Center, Grant No. NAG-1-516, by Boeing Commercial Airplane Company (BCAC), and by Air Line Pilots Association (ALPA).The authors are indebted to Dr. R. L. Bowles, NASA-Langley Research Center, and to Dr. G. R. Hennig, Boeing Commercial Airplane Company, for helpful discussions.     

Optimal abort landing trajectories in the presence of windshear

This paper is concerned with optimal flight trajectories in the presence of windshear. The abort landing problem is considered with reference to flight in a vertical plane. It is assumed that, upon sensing that the airplane is in a windshear, the pilot increases the power setting at a constant time rate until maximum power setting is reached; afterward, the power setting is held constant. Hence, the only control is the angle of attack. Inequality constraints are imposed on both the angle of attack and its time derivative.The performance index being minimized is the peak value of the altitude drop. The resulting optimization problem is a minimax problem or Chebyshev problem of optimal control, which can be converted into a Bolza problem through suitable transformations. The Bolza problem is then solved employing the dual sequential gradient-restoration algorithm (DSGRA) for optimal control problems. Numerical results are obtained for several combinations of windshear intensities, initial altitudes, and power setting rates.For strong-to-severe windshears, the following conclusions are reached: (i) the optimal trajectory includes three branches: a descending flight branch, followed by a nearly horizontal flight branch, followed by an ascending flight branch after the aircraft has passed through the shear region; (ii) along an optimal trajectory, the point of minimum velocity is reached at about the time when the shear ends; (iii) the peak altitude drop depends on the windshear intensity, the initial altitude, and the power setting rate; it increases as the windshear intensity increases and the initial altitude increases; and it decreases as the power setting rate increases; (iv) the peak altitude drop of the optimal abort landing trajectory is less than the peak altitude drop of comparison trajectories, for example, the constant pitch guidance trajectory and the maximum angle of attack guidance trajectory; (v) the survival capability of the optimal abort landing trajectory in a severe windshear is superior to that of comparison trajectories, for example, the constant pitch guidance trajectory and the maximum angle of attack guidance trajectory.Portions of this paper were presented at the IFAC 10th World Congress, Munich, Germany, July 27–31, 1987 (Paper No. IFAC-87-9221).This research was supported by NASA Langley Research Center, Grant No. NAG-1-516, by Boeing Commercial Airplane Company (BCAC), and by Air Line Pilots Association (ALPA). Discussions with Dr. R. L. Bowles (NASA-LRC) and Mr. C. R. Higgins (BCAC) are acknowledged.     

Guidance strategies for near-optimum take-off performance in a windshear

This paper is concerned with guidance strategies for near-optimum performance in a windshear. This is a wind characterized by sharp change in intensity and direction over a relatively small region of space. The take-off problem is considered with reference to flight in a vertical plane.First, trajectories for optimum performance in a windshear are determined for different windshear models and different windshear intensities. Use is made of the methods of optimal control theory in conjunction with the dual sequential gradient-restoration algorithm (DSGRA) for optimal control problems. In this approach, global information on the wind flow field is needed.Then, guidance strategies for near-optimum performance in a wind-shear are developed, starting from the optimal trajectories. Specifically, three guidance schemes are presented: (A) gamma guidance, based on the relative path inclination; (B) theta guidance, based on the pitch attitude angle; and (C) acceleration guidance, based on the relative acceleration. In this approach, local information on the wind flow field is needed.Next, several alternative schemes are investigated for the sake of completeness, more specifically: (D) constant alpha guidance; (E) constant velocity guidance; (F) constant theta guidance; (G) constant relative path inclination guidance; (H) constant absolute path inclination guidance; and (I) linear altitude distribution guidance.Numerical experiments show that guidance schemes (A)–(C) produce trajectories which are quite close to the optimum trajectories. In addition, the near-optimum trajectories associated with guidance schemes (A)–(C) are considerably superior to the trajectories arising from the alternative guidance schemes (D)–(I).An important characteristic of guidance schemes (A)–(C) is their simplicity. Indeed, these guidance schemes are implementable using available instrumentation and/or modification of available instrumentation.Portions of this were presented at the AIAA 24th Aerospace Sciences Meeting, Reno, Nevada, January 6–9, 1986. The authors are indebted to Boeing Commercial Aircraft Company, Seattle, Washington and to Pratt and Whittney Aircraft, East Hartford, Connecticut for supplying some of the technical data pertaining to this study.The authors are indebted to Dr. R. L. Bowles, NASA-Langley Research Center, Hampton, Virginia for helpful discussions. They are also indebted to Mr. Z. G. Zhao, Aero-Astronautics Group, Rice University, Houston, Texas for analytical and computational assistance.This research was supported by NASA-Langley Research Center, Grant No. NAG-1-516. This paper, a continuation of Ref.1, is based in part on Refs. 2–3.     

Optimal take-off trajectories in the presence of windshear

This paper is concerned with optimal flight trajectories in the presence of windshear. With particular reference to take-off, eight fundamental optimization problems [Problems (P1)–(P8)] are formulated under the assumptions that the power setting is held at the maximum value and that the airplane is controlled through the angle of attack.Problems (P1)–(P3) are least-square problems of the Bolza type. Problems (P4)–(P8) are minimax problems of the Chebyshev type, which can be converted into Bolza problems through suitable transformations. These problems are solved employing the dual sequential gradient-restoration algorithm (DSGRA) for optimal control problems.Numerical results are obtained for a large number of combinations of performance indexes, boundary conditions, windshear models, and windshear intensities. However, for the sake of brevity, the presentation of this paper is restricted to Problem (P6), minimax h, and Problem (P7), minimax . Inequality constraints are imposed on the angle of attack and the time derivative of the angle of attack.The following conclusions are reached: (i) optimal trajectories are considerably superior to constant-angle-of-attack trajectories; (ii) optimal trajectories achieve minimum velocity at about the time when the windshear ends; (iii) optimal trajectories can be found which transfer an aircraft from a quasi-steady condition to a quasi-steady condition through a windshear; (iv) as the boundary conditions are relaxed, a higher final altitude can be achieved, albeit at the expense of a considerable velocity loss; (v) among the optimal trajectories investigated, those solving Problem (P7) are to be preferred, because the altitude distribution exhibits a monotonic behavior; in addition, for boundary conditions BC2 and BC3, the peak angle of attack is below the maximum permissible value; (vi) moderate windshears and relatively severe windshears are survivable employing an optimized flight strategy; however, extremely severe windshears are not survivable, even employing an optimized flight strategy; and (vii) the sequential gradient-restoration algorithm (SGRA), employed in its dual form (DSGRA), has proven to be a powerful algorithm for solving the problem of the optimal flight trajectories in a windshear.Portions of this paper were presented at the AIAA Atmospheric Flight Mechanics Conference, Snowmass, Colorado, August 19–21, 1985. The authors are indebted to Boeing Commercial Aircraft Company, Seattle, Washington and to Pratt and Whitney Aircraft, East Hartford, Connecticut for supplying some of the technical data pertaining to this study.This research was supported by NASA-Langley Research Center, Grant No. NAG-1-516. The authors are indebted to Dr. R. L. Bowles, NASA-Langley Research Center, Hampton, Virginia, for helpful discussions.This paper is based in part on Refs. 1–5.     

Maximum survival capability of an aircraft in a severe windshear

This paper is concerned with guidance strategies and piloting techniques which ensure near-optimum performance and maximum survival capability in a severe windshear. The take-off problem is considered with reference to flight in a vertical plane. In addition to the horizontal shear, the presence of a downdraft is assumed.First, six particular guidance schemes are considered, namely: constant alpha guidance; maximum alpha guidance; constant velocity guidance; constant absolute path inclination guidance; constant rate of climb guidance; and constant pitch guidance. Among these, it is concluded that the best one is the constant pitch guidance.Next, in an effort to improve over the constant pitch guidance, three additional trajectories are considered: the optimal trajectory, which minimizes the maximum deviation of the absolute path inclination from a reference value, while employing global information on the wind flow field; the gamma guidance trajectory, which is based on the absolute path inclination and which approximates the behavior of the optimal trajectory, while employing local information on the windshear and the downdraft; and the simplified gamma guidance trajectory, which is the limiting case of the gamma guidance trajectory in a severe windshear and which does not require precise information on the windshear and the downdraft.The essence of the simplified gamma guidance trajectory is that it yields a quick transition to horizontal flight. Comparative numerical experiments show that the survival capability of the simplified gamma guidance trajectory is superior to that of the constant pitch trajectory and is close to that of the optimal trajectory.Next, with reference to the simplified gamma guidance trajectory, the effect of the feedback gain coefficient is studied. It is shown that larger values of the gain coefficient improve the survival capability in a severe windshear; however, excessive values of the gain coefficient are undesirable, because they result in larger altitude oscillations and lower average altitude.Finally, with reference to the simplified gamma guidance trajectory, the effect of time delays is studied, more specifically, the time delay ₁ in reacting to windshear onset and the time delay ₂ in reacting to windshear termination. While time delay ₂ has little effect on survival capability, time delay ₁ appears to be critical in the following sense: smaller values of ₁ correspond to better survival capability in a severe windshear, while larger values of ₁ are associated with a worsening of the survival capability in a severe windshear.This research was supported by NASA-Langley Research Center, Grant No. NAG-1-516, and by Boeing Commercial Airplane Company.     

Decomposition technique and optimal trajectories for the aeroassisted flight experiment

The aeroassisted flight experiment (AFE) refers to a spacecraft to be launched and then recovered by the space shuttle in 1994. It simulates a transfer from a geosynchronous Earth orbit (GEO) to a low Earth orbit (LEO). Specifically, the AFE spacecraft is released from the space shuttle and is accelerated by means of a solid rocket motor toward Earth, so as to achieve atmospheric entry conditions close to those of a spacecraft returning from GEO. Following the atmospheric pass, the AFE spacecraft ascends to the specified LEO via an intermediate parking Earth orbit (PEO). The final maneuver includes the rendezvous with and the capture by the space shuttle. The entry and exit orbital planes of the AFE spacecraft are identical with the orbital plane of the space shuttle. In this paper, with reference to the AFE spacecraft, an actual GEO-to-LEO transfer is considered and optimal trajectories are determined by minimizing the total characteristic velocity. The optimization is performed with respect to the time history of the controls (angle of attack and angle of bank), the entry path inclination and the flight time being free. Two transfer maneuvers are considered: (DA) direct ascent to LEO; (IA) indirect ascent to LEO via PEO. While the motion of the AFE spacecraft in a 3D-space is described by a system of six ODEs, substantial simplifications are possible if one exploits these facts: (i) the instantaneous orbital plane is nearly identical with the initial orbital plane; (ii) the bank angle is small; and (iii) the Earth's angular velocity is relatively small. Under these assumptions, the complete system can be decoupled into two subsystems, one describing the longitudinal motion and one describing the lateral motion. The angle of attack history, the entry path inclination, and the flight time are determined via the longitudinal motion subsystem; in this subsystem, the total characteristic velocity is minimized subject to the specified LEO requirement. The angle of bank history is determined via the lateral motion subsystem; in this subsystem, the difference between the instantaneous bank angle and a constant bank angle is minimized in the least square sense subject to the specified orbital inclination requirement. It is shown that both the angle of attack and the angle of bank are constant. This result has considerable importance in the design of nominal trajectories to be used in the guidance of AFE and AOT vehicles.     

Optimal Starting Conditions for the Rendezvous Maneuver, Part 1: Optimal Control Approach

We consider the three-dimensional rendezvous between two spacecraft: a target spacecraft on a circular orbit around the Earth and a chaser spacecraft initially on some elliptical orbit yet to be determined. The chaser spacecraft has variable mass, limited thrust, and its trajectory is governed by three controls, one determining the thrust magnitude and two determining the thrust direction. We seek the time history of the controls in such a way that the propellant mass required to execute the rendezvous maneuver is minimized. Two cases are considered: (i) time-to-rendezvous free and (ii) time-to-rendezvous given, respectively equivalent to (i) free angular travel and (ii) fixed angular travel for the target spacecraft. The above problem has been studied by several authors under the assumption that the initial separation coordinates and the initial separation velocities are given, hence known initial conditions for the chaser spacecraft. In this paper, it is assumed that both the initial separation coordinates and initial separation velocities are free except for the requirement that the initial chaser-to-target distance is given so as to prevent the occurrence of trivial solutions. Analyses performed with the multiple-subarc sequential gradient-restoration algorithm for optimal control problems show that the fuel-optimal trajectory is zero-bang, namely it is characterized by two subarcs: a long coasting zero-thrust subarc followed by a short powered max-thrust braking subarc. While the thrust direction of the powered subarc is continuously variable for the optimal trajectory, its replacement with a constant (yet optimized) thrust direction produces a very efficient guidance trajectory: Indeed, for all values of the initial distance, the fuel required by the guidance trajectory is within less than one percent of the fuel required by the optimal trajectory. For the guidance trajectory, because of the replacement of the variable thrust direction of the powered subarc with a constant thrust direction, the optimal control problem degenerates into a mathematical programming problem with a relatively small number of degrees of freedom, more precisely: three for case (i) time-to-rendezvous free and two for case (ii) time-to-rendezvous given. In particular, we consider the rendezvous between the Space Shuttle (chaser) and the International Space Station (target). Once a given initial distance SS-to-ISS is preselected, the present work supplies not only the best initial conditions for the rendezvous trajectory, but simultaneously the corresponding final conditions for the ascent trajectory.     

Aircraft Take-Off in Windshear: A Viability Approach

This paper is devoted to the analysis of aircraft dynamics during tage-off in the presence of windshear. We formulate the take-off problem as a differential game against Nature. Here, the first player is the relative angle of attack of the aircraft (considered as the control variable) and the second player is the disturbance caused by a windshear. We impose state constraints on the state variables of the game, which represents aircraft safety constraints (minimum altitude, given altitude rate). By using viability theory, we address the question of existence of an open loop control assuring a viable trajectory (i.e. satisfying the state constraints) no matter the disturbance is, i.e. for all admissible disturbances causeed by the windshear. Through numerical simulations of the viability kernel algorithm, we demonstrate the capabilities of this approach for determining safe flight domains of an aircraft during take-off within windshear.     

Optimal trajectories for aeroassisted,coplanar orbital transfer

This paper considers both classical and minimax problems of optimal control which arise in the study of aeroassisted, coplanar orbital transfer. The maneuver considered involves the coplanar transfer from a high planetary orbit to a low planetary orbit. An example is the HEO-to-LEO transfer of a spacecraft, where HEO denotes high Earth orbit and LEO denotes low Earth orbit. In particular, HEO can be GEO, a geosynchronous Earth orbit.The basic idea is to employ the hybrid combination of propulsive maneuvers in space and aerodynamic maneuvers in the sensible atmosphere. Hence, this type of flight is also called synergetic space flight. With reference to the atmospheric part of the maneuver, trajectory control is achieved by means of lift modulation. The presence of upper and lower bounds on the lift coefficient is considered.Within the framework of classical optimal control, the following problems are studied: (P1) minimize the energy required for orbital transfer; (P2) minimize the time integral of the heating rate; (P3) minimize the time of flight during the atmospheric portion of the trajectory; (P4) maximize the time of flight during the atmospheric portion of the trajectory; (P5) minimize the time integral of the square of the path inclination; and (P6) minimize the sum of the squares of the entry and exit path inclinations. Problems (P1) through (P6) are Bolza problems of optimal control.Within the framework of minimax optimal control, the following problems are studied: (Q1) minimize the peak heating rate; (Q2) minimize the peak dynamic pressure; and (Q3) minimize the peak altitude drop. Problems (Q1) through (Q3) are Chebyshev problems of optimal control, which can be converted into Bolza problems by suitable transformations.Numerical solutions for Problems (P1)–(P6) and Problems (Q1)–(Q3) are obtained by means of the sequential gradient-restoration algorithm for optimal control problems. The engineering implications of these solutions are discussed. In particular, the merits of nearly-grazing trajectories are considered.This research was supported by the Jet Propulsion Laboratory, Contract No. 956415. The authors are indebted to Dr. K. D. Mease, Jet Propulsion Laboratory, for helpful discussions. This paper is a condensation of the investigation reported in Ref. 1.     

Optimal Trajectories and Guidance Schemes for Ship Collision Avoidance

The best strategy for collision avoidance under emergency conditions is to maximize wrt the controls the timewise minimum distance between the host ship and the intruder ship. In a restricted waterway area, two main constraints must be satisfied: the lateral deviation of the host ship from the original course is to be contained within certain limits; the longitudinal distance covered by the host ship is to be subject to a prescribed bound. At the maximin point of the encounter, the time derivative of the relative distance vanishes; this yields an inner boundary condition (orthogonality between the relative position vector and the relative velocity vector) separating the main phases of the maneuver: the avoidance and recovery phases. In this way, the optimal trajectory problem (a Chebyshev problem) can be converted into a Bolza problem with an inner boundary condition. Numerical solutions are obtained via the multiple-subarc sequential gradient-restoration algorithm (SGRA). Because the optimal trajectory is not suitable for real-time implementation, a guidance scheme approximating the optimal trajectory in real time is to be developed. For ship collision avoidance, the optimal trajectory results show that the rudder angle time history has a bang-bang form characterized by the alternation of saturated control subarcs of opposite signs joined by rapid transitions. Just as the optimal trajectory can be partitioned into three phases (avoidance phase, recovery phase, steady phase), a guidance trajectory can be constructed in the same way. For the avoidance and recovery phases, use of decomposition techniques leads to an algorithm computing the time lengths of these phases in real time. For the steady phase, a feedback control scheme is used to maneuver the ship steadily. Numerical results are presented. Portions of this paper were presented by the senior author at the 13th International Workshop on Dynamics and Control, Wiesensteig, Germany, 22-26 May 2005, in honor of George Leitmann. This research was supported by NSF Grant CMS-02-18878.     

Guidance Trajectories for Spacecraft Rendezvous

In a previous paper of Miele et al. (J. Optim. Theory Appl. 132(1), 2007), we employed the single-subarc sequential gradient-restoration algorithm to optimize the three-dimensional rendezvous between a target spacecraft in a planar circular orbit and a chaser spacecraft with an initial separation distance and separation velocity. The achieved continuous solutions are characterized by two, three, or four subarcs depending on the performance index (time, fuel) and the constraints. In this paper, based on the solutions in Miele et al. (J. Optim. Theory Appl. 132(1), 2007), we employ the multiple-subarc sequential gradient-restoration algorithm to produce pieced guidance trajectories implementable in real time via constant control components. In other words, in this investigation, we force the controls to behave as parameters in each subarc. With the above understanding, we investigate four problems: (P1) minimum time, fuel free; (P2) minimum fuel, time free; (P3) minimum time, fuel given; (P4) minimum fuel, time given. Problem P1 results in a two-subarc solution, each subarc with constant controls: a max-thrust accelerating subarc followed by a max-thrust braking subarc. Problem P2 results in a four-subarc solution, each subarc with constant controls: an initial coasting subarc, followed by a max-thrust braking subarc, followed by another coasting subarc, followed by another max-thrust braking subarc. Problems P3 and P4 result in two, three, or four-subarc solutions depending on the performance index and the constraints, albeit with constant controls in each subarc. For each of the problems studied, the performance index of the multiple-subarc pieced guidance trajectory approximates well the performance index of the single-subarc continuous optimal trajectory of Miele et al. (J. Optim. Theory Appl. 132(1), 2007) as well as the performance index of the multiple-subarc pieced optimal trajectory: the pairwise relative differences in performance index are less than 1/100. This research was supported by NSF under Grant CMS-0218878.     

Abort landing in windshear: Optimal control problem with third-order state constraint and varied switching structure

Optimal abort landing trajectories of an aircraft under different windshear-downburst situations are computed and discussed. In order to avoid an airplane crash due to severe winds encountered by the aircraft during the landing approach, the minimum altitude obtained during the abort landing maneuver is to be maximized. This maneuver is mathematically described by a Chebyshev optimal control problem. By a transformation to an optimal control problem of Mayer type, an additional state variable inequality constraint for the altitude has to be taken into account; here, its order is three. Due to this altitude constraint, the optimal trajectories exhibit, depending on the windshear parameters, up to four touch points and also up to one boundary arc at the minimum altitude level. The control variable is the angle of attack time rate which enters the equations of motion linearly; therefore, the Hamiltonian of the problem is nonregular. The switching structures also includes up to three singular subarcs and up to two boundary subarcs of an angle of attack constraint of first order. This structure can be obtained by applying some advanced necessary conditions of optimal control theory in combination with the multiple-shooting method. The optimal solutions exhibit an oscillatory behavior, reaching the minimum altitude level several times. By the optimization, the maximum survival capability can also be determined; this is the maximum wind velocity difference for which recovery from windshear is just possible. The computed optimal trajectories may serve as benchmark trajectories, both for guidance laws that are desirable to approach in actual flight and for optimal trajectories may then serve as benchmark trajectories both for guidance schemes and also for numerical methods for problems of optimal control.This paper is dedicated to Professor George Leitmann on the occasion of his seventieth birthday.     

Optimal Control of a Ship for Course Change and Sidestep Maneuvers

We consider a ship subject to kinematic, dynamic, and moment equations and steered via rudder under the assumptions that the rudder angle and rudder angle time rate are subject to upper and lower bounds. We formulate and solve four Mayer problems of optimal control, the optimization criterion being the minimum time.Problems P1 and P2 deal with course change maneuvers. In Problem P1, a ship initially in quasi-steady state must reach the final point with a given yaw angle and zero yaw angle time rate. Problem P2 differs from Problem P1 in that the additional requirement of quasi-steady state is imposed at the final point.Problems P3 and P4 deal with sidestep maneuvers. In Problem P3, a ship initially in quasi-steady state must reach the final point with a given lateral distance, zero yaw angle, and zero yaw angle time rate. Problem P4 differs from Problem P3 in that the additional requirement of quasi-steady state is imposed at the final point.The above Mayer problems are solved via the sequential gradient-restoration algorithm in conjunction with a new singularity avoiding transformation which accounts automatically for the bounds on rudder angle and rudder angle time rate.The optimal control histories involve multiple subarcs along which either the rudder angle is kept at one of the extreme positions or the rudder angle time rate is held at one of the extreme values. In problems where quasi-steady state is imposed at the final point, there is a higher number of subarcs than in problems where quasi-steady state is not imposed; the higher number of subarcs is due to the additional requirement that the lateral velocity and rudder angle vanish at the final point.     

Optimal Starting Conditions for the Rendezvous Maneuver,Part 2: Mathematical Programming Approach

In a companion paper (Part 1, J. Optim. Theory Appl. 137(3), [2008]), we determined the optimal starting conditions for the rendezvous maneuver using an optimal control approach. In this paper, we study the same problem with a mathematical programming approach. Specifically, we consider the relative motion between a target spacecraft in a circular orbit and a chaser spacecraft moving in its proximity as described by the Clohessy-Wiltshire equations. We consider the class of multiple-subarc trajectories characterized by constant thrust controls in each subarc. Under these conditions, the Clohessy-Wiltshire equations can be integrated in closed form and in turn this leads to optimization processes of the mathematical programming type. Within the above framework, we study the rendezvous problem under the assumption that the initial separation coordinates and initial separation velocities are free except for the requirement that the initial chaser-to-target distance is given. In particular, we consider the rendezvous between the Space Shuttle (chaser) and the International Space Station (target). Once a given initial distance SS-to-ISS is preselected, the present work supplies not only the best initial conditions for the rendezvous trajectory, but simultaneously the corresponding final conditions for the ascent trajectory.     

Wind identification along a flight trajectory,part 3: 2D-dynamic approach

This paper deals with the identification of the wind profile along a flight trajectory by means of a two-dimensional dynamic approach. In this approach, the wind velocity components are computed as the difference between the inertial velocity components and the airspeed components. The airspeed profile as well as the nominal thrust, drag, and lift profiles are obtained from the available DFDR measurements. The actual values of the thrust, drag, and lift are assumed to be proportional to the respective nominal values via multiplicative parameters, called the thrust, drag, and lift factors. The thrust, drag, and lift factors plus the inertial velocity components at impact are determined by matching the flight trajectory computed from DFDR data with the flight trajectory available from ATCR data. This leads to a least-square problem which is solved analytically under the additional requirement of closeness of the multiplicative factors to unity. Application of the 2D-dynamic approach to the case of Flight Delta 191 shows that, with reference to the last 180 sec before impact, the values of the multiplicative factors were 1.09, 0.84, and 0.89; this implies that the actual values of the thrust, drag, and lift were 9% above, 16% below, and 11% below their respective nominal values. For the last 60 sec before impact, the aircraft was subject to severe windshear, characterized by a horizontal wind velocity difference of 123 fps and a vertical wind velocity difference of 80 fps. The 2D-dynamic approach is applicable to the analysis of windshear accidents in take-off or landing, especially for the case of older-generation, shorter-range aircraft which do not carry the extensive instrumentation of newer-generation, longer-range aircraft. The same methodology can be extended to the investigation of aircraft accidents originating from causes other than windshear (e.g., icing, incorrect flap position, engine malfunction), above all if its precision is further increased by combining the 2D-dynamic approach and the 2D-kinematic approach.     

Aircraft control for flight in an uncertain environment: Takeoff in windshear

The design of the control of an aircraft encountering windshear after takeoff is treated as a problem of stabilizing the climb rate about a desired value of the climb rate. The resulting controller is a feedback one utilizing only climb rate information. Its robustness vis-a-vis windshear structure and intensity is illustrated via simulations employing four different windshear models.Notations ARL aircraft reference line - D drag force, lb - g gravitational force per unit mass=const, ft sec^–2 - h vertical coordinate of aircraft center of mass (altitude), ft - L lift force, lb - m aircraft mass=const, lb ft^–1 sec² - O mass center of aircraft - S reference surface, ft² - t time, sec - T thrust force, lb - V aircraft speed relative to wind-based reference frame, ft sec^–1 - V _e aircraft speed relative to ground, ft sec^–1 - W _x horizontal component of wind velocity, ft sec^–1 - W _h vertical component of wind velocity, ft sec^–1 - x horizontal coordinate of aircraft center of mass, ft - relative angle of attack, rad - relative path inclination, rad - _e path inclination, rad - thrust inclination, rad - air density=const, lb ft² sec² Dot denotes time derivative.     

Optimal Control of a Ship for a Course-Changing Maneuver

The steering control of a ship during a course-changing maneuver is formulated as a Bolza optimal control problem, which is solved via the sequential gradient-restoration algorithm (SGRA). Nonlinear differential equations describing the yaw dynamics of a steering ship are employed as the differential constraints, and both amplitude and slew rate limits on the rudder are imposed. Two performance indices are minimized: one measures the time integral of the squared course deviation between the actual ship course and a target course; the other measures the time integral of the absolute course deviation. Numerical results indicate that a smooth transition from the initial set course to the target course is achievable, with a trade-off between the speed of response and the amount of course angle overshoot.     

Multiple-Subarc Gradient-Restoration Algorithm, Part 2: Application to a Multistage Launch Vehicle Design

In Part 1 (see Ref. 2), a multiple-subarc gradient-restoration algorithm (MSGRA) was developed with the intent of enhancing the robustness of gradient-restoration algorithms and also enlarging the field of applications. Indeed, MSGRA can be applied to optimal control problems involving multiple subsystems as well as discontinuities in the state and control variables at the interface between contiguous subsystems.In Part 2 (this paper), MSGRA is applied to compute the optimal trajectory for a multistage launch vehicle design, specifically, a rocket-powered spacecraft ascending from the Earth surface to a low Earth orbit (LEO). Single-stage, double-stage, and triple-stage configurations are considered. For multistage configurations, discontinuities in the mass occur at the interfaces between consecutive stages.The numerical results show that, given the current levels of the engine specific impulse and spacecraft structural factor, the single-stage version is not feasible at this time, while the double-stage and triple-stage versions are feasible. Further increases in the specific impulse and decreases in the structural factor are needed if the single-stage configuration has to become feasible.Also, the numerical results show that the optimal trajectory requires initially maximum thrust, followed by modulated thrust so as to satisfy the maximum acceleration constraint, followed by nearly zero thrust for coasting flight, followed by a final burst with moderate thrust so as to increase the spacecraft velocity to the circular velocity needed for LEO insertion. The above properties of the optimal thrust time history are useful for developing the guidance scheme approximating in real time the optimal trajectory for a launch vehicle design.     

Abort landing in the presence of windshear as a minimax optimal control problem,part 2: Multiple shooting and homotopy

In Part 1 of the paper (Ref. 2), we have shown that the necessary conditions for the optimal control problem of the abort landing of a passenger aircraft in the presence of windshear result in a multipoint boundary-value problem. This boundary-value problem is especially well suited for numerical treatment by the multiple shooting method. Since this method is basically a Newton iteration, initial guesses of all variables are needed and assumptions about the switching structure have to be made. These are big obstacles, but both can be overcome by a so-called homotopy strategy where the problem is imbedded into a one-parameter family of subproblems in such a way that (at least) the first problem is simple to solve. The solution data to the first problem may serve as an initial guess for the next problem, thus resulting in a whole chain of problems. This process is to be continued until the objective problem is reached.Techniques are presented here on how to handle the various changes of the switching structure during the homotopy run. The windshear problem, of great interest for safety in aviation, also serves as an excellent benchmark problem: Nearly all features that can arise in optimal control appear when solving this problem. For example, the candidate for an optimal trajectory of the minimax optimal control problem shows subarcs with both bang-bang and singular control functions, boundary arcs and touch points of two state constraints, one being of first order and the other being of third order, etc. Therefore, the results of this paper may also serve as some sort of user's guide for the solution of complicated real-life optimal control problems by multiple shooting.The candidate found for an optimal trajectory is discussed and compared with an approximate solution already known (Refs. 3–4). Besides the known necessary conditions, additional sharp necessary conditions based on sign conditions of certain multipliers are also checked. This is not possible when using direct methods.An extended abstract of this paper was presented at the 8th IFAC Workshop on Control Applications of Nonlinear Programming and Optimization, Paris, France, 1989 (see Ref. 1).This paper is dedicated to Professor Hans J. Stetter on the occasion of his 60th birthday.