Terminal Guidance with Attack Time and Impact Angle Constraints for Hypersonic Vehicle in Dive Phase

The terminal guidance problem of hypersonic gliding missiles under multiple constraints was considered regarding a coordinated attack in time and space. In this paper, a three-dimensional terminal guidance method based on the real-time updating of parameters is proposed with consideration to the attack time and terminal angle. By this method, the three-dimensional motion was decomposed into the horizontal and vertical plane. In each plane, a real-time solution of the optimal trajectory yielded angle commands that satisfied the terminal angle constraints. Additionally, the time was controlled by changing the flight path angle in the front section of the vertical plane. The results obtained by numerical simulation indicated the good applicability and robustness of this method.


Introduction
The cooperative combat of missile clusters has received much attention in recent years because it achieves good robustness and concealment.It can effectively improve the missiles' penetration capability.Furthermore, hypersonic gliding missiles have the advantages of long range, fast flight speed, large lift-drag ratio, and high maneuverability.To take full advantage of the respective advantages of the hypersonic gliders and clusters and to realize synergistic attacks with gliding missile clusters, it is necessary to solve the terminal guidance problems under multiple constraints in the same way as in [1].
In terminal guidance that considers multiterminal constraints, the guidance law, which considers the attack time and terminal angle constraints, plays an important role in the coordinated guidance of guided weapons and has achieved certain results.Various studies have been conducted on guidance with consideration to the impact angle [2][3][4][5][6].In this paper, related studies with regard to the attack time and impact angle will be discussed.For example, a twodimensional guidance law, which can satisfy the expected time and terminal angle, has been proposed by Si and Song to address the problem of UAV homing [7].This study investigated the feedback loop of the angle and the additional instruction of the time constraint.Lee et al. designed proper path points and angles by implementing a geometric method [8] and designed the guidance law of the terminal angle constraint.In combination with the attack time guidance law based on proportional guidance, the design of the terminal guidance law satisfies the angle and time constraints [8].Jung and Kim proposed a two-dimensional guidance law that controls the attack time and terminal angle [9].The landing angle requirement was satisfied when the drop point precision was zero.Then, the attack time was controlled by the additional control instruction with the feedback of the impact time error.Jeon et al. used the maximum principle to deduce the optimal guidance law with the shortest time as the performance index and strike the fixed target at a certain angle [10].Based on the Lyapunov stability theory, a three-dimensional missile guidance law with attack angle control was designed by Song and Shin [11].The time error was compensated by maneuvering the flight to realize time control.The abovementioned terminal guidance methods considered multiterminal constraints that assumed two-dimensional plane motion, constant velocity, and disregarding the aerodynamic drag.However, it is difficult to satisfy such assumptions during the actual flight process.
This paper presents a three-dimensional terminal guidance law based on real-time parameter updating.The proposed method decomposes the three-dimensional motion to the horizontal and vertical planes and uses the optimal control method to obtain the optimal trajectory in real time by combining the current motion state and terminal constraint information.The obtained angle instructions satisfied the terminal constraints, and the process error was eliminated through the real-time parameter updating.Moreover, we propose to control the attack time by changing the flight path angle in the forepart of the vertical plane.The effectiveness of the method was verified by the results of numerical simulation.Additionally, the results revealed the high precision of the landing point, terminal angle, and attack time.

Problem Formulation
The earth fixed coordinate system O_xyz is shown in Figure 1.The origin of the coordinate system is the target point, which is approximately still.The X axis points refer to the east, while the Y axis refers to the north.The Z axis represents the vertical upward direction in the right-hand coordinate system.Various relevant angles are defined below.The azimuth angle is measured from the positive X axis to the projection of the missile target line in the horizontal plane, and its range is −π ≤ θ ≤ π.The elevation angle is measured from the horizontal plane to the missile target line.When the missile target line is higher than the horizontal plane, the elevation angle is positive.The path angle γ is the angle between the velocity variable and the horizontal plane.This angle is positive when the velocity variable is higher.The heading angle ψ is measured from the positive Y axis in the clockwise direction to the horizontal projection.
In the above O_xyz coordinate system, the threedimensional motion equation is expressed by (1).The terms D and L represent the aerodynamic drag and lift, respectively.Both of them are related to the angle of attack, height, and velocity.The terms including x, y, z represent the centroid of the gliding missile, and σ represents the bank angle.
Guidance is designed to determine the appropriate attack angle and bank angle curves, such that the aircraft can strike the target according to the guidance law.Once these two commanded angle curves are designed, the entire movement of the aircraft can be determined.
In O_xyz, the ideal placement, attack time, and terminal angle constraints can be expressed as follows: The subscript f represents the value of the variable at the terminal.

Design of Guidance Law
By decomposing the three-dimensional motion into two planes, we can get the horizontal and vertical trajectory equations as (3) and (4).Additionally, s is the distance travelled in the horizontal plane.Therefore, s f = 0.
3.1.Horizontal Plane.In [6], the method of maneuvering in the horizontal plane was used to control the attack time.International Journal of Aerospace Engineering However, for the hypersonic gliding missile, its high altitude and small air density resulted in the weak control ability of the path angle.If the missiles were maneuvered in the horizontal plane, this would increase the pressure of the path angle control and could also increase the drop point error and terminal angle error.Taking this possibility into account, the time control of the hypersonic gliding missile in the form of lateral maneuvering was not considered in this study.The guidance law satisfying the terminal angle constraint was only considered in the horizontal plane.From (3), we can obtain the following relationship: Let where x i ′ represents x i , which is derived from y.Then, (3) can be rewritten as follows: x 20 = tan ψ 0 , This is an optimal control problem.By considering that y is an independent variable, the problem can be treated as the mathematical design of a curve that satisfies the endpoint's position and slope, as shown in Figure 2.
In the above optimal control problem, the relationship between the control variable u and the commanded heading angle can be expressed as follows: Owing to the limitation of the aerodynamic force, the size of the admissible control u varies with height and speed.Therefore, it is difficult to determine the scope of the admissible control and hard to obtain an analytical solution.By ignoring the size of u, the performance index is expressed by (9).Then, the optimal control problem can be solved by the variational method.
The solution of the above optimal control problem is expressed as follows: where a 1 and a 2 are related to the state and constraint of the terminal point and can be solved as follows: Guidance Law of Terminal Angle Constraints.The trajectory equation in the vertical plane can also be expressed as follows: where z i ′ denotes z i for the derivation of y.This problem could be solved by the optimal control method.However, the difference is that in the vertical plane the last stage is where the angle change is mostly concentrated.In the terminal period of the guidance, its height decreases drastically, which results in greater air density.Therefore, the angle variation ability is stronger in the low area.We may assume that the control variable is expressed as follows: where the constant term k ensures that there is always a certain angle change tendency at high altitude and that the The control variable u 2 is given in (13) without considering the performance index.If the above equation is solvable, the terminal constraints can be satisfied.By considering the influence of the actual height, density, and other changes, we designed the control variable directly as expressed by (13).
Let us consider parameter k as a known value and substitute the related initial conditions and terminal constraints.Then, parameters b 1 and b 2 in the upper equation can be solved by the integral.
The value of the constant parameter k can be approximated as follows: where L is the lift determined by height, speed, and mass, and η is a proportional coefficient, which is less than 1.
The relationship between the above control variable u 2 and the commanded angle can be obtained as follows: The above method solves the relative instruction angle through the present state and terminal constraints of the aircraft.Then, the attack angle a and bank angle σ can be obtained as described in [7].The next period state of the aircraft can be obtained by (1) such that the parameters are updated repeatedly.Each step is optimal under the condition of the control variable, which eliminates various errors that occur during the process.

Correction of Guidance Law with Consideration to
Attack Time.As shown in Figure 3 and according to (13), the longitudinal trajectory is almost straight in the early stages of guidance.The rapid dive is only realized in the last stage.During this period, the main change of the path angle is completed.In the early stage of guidance, the path angle changes little, the height descends slowly, and the aerodynamic resistance is small.Therefore, the speed remains almost unchanged.According to the trajectory equation in the longitudinal plane, the following relationship holds: In (17), the term v represents the effective average speed in the z direction, and the term T represents the total time.We can see that the key to control the attack time is to limit the size of the effective average speed v.
From (18) and ( 19), the following conclusions can be drawn: (1) The average velocity variation caused by the variation of the path angle γ is large.It is assumed that the path angle γ decreases continuously during the entire guidance process, and that the initial speed is 2800 m/s.If the initial path angle is reduced from 0 to −0.1 °, the effective average velocity v will increase by approximately 48.87 m/s (2) The increase of the drag caused by the changes of the path angle γ is small.It is noted that as the path angle decreases, the decreasing rate of the height accelerates, and the mean drag increases.In fact, because of the low atmospheric density and low resistance, the increase of the actual resistance is limited.
According to (18), the resistance assigned to the z direction is multiplied by the sin γ ratio such that its influence is further reduced (3) The change of the path angle γ also results in the increase of the gravity acceleration components.Equation (19) expresses that the acceleration effect of the gravity component increases with the decrease of the path angle γ.This further reduces the deceleration influence of the resistance.The acceleration assigned to the z direction is also multiplied by the sin γ ratio Thus, according to such conclusions, the path angle γ in the front section of the vertical plane has an important role 4 International Journal of Aerospace Engineering in time control.It can significantly increase the effective average speed and reduce the attack time.This method is more effective in comparison with the time control of the maneuvering method in the horizontal plane.Moreover, this method does not only aim to decrease the load of the direction deflection but also to increase the control capability of path angle γ through height reduction.It is possible to design the path angle to control the time effectively.Assume that the changes of path angle γ discussed in this section are as follows:

0, others, 20
where i is the adjustable parameter of the path angle.In this case, a height reduction strategy is adopted.The path angle γ is adjusted with parameter i when it is less than zero and its absolute value is smaller than the current elevation angle ϕ.
In any other case, the variability of γ is equal to zero.The value 10 −5 is the order of the experience value of the path angle control ability at the current time.
For this adjustment, we need to adjust the aircraft such that γ > 0 and γ < 0 before the guidance starts.To formally start the guidance process, we can use the following strategies to adjust the missiles: Step 2. The maximum lift L m of the current aircraft can be calculated.Let σ = 0, L = L m .
Step 3. Repeat the above process until the aerodynamic capability satisfies γ = 0, and then start the guidance.
The attack time varies with parameter i.Here, time prediction is adopted to solve i by dichotomy.The time prediction method is a simulation of the subsequent flight process.As shown in Figure 4, the main process takes into account factors such as the renewal of parameters during the integral process and the maximum lift limit.However, the limit of the maximum angle variation and the intermediate process error are not considered.
In Figure 5, the virtual line is the predicted time result, while the real line is the actual flight time of the glider.It can be seen that the curve of the forecasted result is close to the actual value curve and that the precision is high.In other words, the result can provide good support to the selection of parameters.At the same time, we can see that the actual flight time increases when parameter i increases.Obviously, there exists an upper time limit.
As shown in Figure 4, the flow of the entire guidance process was determined.With the exception of parameter i, which must be determined by the attack time, the other parameters can be calculated by the missile's current state of motion.

Simulation Analysis
Figure 4 shows the basic flow of the entire guidance process.During the simulation process, the aerodynamic data of the aircraft were based on CAV-H.The guidance step time was 0.1 s.Additional relevant data are listed in Tables 1 and 2.
Let the angle between the initial heading angle and the projectile connection be Δψ 0 .The angle between the terminal heading angle and the initial moment of the projectile connection is Δψ f .The relevant equations are expressed as follows: where θ 0 is the azimuth of the glider at the initial moment.
Based on the data listed in the above tables, the main performance of the method was investigated without considering the time constraints and by changing one or two initial parameters.Figure 6 shows the variation of the dropping error.The dotted line in the graph represents the dropping error that changed with Δψ 0 when Δψ f = −3 °.The error was smaller in the range of −18 °< Δψ 0 < 18 °.Beyond this range, the error increased sharply.The real line represents the terminal error that changed with Δψ f when Δψ 0 = −3 °.This error was smaller in the range of −26 °< Δψ f < 18 °.Therefore, this method should generally be used when the difference between the initial heading angle and the terminal heading angle does not exceed 20 °.Moreover, this method demonstrates that the hypersonic gliding aircraft has the characteristic of weak lateral maneuvering, owing to the high altitude.
The error decreasing with terminal path angle γ f is shown in Figure 7.When γ f ≥ −76 °, the guidance is precise.When γ f < −76 °, the error increases owing to the lack of aerodynamic capability.In other words, the guidance method is not suitable for a terminal path angle constraint of less than −76 As shown in Figure 8, the dropping error changes the time control parameter i.When the parameter i > −68, the error is small.If the required time is further reduced, that is, when i is less than 68, the error will increase remarkably.
Figure 9 shows the velocity-time curve with different i parameter values.From the end of the curve, it can be seen that the decrease of the speed curve slows down, while the effective average speed v increases.Therefore, the attack time is less, which agrees with the results obtained by the previous analysis.
The entire process was simulated without changing the parameters listed in Tables 1 and 2. As shown in Figure 10, the three-dimensional trajectory is in the top of the graph, and the projection in the horizontal plane inclines to the straight line.

6
International Journal of Aerospace Engineering The attack curve of the angle a is shown in Figure 11.Before the guidance began, the glider flew at a high attack angle, which allowed sufficient lift to maintain the altitude.Then, the guidance process started.As the height decreased, the atmospheric density increased.At this time, only a small attack angle was required to satisfy the lift requirement.Therefore, the attack angle gradually decreased.At the end of the guidance, the path angle had to be reduced sharply, and the required lift was negative.Therefore, the attack angle became negative at a certain rate.
Figure 12 shows the change in the bank angle over time.The bank angle was reduced to zero at a certain angular rate before the guidance was enabled.After the guidance began, the bank angle experienced a drastic change, which was related to the change of the attack angle from positive to negative and to the lift change direction.Therefore, a sign change was required for the bank angle.
The simulation result is presented in Table 3.The placement accuracy, terminal angle, and attack time reached a good level.The position error was within 0.01 m; the angle error was within 0.1 °, and the time error was zero.The method was effective and satisfied the constraints of landing accuracy, terminal angle, and attack time.

Robustness Analysis
The guidance process of gliding missiles is affected by many types of disturbance factors, with the main ones being as follows: starting point state error, atmospheric model error, aerodynamic coefficient error, aircraft quality error, observation information error, and aircraft structural design error.The main factors considered in this section are the following: atmospheric model error, observation information error, and aerodynamic coefficient error.Relevant information regarding the error types and sizes is presented in Table 4.All of the error types listed in the table had a Gaussian distribution.
To verify the robustness of this method, 1000 sets of random errors in the simulation were selected for analysis according to the size and distribution of the error presented in Table 4.    15 and 16.It can be seen that the heading angle had good accuracy and that its error was less than 0.1 °.The error of the path angle was slightly larger but lower than 1 °.The dropping error is shown in Figure 17, where it can be seen that, in most cases, the error was less than 5 m, and the end precision was high.The error distribution of the attack time is shown in Figure 18, where it can be seen    International Journal of Aerospace Engineering that the attack time error was approximately 0.6 s, which indicates that the time was controlled effectively.

Conclusion
This paper presented a terminal guidance method based on real-time parameter updating.The proposed method satisfies the terminal angle and attack time constraints.The performance and accuracy of the method were investigated through simulation analysis.This method decomposes the threedimensional motion of the aircraft into horizontal and vertical planes.And it designs the instantaneous orbit based on the current state of the aircraft, terminal constraints, and flight characteristics.Finally, the proposed method obtains the aircraft's control instructions.Time is controlled by changing the flight angle in the vertical plane.The simulation results indicated that this method could solve the terminal guidance problem with regard to the terminal angle and attack time.The results revealed that the placement, terminal angle, and flight time had good accuracy and that the method had good robustness.

Figure 1 :
Figure 1: Coordinate system and related angles.

Figure 3 :
Figure 3: Division of vertical guidance stage.

Figure 4 :
Figure 4: Terminal guidance process based on real-time parameter updating.

Figure 5 :igure 7 :Figure 8 :
Figure 5: Predicted time and actual time with parameter i

Figure 9 :
Figure 9: Different velocity curves with i.

Table 2 :
Initial state of missile.

Table 4 :
Type and size of error distribution.