Abstract

A novel standard trajectory design and tracking guidance used in the multiple active leap maneuver mode for hypersonic glide vehicles (HGVs) is proposed in this paper. First, the dynamic equation and multiconstraint model are first established in the flight path coordinate system. Second, the reference drag acceleration-normalized energy (D-e) profile of the multiple active leap maneuver mode is quickly determined by the Newton iterative algorithm with a single design parameter. The range to go error is corrected by the drag acceleration profile update algorithm, and the drag acceleration error of the gliding terminal is corrected by the aerodynamic parameter estimation algorithm. Then, the reference drag acceleration tracking guidance law is designed based on the prescribed performance control method. Finally, the CAV-L vehicle model is used for numerical simulation. The results show that the proposed method can satisfy the design requirements of drag acceleration under multiple active leap maneuver modes, and the reference drag acceleration can be tracked precisely. The adaptability and robustness of the proposed method are verified by the Monte Carlo simulations under various combined deviation conditions.

1. Introduction

An HGV is a kind of vehicle that achieves long-range nonballistic re-entry by performing unpowered gliding motion in near space. HGVs possess many characteristics, such as high flight speeds, strong breakthrough maneuver abilities, and wide coverage of airspaces, making them a popular topic in the current research and development of new equipment [1, 2]. However, with the rapid development of relevant technologies in air and missile defense systems and the gradual improvement of equipment systems, there have been certain breakthroughs in trajectory tracking and prediction methods for the gliding phase of HGVs [37], as well as defensive interception technologies [8, 9]. These advancements pose a certain threat and challenge to the survivability of HGVs. Therefore, researching the utilization of the high lift-to-drag ratio characteristics of HGVs to achieve multiple extended range maneuver within the atmosphere can reduce the tracking and prediction accuracy of defense systems as well as interception probabilities [10]. This is important for enhancing the survivability and breakthrough capability of HGVs.

Regarding the trajectory design for the gliding phase of HGVs, existing trajectory design methods include standard trajectory guidance, predictive correction guidance, and closed-loop optimal guidance [11]. The standard trajectory guidance law is more developed and widely used. It is typically implemented by selecting specific flight profiles, such as the drag acceleration-velocity (D-V) [12], drag acceleration-normalized energy (D-e) [13], and altitude-velocity (H-V) [14]. Regarding the drag acceleration-energy profile, Gifty et al. [15] designed a drag acceleration curve with an analytical range form by using a multi-segment-modulated cubic polynomial. Sagliano and Mooij [16] designed the drag acceleration curve by the convex optimization method, and Zhang et al. [17] designed the quadratic polynomial and linear drag acceleration curve. The aforementioned trajectory design methods have low computational complexities but do not account for the requirements for gliding phase leap maneuver.

Aiming at the problem of leap maneuver trajectory design for the gliding phase of HGVs, Li et al. [18] used a stochastic gradient particle swarm optimization algorithm to design the drag acceleration curve. However, their method has high computational complexity. Regarding the glide maneuver trajectory design problem, Zhang et al. [19] divided the glide phase maneuver into two strategies, longitudinal maneuver and lateral maneuver. They used sequential quadratic programming methods to design five different maneuver trajectory modes, including longitudinal serpentine maneuvers. In relation to lateral maneuver strategy, Zhu et al. [20] incorporated lateral maneuver paths into the trajectory model and optimized the trajectory. Liang et al. [21] introduced the concept of maneuver coefficient and designed a three-dimensional trajectory planning method for decoupling the longitudinal and lateral planes by using the maneuver coefficient to describe the strength of lateral maneuver. For the longitudinal maneuver strategy, Tan et al. [22] used the adaptive pseudospectral method to obtain the optimal dive trajectory considering the path constraints. An et al. [23] applied predictive correction guidance methods to longitudinal trajectory design and obtained the gliding phase ballistic trajectory form through offline calculations. Zhu et al. [24] obtained the optimal maneuver form for the dive section based on the optimal control method with minimum energy consumption and line-of-sight angular velocity as the performance index in the dive stage. Aiming at the planar tripartite pursuit problem model, Hu et al. [25] solved the optimal breakout strategy based on differential countermeasure theory, which realized the breakout strike mission requirements to a certain extent. For the case of antidefense of two interceptors, Shen et al. [26] solved the surprise ballistic design problem by the second-order cone planning method and obtained an antidefense strategy based on the initial line-of-sight angle of the interceptors. Jiang et al. [27] studied a hypersonic gliding vehicle antidefense method based on deep reinforcement learning, which abstracts the confrontation process into a generalized three-body confrontation optimization problem and generates a breakout guidance law and the corresponding breakout trajectory through data training. From the perspective of escaping the KKV maneuver interception range, Liu et al. [28] first analysed the feasibility of anti-head-on interception for hypersonic glide vehicle and then designed a maneuver strategy in resolved form based on the state of the KKV at the moment of its separation from the booster stage. However, the abovementioned maneuver trajectory design method has insufficient adaptability to online trajectory design due to the real-time performance of the optimization algorithm and the difficulty of adjusting the maneuver state.

To address the design issues of the gliding phase guidance method for vehicle and to achieve controllable active leap maneuver during the gliding flight process, it is necessary to track the reference standard states designed based on optimization conditions considering process and terminal constraints. This involves generating high-precision tracking guidance commands to eliminate deviations between the actual trajectory and the standard trajectory. Mease and Kremer applied the differential geometric feedback linearization theory [29] and the evolved acceleration guidance logic for entry (EAGLE) [30] method for tracking the drag acceleration standard profile of spaceplanes, achieving excellent control results. The traditional proportional-differential guidance law [31] has limited adaptability, and the guidance parameters need to be reset for different tasks. To improve the guidance parameter sensitivity, Liang et al. [32] used the incremental nonlinear dynamic inversion method to track the design reference state, only calculating the needed guidance command increments at a given time. However, since the dynamic inversion method requires an accurate mathematical model of the object, it is sensitive to deviations. To improve the robustness of the guidance law, the sliding mode control method has been introduced into the design of the gliding tracking guidance law. An et al. [23] designed a finite-time convergent sliding mode guidance law according to the tracking error, and Li et al. [33] designed a sliding mode guidance law combining the global integral sliding mode surface and the exponential form reaching law. The above method can achieve good state tracking ability through reasonable parameter design under standard conditions, but there are some shortcomings, such as repeated formulation of guidance parameters and the inability to constrain dynamic performance. To solve the above problems, Bechlioulis and Rovithakis [34] and Chandramohan and Calise [35] studied the prescribed performance control method. This method introduces the prescribed performance function and the error conversion function to constrain the transient performance and steady-state performance of the system at the same time. This method transforms the tracking error problem into a uniformly bounded problem of the prescribed performance function, which reduces the sensitivity of the actual system to the control parameters. In addition, linear quadratic regulator (LQR) [36], optimal control, adaptive control, and other methods are also widely used in the study of standard trajectory tracking laws.

In this paper, an analytical design method of drag accelerations and a robust gliding guidance law that can realize multiple active leap maneuver are studied. First, in the drag acceleration-normalized energy (D-e) profile, the maneuver acceleration in the form of a polyline is rapidly designed by single-parameter iteration. Under the premise of satisfying the terminal and path constraints, the longitudinal leap maneuver is realized by using the change in drag acceleration in the form of a polyline. Second, to ensure that the standard reference drag acceleration meets the requirements of the terminal state constraint under the condition of deviations, the drag acceleration profile update algorithm is used for rapid correction. Afterward, the reference drag acceleration tracking guidance law is designed based on the sliding mode and prescribed performance control method. The tracking error convergence problem is transformed into a uniformly bounded problem, which improves the robustness of the gliding guidance law. Finally, through numerical simulation and the Monte Carlo calculations, the proposed method is verified under nominal conditions and disturbance conditions.

2. Gliding Phase Formulation

2.1. Coordinate System

Since the descriptions of the state parameters such as the force, position, and velocity of an HGV are usually in different coordinate systems, the dynamic model is established in the flight path coordinate system to ensure that analyses and calculations are conducted under the same space-time reference.

The earth-centered inertial coordinate system is a geocentric coordinate system with its origin at the center of the earth. The axis is the intersection of equatorial and ecliptic planes. The axis is the rotation direction of the earth. The axis satisfies the right-hand rule. The origin of the coordinate system is the center of the vehicle mass. The axis is the direction of the flight velocity. The axis is located in the plane formed by the velocity vector and the vector from the center of the earth to the vehicle, perpendicular to the axis , and positive in the upward direction. The axis satisfies the right-hand rule, as shown in Figure 1.

2.2. Dynamic Equations

It is assumed that the vehicle has no power, and the mass is constant during the gliding phase. The sideslip angle is zero. The earth is assumed to be a uniform sphere rotating around its own axis, without considering the nonspherical perturbation caused by the earth’s oblateness. The atmosphere is assumed to be stationary relative to the earth and uniform at the same altitude. The 3DOF dynamic equation of the vehicle is established in the flight path coordinate system as follows:

The following six state variables are included in the above dynamic equations: . denotes the distance from the center of the earth to the center of mass of the vehicle. is the longitude. is the latitude. is the earth-relative velocity. and are the flight path angle and velocity heading angle measured from the north in a clockwise direction, respectively. is the bank angle, which is one of the two control variables. is the gravitational acceleration with a value of 9.8 m/s2. are the additional terms caused by the rotation of the earth, and their sizes are represented as follows: where is the rotational angular velocity of the earth with a value of .

The terms and are the aerodynamic drag and lift, respectively, and they are expressed as follows: where and are drag and lift coefficients that are functions of and the Mach number . is the atmospheric density and can be calculated by , where is the atmospheric density at sea level and is the scale altitude. is the altitude, and is the radius of the earth. is the reference area, and is the vehicle mass. Since and are directly related to the angle of attack of the vehicle and the bank angle can control the component of the lift in the longitudinal plane, it can be seen from Eq. (1) that and can be used as the control variables of the glide trajectory design.

2.3. Constraint Analysis

During the gliding phase of vehicle flight, due to the long-term flight in the thin atmosphere, energy dissipation will be caused due to aerodynamic forces, which will produce strong mechanical and thermal effects, such as heat flow, overload, and dynamic pressure. At the same time, to meet the wide range of maneuver of the vehicle within the atmosphere, quasiequilibrium gliding should also be considered in the design of the profile to ensure flight capabilities during the maneuver process [37]. Therefore, a variety of path constraints must be satisfied in the trajectory design of the vehicle to achieve the relevant state variables within the acceptable envelope range. In summary, the constraint conditions are as follows: where is a constant, is the radius of the vehicle nose, is the reference velocity, is the maximum heating rate, is the maximum load factor, and is the maximum dynamic pressure.

In addition, the trajectory design for the gliding phase also needs to satisfy the terminal constraints, including the terminal altitude , terminal velocity , and terminal remaining flight range .

3. Multiple Leap Maneuver Flight Profile Design

In this paper, to achieve the goal of an active leap maneuver within the longitudinal flight profile during the gliding phase, the standard maneuver reference drag acceleration curve is designed in the drag acceleration-normalized energy ratio (D-e) profile. First, the gliding corridor is determined in the angle of attack flight profile. Then, the reference standard maneuver drag acceleration is determined based on the analytic form of the single-parameter mapping of the range to go. Finally, the algorithm for updating and preprocessing the standard reference maneuver drag acceleration is provided.

3.1. Profile Design of the Angle of Attack

The gliding corridor needs to be determined to obtain the profile in advance. The design of needs to account for both the aerodynamic heating and range. In the initial stage of the gliding phase, the heating flow is the main constraint condition. At this time, a larger magnitude of should be adopted to reduce aerodynamic heating. When the heating flow is no longer the main constraint, should be reduced, and the range is increased near corresponding to the maximum lift-drag ratio. In engineering implementation, the envelope can be designed as a piecewise linear function of velocity [38] as follows: where and are the maximum angle of and the corresponding of the maximum lift-to-drag ratio, respectively, and and are two given velocities.

3.2. Determination of the Gliding Flight Corridor

During the gliding phase, the vehicle can achieve changes in the altitude by altering the magnitude of drag. Since the energy that can characterize the change in the flight state during the gliding phase decreases monotonically, the energy is a more appropriate trajectory design-independent variable than time. Therefore, in this paper, the glide trajectory is designed in the D-e profile. The energy of the vehicle during flight is expressed as follows:

To improve the computational efficiency, the energy is normalized as follows: where is the normalized energy of the vehicle, is the initial energy of the gliding phase, and is the terminal energy of the gliding phase.

Based on the determination of in Section 3.1, the gliding flight corridor of the drag acceleration-normalized energy (D-e) profile can be obtained according to the maximum dynamic pressure , maximum overload , maximum heating flow , and quasiequilibrium glide constraints given in Eq. (4) [39]: where is the geocentric distance, which can be taken as a function of the linear change with the normalized energy as follows: where is the initial geocentric distance of the gliding phase and is the terminal geocentric distance of the gliding phase.

Combining the above constraints, the upper and lower boundaries of the drag acceleration-normalized energy (D-e) profile can be obtained as follows:

and represent the upper and lower boundaries of the D-e profile gliding corridor, respectively.

3.3. Drag Acceleration Profile Design of Multiple Leap Maneuver Mode
3.3.1. Reference Drag Acceleration Profile Planning Algorithm

To enhance the longitudinal maneuver flexibility of the vehicle and achieve multiple leap maneuver, the form of the single-stage maneuver drag acceleration is designed as follows: where and are the normalized energy parameters corresponding to the design nodes of the initial and end drag acceleration, respectively. is the base drag acceleration corresponding to the design node . is the additional drag acceleration. The drag acceleration form of a single-segment maneuver is composed of an oblique straight segment, a horizontal segment, and an oblique straight segment. It changes the original drag acceleration through additional drag acceleration to meet the longitudinal maneuver of vehicle. The single-stage maneuver drag acceleration profile is shown in Figure 2.

A flexible maneuver configuration can be achieved by introducing multiple single-stage maneuvers during the gliding flight phase. The design nodes of the multistage maneuver drag acceleration profile can be expressed in the form of a set as follows:

In the design of multistage maneuver acceleration profiles, the base drag acceleration and the additional drag acceleration of the th maneuver are expressed as and . When , , and are known, the distribution design of the multisegment maneuver drag acceleration profile can be further carried out. In this paper, takes the three-segment drag acceleration profile form as the value basis as follows: where and are the drag acceleration of the initial and terminal of the gliding phase, respectively. is the drag acceleration that satisfies the constraint of the remaining flight range . The multistage maneuver process is limited to the horizontal segment of drag acceleration in Eq. (14). The elements in the design node set satisfy the following constraint:

The preset maneuver corridor is , and the drag acceleration is taken as , which is the benchmark. The estimated range to go based on the single-stage maneuver drag acceleration mode is calculated as follows:

According to the drag acceleration profile analytical predictive correction method, in the form of single-segment maneuver drag acceleration has analytical solution [40]:

The estimated range to go based on the multistage maneuver drag acceleration mode is calculated as follows:

In Eq. (17), a total of stages of maneuver drag acceleration are introduced in the design, and the preset range to go requirement is realized through a single-parameter iterative search [41]. In this paper, the Newton iterative method [42] is applied to solve for the profile design parameter as follows:

3.3.2. Reference Drag Acceleration Profile Update Algorithm

Due to the influence of the aerodynamic uncertainty during the gliding phase, the drag coefficient of the vehicle needs to be corrected. The corrected terminal drag acceleration at the gliding terminal time is as follows [43]: where is the drag acceleration coefficient deviation and is the predetermined drag acceleration at the gliding terminal time before updating.

In the process of gliding maneuver, due to the simplified calculation of the estimated range to go and the assumption of the large arc of the expected flight range, the drag acceleration profile needs to be updated to correct the range deviation in the actual flight process. In the multistage maneuver acceleration profile, based on the maneuver drag acceleration form designed in Eq. (17), the profile update nodes are chosen as follows to minimize the impact on the maneuver process: where is the update node additional item and satisfies .

The profile is updated in the form of three segment broken lines represented by Eq. (14). It is known that the remaining design nodes , , and are not updated. The normalized energy nodes of the updated drag acceleration profile are reselected as follows: where avoids large changes in drag acceleration after profile updating due to the small energy interval by increasing the scaling factor of the energy interval about the additional drag acceleration.

The updated reference drag acceleration profile maintains the form of Eq. (14).

3.3.3. Reference Drag Acceleration Profile Preprocessing Algorithm

Using a piecewise linear drag acceleration form, Eq. (14) may lead to first-order discontinuity in the design profile, which could lead to significant abrupt changes in the guidance command. Therefore, it is necessary to implement a smooth transition process near the segment endpoints or the profile design nodes to address this issue.

The set of preset values can be expressed as follows:

The smooth transition profile is expressed as follows: where is the design node and . The value of is small. The value of requires that the smoothed profile has little effect on the estimated range to go and satisfies the following constraints:

Eq. (24) ensures that the smooth transition trajectory does not exceed other design nodes and maintains the main characteristics of the broken line reference drag acceleration curve.

4. Maneuver Drag Acceleration Tracking Guidance Law

When the drag acceleration profile design is obtained, it is necessary to design the tracking guidance law to track the design profile. In this paper, the bank angle is chosen as the control variable to ensure that the actual drag acceleration equals the corresponding standard reference drag acceleration of the profile. In this section, the method based on the sliding mode and preset performance control is used to design the standard reference drag acceleration tracking guidance law to ensure that the vehicle satisfies multiple constraints at the gliding terminal time.

The error conversion function is defined as follows: where and are the convergence error envelopes. Moreover, and represent the actual error convergence boundaries, which are designed to be in an exponential convergence form as follows: where and are the initial and final values of the error envelope, respectively, and is the convergence rate factor.

In Eq. (25), is the sliding mode variable, which can be expressed as follows: where is a sliding mode design parameter greater than zero. The reference drag acceleration tracking error and its derivative are expressed as follows:

Reference [44] has demonstrated that when the sliding mode variable tends to zero, will converge to zero. Taking further derivative of the drag acceleration tracking error yields as follows: where is the virtual control variable, , and is the uncertainty. Additionally, and can be determined by the following equations:

The derivative of the conversion error defined in Eq. (25) can be obtained as follows: where , , and are expressed as follows:

Substituting Eq. (29) and Eq. (33) into Eq. (31) yields the following expression:

Considering the uncertain disturbance term, is designed as follows: where is the balance term, which is used to counteract the known term in Eq. (34). The construction form of is as follows:

By substituting Eq. (35) and Eq. (36) into Eq. (34), the following simplification can be obtained:

The uncertain disturbance is passed through a first-order low-pass filter to obtain the following equation: where is the time constant term and satisfies [45]. From the first-order low-pass filter, the error estimation can be obtained as follows:

The variable is designed as follows: where is the design parameter and satisfies [45]. By substituting Eq. (40) into Eq. (39), the following expression can be obtained:

Both sides of the above equation can be integrated to obtain the following expression: where is the initial value of estimation error. By substituting Eq. (42) into Eq. (39), can be simplified as follows:

According to Eq. (35) and considering the amplitude constraint, the magnitude of the control variable can be obtained as follows:

It is proved in Reference [45] and Reference [46] that the conversion error function is uniformly bounded by the guidance law (Eq. (44)) under the condition that the interference error is bounded, the sliding mode surface is guaranteed to be within the error envelope, and the tracking error approaches zero.

5. Simulation Verification and Analysis

The CAV-L vehicle model is taken as the simulation object. Its mass is 816.48 kg, and the reference area is 0.32258 m2. The lift-drag coefficient is provided in Reference [43]. The design parameters of the angle of attack profile are as follows: , , , and . The design parameters of the drag acceleration profile of the two-stage longitudinal leap maneuver are selected to be and . The design nodes are set as follows:

The guidance parameters are set to , , , , , and . The constraints for the dynamic pressure, load factor, and heat flux during the gliding phase are , , and , respectively. The target location is 115°E and 32°N.

5.1. Nominal Mission Simulation

To verify the effectiveness of the designed reference drag acceleration profile and the tracking guidance law under multiple active leap maneuver modes, two sets of different initial conditions are selected for simulation analysis in this section. The initial state and the terminal constraints of the gliding phase are shown in Tables 1 and 2, respectively.

Simulations were conducted using the drag acceleration profile design method and tracking guidance algorithm proposed in this article. The results are as follows. Figure 3 is the reference drag acceleration and actual drag acceleration curve under the two-stage longitudinal leap maneuver modes in the D-e profile. It can be seen from the figure that the turning point of the reference drag acceleration profile is smooth and the tracking error of the actual drag acceleration is very small. This indicates that the vehicle can effectively control drag acceleration and maintain a stable flight state during the gliding phase.

The variations in the vehicle angle of attack and command bank angle are shown in Figures 4 and 5, respectively. Figures 6 and 7 present the three-dimensional trajectories and velocity-altitude curve during the gliding phase, respectively. It can be concluded that the algorithm studied in this paper can guide the vehicle to the terminal expected state and satisfy the multiple constraint conditions. The guidance algorithms exhibit strong adaptability. Further analysis reveals that the variation in the reference drag acceleration designed in Figure 3 can achieve the longitudinal leap maneuver shown in Figure 6. When the drag acceleration significantly changes, the high-precision tracking of the drag acceleration can be realized by implementing the rapid large-scale bank angle.

Table 3 shows the terminal actual state of the gliding phase. The terminal altitude error is less than 220 m, the terminal range to go error is less than 100 m, and the velocity error is less than 10 m/s, which can accurately satisfy the constraints of the gliding handover.

5.2. Monte Carlo Simulation

To verify the robustness of the design method to the initial state deviation and aerodynamic parameter uncertainties, the Monte Carlo simulations are carried out. The initial state of the gliding phase simulation of the vehicle is based on the mission 1 state in Table 1, and the range of disturbance deviations is selected based on the principles, as shown in Table 4.

The Monte Carlo simulation results of 200 dispersed cases are shown in Figures 811.

From the calculation results of Figures 8 and 9, it can be seen that due to the influence of aerodynamic parameter deviations, there is a certain deviation between the actual drag acceleration and the reference resistance acceleration in the initial stage of gliding. The gliding guidance law studied in this paper can quickly converge to the design profile of the reference drag acceleration in finite time and maintain high-precision tracking. The drag acceleration state of the reference terminal is corrected by the online estimation of the aerodynamic coefficient to ensure that the vehicle meets the terminal state handover constraint requirements. The terminal positions during the gliding phase under comprehensive deviation conditions are shown in Figure 10. The actual positions are all located on the circle that satisfies the remaining gliding terminal range constraints. The statistical simulation results show that the mean and variance of are 19986.8 m and 75.3 m, respectively, and the mean and variance of are 18032.4 m and 169.6 m, respectively. The Monte Carlo simulation results demonstrate that the design method of the drag acceleration profile and the gliding guidance algorithm for the multiple leap maneuver mode proposed in this paper are reliable and effective, with strong adaptability.

6. Conclusions

In this paper, a fast design algorithm of drag accelerations that can be used to achieve multiple maneuvers in the gliding phase and a gliding guidance algorithm with a strong anti-interference ability are proposed. Through single-parameter analytical iterations, the reference drag acceleration profile of multiple leap maneuvers satisfying process and terminal constraints can be quickly obtained, which has the potential for online planning applications. The drag acceleration tracking algorithm based on the sliding mode and prescribed performance control has a strong anti-interference ability. The simulation results show that the drag acceleration profile design method and tracking guidance algorithm in this paper have good adaptability and robustness based on different missions and disturbance deviation conditions and have great engineering application potential.

Data Availability

The data that support the findings of this study are available from the corresponding author upon reasonable request.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Acknowledgments

This study was funded by the National Natural Science Foundation of China (grant no.: 62101579).