Research on Speed Scheme for Precise Attack of Miniature Loitering Munition

The loitering munitions are advanced new ammunitions, which reﬂect both the characteristics of UAV and missile and have their own novel characteristics. This paper analysed the characteristics of the attack process of a small loitering munition and proposed a speed scheme suitable for precise strike of loitering munitions based on the concept of traditional ammunition’s speed-thrust schemes. Then, the boundary conditions and the existence of the solution have been discussed. Finally, ﬂight test results showed that the scheme was eﬀective.


Introduction
A loitering munition (LM), as a new member of the precisionguided munitions [1][2][3][4][5], fits in the niche between cruise missiles and unmanned combat aerial vehicles (UCAVs) sharing characteristics with both. In the cruise section, the LM mainly conducts flight and mission target search in the form of a UAV, but in the launch section and the final guidanceattacking section, the LM works as a missile [5][6][7]. After launch, the original folded wings are unfolded, and the aircraft is converted from projectile form to drone form, and the configuration of the loitering munitions does not change again after the conversion of the machine finished [1,[8][9][10]. erefore, it is a challenge for a LM to adapt to the control requirements of different mission phases under the same aerodynamic shape. e typical working process for a loitering munition includes pneumatic launch, loitering and gathering intelligence, locking on target, aborting mission, or precise attack. Some mature products of small LM are shown in Figure 1, and the working process diagram for a LM can be seen in Figure 2 [11].
Traditional precision-guided munitions include various types of guided missile, guided bomb, guided projectile, and cruise missile [2,[12][13][14]. e design of the propulsion system and guidance control system occupies important positions in the system design, but they are kept independent as usual. e main goal of propulsion design is to achieve effective range and reach the final flight speed to meet the requirements of the attack mission [15]. For example, the antitank missile uses a rocket engine as the propulsion to fly at high speed, and its main speed schemes include the following three types: boost + glide flight, boost + acceleration cruising + glide flight, and acceleration boost + constant velocity (boost + cruising flight). A propulsion map and speed of typical high-speed kinetic energy antitank missile are shown in Figure 3, and in general, while designing the guidance algorithms ,it will be assumed that the flight speed is constant, especially at the terminal guidance phase [16][17][18], or the flight speed is time varying passively affected by gravity, propulsion, and aerodynamic forces [12,[19][20][21]. At the same time, the research topic of speed is mainly focused on the cruise flight process [22,23] or is considered as a terminal constraint at some special flight mission, for example, a rendezvous mission to a tanker [24].
Because the pneumatic layout and aerodynamic characteristics of LMs are closer to UAV than to missiles [25,26], it is unreasonable to assume that the speed is constant or ignore the dynamic adjustable characteristics of speed just like a traditional precision-guided munition in designing the guidance and control algorithms in the terminal guidance phase. erefore, it is of great significance to carry out research on the speed scheme at the terminal attack stage for a LM.
is paper is concerned with the problem of designing speed-trajectory characteristics of the LM within the attack phase. First of all, the speed requirement of guidance and control design in the LM's precision attacking process is discussed. en, a speed scheme is proposed to improve the accuracy of a LM strike. In the following sections, the boundary condition for the establishment of the scheme is    given, and the existence conditions of the system solution is analysed. Finally, a flight test has been implemented to verify the theoretical rationality and practical value of the proposed scheme, where a method for parameter acquisition was given at the same time.
It should be noted that the main topic of this paper is to analyse the dynamic characteristics of thrust and velocity, mainly from the theoretical analysis to solve the problems in practical applications, so the specific implementation of propulsion is not discussed.

Equations of Aircraft Motion.
As a novel precisionguided munition, LM in the form of small UAV could use the low-cost strap-down seeker [12,27,28] to achieve target capture and tracking during autonomous attack mission, so we discuss problems based on strap-down seeker against a stationary target. e flight model was expressed by differential equations in body-fixed coordinate which was defined as follows: the origin of coordinate system is located at the LM centroid. X axis is aligned with the longitudinal axis of LM; positive direction points to the LM's nose. Y axis is perpendicular to the longitudinal axis and lies in horizontal plane; positive direction points to right. Z axis is perpendicular to XY plane; positive direction points down. is body-fixed coordinate system is illustrated in Figure 4.
In order to introduce the speed scheme clearly, only the longitudinal motion of a LM has been discussed. e force analysis of a LM is shown in Figure 5, and the homing guidance geometry of missile-target engagement of the flight process is shown in Figure 6.
L is the lift force (aerodynamic force along z direction), D is the drag force (aerodynamic force along x direction), T is the thrust of engine, m is the mass of LM, g is the gravitational acceleration, and θ is the pitch angle (Euler angle along y axis) of LM.
R is the instantaneous separation between LM and target, V is the speed of the loitering munition in x direction in body-coordinate, V Γ is the speed of the loitering munition in z direction in body-coordinate, η is the angle made by the line of sight with respect to the fixed reference, and q is the pitch rate (angular rate along y axis).
Assume that V Γ is small and θ is close to η, and the thrust T is assumed to act through the LM centroid and along the x axis. en, from the engagement geometry and the force analysis of the aircraft, we can see the following. e relative motion relationship between the LM and the target can be described as the following equations: e dynamic equations of the LM can be written as follows: where I y is the LM's moment of inertia in y direction. e equations of aerodynamic force and moment can be written as follows: whereM is the aerodynamic moment component in y direction, C D is the drag coefficient, C L is the lift coefficient, C m is the pitching moment coefficient, Q is the dynamic pressure, S is the reference area, and c A is the mean aerodynamic chord.
In the above expressions, the pitching moment coefficient can be expressed as where C mL is the moment derivative with respect to lift, C mδ e is the moment derivative with respect to control-surface deflection, δ e is the control-surface deflection, and C mq is the moment derivative with respect to pitch rate. e dynamic pressure can be expressed as where ρ is the density of air. en, the equations describing the LM's motion for the attacking stage can be expressed as e state variables of the system during attacking can be selected as e control variables are Equation (6) can be organized as is is the system equation that describes the motion of the aircraft at the attacking phase.

A Speed Scheme Suitable for a LM with Strap-Down Seeker
for the Attacking Stage. For a LM with strap-down seeker, it is important to keep flight stabilized to get an ideal target tracking feedback, which means there should be no excessive normal shifting motion during attacking. e desired normal acceleration should be set: According to the 4th expression of equation (9), this condition can be described as It can be solved from the above equation that the expected speed during the attacking process is where the subscript "tc" indicates that the value of the variable is an expected value which depends on the attacking process.
When the LM is flying at this expected speed, it means According to the 3rd expression of equation (9), this condition can be described as By substituting equation (12) into equation (14), the equilibrium thrust required should be A speed scheme that satisfies strap-down LM precision attack requirements could be defined by equations (14) and (15), and it can be formed as

Boundary Conditions for the Solution.
From equation (16), we can see that the parameters affecting the equations can be divided into two parts. We named them the static parameters and the dynamic parameters. e static parameters mainly include the gravitational acceleration, LM's weight, the density of air, and the wing reference area. e dynamic parameters mainly include lift coefficient, drag coefficient, pitch rate, and pitch angle. In this paper, the static parameters were mainly regarded as known conditions and only the dynamic parameters have been analysed.

(1) Pitch Angle Limited Value Ensures at the rust Has a
Solution. e solution of the speed scheme requires where T max is the maximum thrust of the engine. It can be seen from equation (16) that V tc > 0 is obvious; for T tc , there is e design of the propulsion system should naturally meet the upper boundary requirements, so this paper studied the lower boundary of the pitch angle only.
According to (18), when given a certain drag coefficient, there is at means where the subscript "lb" means this pitch angle is the lower boundary for speed scheme to have a positive thrust. erefore, in order to meet the requirements of this scheme, the relationship between the limit value of the pitch angle and the flight speed during the attacking must obey the expression above, and it can be represented by a curve, as shown in Figure 7: According to equation (20), for different LMs or just for different flight status, if the drag coefficient CD is changing, the above curve in Figure 7 should be turned to a surface, which satisfies the relationship between the pitch angle value and the drag coefficient CD, as shown in Figure 8.
It can be seen from Figure 8 that while the speed is constant, the greater the drag coefficient the steeper the pitch angle, and while the drag coefficient is constant, the greater the speed the steeper the pitch angle.
(2) Balance between Flight Speed and Pitch Angle. While the LM's lift coefficient is constant, it can be seen from the scheme equation (16) that the larger the pitch angle (absolute value) is, the smaller the flight speed that satisfies the normal equilibrium condition. According to the 1 st expression of equation (16), a balance pitch can be calculated by equation (21), and the balance map of flight speed and pitch angle is shown in Figure 9.
where the subscript "balance" means this value is specific for the variable (pitch angle) to meet the speed scheme's expected speed V tc .

Mathematical Problems in Engineering
In addition, when CL changes within a certain range with flight conditions, the relationship between balance pitch angle, CL, and speed is shown in Figure 10: As can be seen from the above figure, in order to obtain a large attacking speed and obtain a steep attacking angle, the LM should have a small lift coefficient.
(3) Valid Pitch Angles according to Prerequisites of Lower Boundary and the Balance Relationship. For a given set of aerodynamic characteristics [CD, CL], it can be known from the speed scheme that the pitch angle during attacking must satisfy both the lower boundary and the balance relationship we have given. We named these pitch angles "valid pitch angles." is requirement is described in Figure 11.
When the aerodynamic characteristics [CD, CL] change within a certain range, valid pitch angles that satisfy the solution of the speed scheme under different conditions can be obtained, as shown in Figure 12.
erefore, the valid pitch angles during attacking phase with respect to certain aerodynamic parameters can be expressed as where the subscript "valid" means the pitch angle satisfies both the lower boundary and the balance relationship. According to equation (22), if the lift coefficient and the drag coefficient of a LM during attacking are known, we can calculate the boundary condition and the balance points of the pitch angle and then set the attack conditions according to a valid pitch angle. At the same time, the thrust needed can be calculated according to the speed scheme.

e Existence and Uniqueness of Solution to the Motion of the LM.
Having determined the boundary conditions of the solution, we only know that the scheme has the application basis with the corresponding conditions. In this section, the LM's motion reliability during the attacking flight was confirmed. We attributed this problem to a uniqueness problem of the system solution [29].
When the design of the control system is completed, the ideal attitude control loop can be equivalent to an inertia unit with a time constant of T θ . Assume that the flight process is stable and the guidance process is implemented by the tracking guidance law, that is, the pitch rate is 0 (q≈0). Combine the proposed speed scheme equation (16) with the LM motion equation (9), and a new description of the system can be given: System (23) can be linearized at any working point and expressed as follows: Based on Taylor's expansion, the system characteristic matrix can be given as follows: e index of the variable "e" identifies the value at the working point, for example, flying according to the speed scheme is a special case.
It can be seen from the system characteristic matrix that the elements in the fourth column are all zero. at means the normal velocity V Γ does not affect the condition of the 6 Mathematical Problems in Engineering  system. So, in order to simplify the analysis of the system, the influence of the normal velocity can be ignored, and a simplified system characteristic matrix can be given: So, we get According to the compression mapping principle [30], when the system parameters meet the condition the system should have a unique solution.

Scheme Verification
After analysing the boundary conditions and the uniqueness of the solution, requirements of the program application have been grasped basically. So, we carry out the verification work next. First of all, the steps of the test were introduced. In this work, we found that there were many different methods to obtain the parameters the test needs, and different methods are different in difficulty and accuracy. erefore, a method to get the lift coefficient and drag coefficient was introduced first, and then the experimental flow of scheme verification was given.

Acquisition of Parameters and the Scheme Testing Process.
e thrust was calculated according to the equation (16), wherein the LM weight m, the density of air, the wing reference area, and the gravitational acceleration can be regarded as known parameters. e pitch angle and the pitch rate could be measured with sensors during the flight. erefore, the lift coefficient and the drag coefficient were the key parameters for the realization of the speed and thrust scheme, which needs to be acquired for the verification of flight test. e following is a further explanation of the acquisition of the lift coefficient and the drag coefficient before the verification process.
For small-sized aircraft, obtaining aerodynamics data through wind tunnel test has the disadvantages of complicated experiment and high cost. Especially for occasions requiring rapid prototyping and special control performance verification, traditional wind tunnel test cannot meet low cost and rapid evaluation in engineering application. e commonly used methods also include the aerodynamic parameter identification method based on flight test and the means of simulation calculation [31][32][33][34]. e main disadvantage of the simulation calculation is that the accuracy cannot be reliably guaranteed, and the result needs to be evaluated and revised using a series of experimental data.
In order to reduce the cost of experiments and quickly obtain parameters closer to the actual situation than simulation calculations and in order to ensure the safety of the test, this paper adopted the aerodynamic parameter identification method of the aircraft balance point based on an object understanding calculation method.
According to equation (9), the formulas we introduced to calculate the lift coefficient and the drag coefficient are C L � V fly data q fly data − _ V Γfly data + g cos θ fly data (1/2)ρV 2 fly data (S/m) , where the subscript "fly_data" indicates that the specific values should be obtained from a fly test. In order to avoid the situation that the singular matrix cannot be solved and in order to obtain a more accurate solution using experimental data with measurement errors, the above equations were extended to overdetermined equations.
Rewrite the above expressions as where 8

Mathematical Problems in Engineering
V fly data (2)q fly data (2) − _ V Γfly data (2) + g cos θ fly data (2) · · · V fly data (i)q fly data (i) − _ V Γfly data (i) + g cos θ fly data (i) In the above expressions, "i" is the number of equations in the overdetermined equations. Because solving the overdetermined equation is not the focus of this article, it will not be described in detail [35,36]. e parameter calculation and scheme test flow are given in Figure 13: e following is a detailed description of the verification flight test process.

Prophase Flight Test for Data Collection.
e prophase flight test was designed to collect data for parameter calculations. e UAV used for flight test is shown in Figure 14, and the parameters of the UAV are shown in Table 1.
A simple data acquisition and flight control module was designed. e hardware system is shown in Figures 15 and  16, and the device list is shown in Table 2.
According to the equations (32), (34), (35), and (37), the data needed to be collected include speed, pitch angle, thrust of the engine (in the experiment, the output percentage of the engine controller was collected and then converted into the actual thrust), pitch rate, and accelerations in x direction and z direction in body-coordinate frame. e data storage frequency was 50 Hz and the duration was 2 s. e obtained experimental data are shown in the following figures: Figure 17 shows the output percentage of the engine controller, Figure 18 shows the acceleration in z direction in body-coordinate, Figure 19 shows the acceleration in x direction in body-coordinate, Figure 20 shows the pitch rate, Figure 21 shows the pitch angle, and Figure 22 shows the flight speed.

Calculation of Lift and Drag Coefficients.
e lift coefficient and the drag coefficient were calculated according to equations (30) and (31) after the experimental data were obtained, and the calculation results are shown in Table 3.

Boundary Condition Calculation.
According to equations (20)- (22) and the results of the previous section (see Table 3), the valid pitch angles could be calculated. e valid pitch angle map and the critical point of the valid pitch angle can be seen in Figure 23.
It can be seen that when the pitch angle is in the range of θ valid ∈ [−25 ∘ , 0], there will be a valid solution of the system according to the scheme proposed in this paper, and the attacking mission could be performed.
According to the previous flight test experience, the target attack range was 100-150 m after the target was clear in the field of view limited by the camera performance (the relationship between flight altitude, camera performance, and pitch angle is complex and is not the focus of this article, so it will not be described in detail). So, the reasonable range of attack altitude should be controlled within the following range: Mathematical Problems in Engineering where the subscript "attack" means the value is at the time the attack started. at is,

Judgement of the Existence and Uniqueness of Solutions to the Motion.
According to the experimental data analysis, the valid pitch angle was in the range of θ valid ∈ [−25 ∘ , 0], the flight speed was in the range of V ∈ [0, 23], and the pitch angle response time was T θ � 0.2 s. According to equation (27), sin min θ valid .
(40) So, the system will have a unique solution when at is, It can be seen from the above calculation result that there would be a unique state trajectory in the attacking process before the distance to the target is less than 2.2819 meters.
is limit distance satisfied the conditions of the test; in other words, the speed scheme proposed in this paper was feasible.

Control System Block Build according to the Speed Scheme.
e logic of the flight test: use the scheme to calculate an expected thrust during the attacking phase and cooperate with the attitude control module to complete the impact task against a nonmoving target.
e control system constructed in this paper consisted of an attitude control module and a speed control module. e attitude control module was realized by a classical PID controller. e speed control module was formed according to the speed scheme. e system block diagram is shown in Figure 24. e speed controller parameters were calculated according to equation (16) and are shown in Table 4.

Test Implementation Method.
e UAV was the same one used in data acquisition (see Figure 14), and the flight process can be seen in Figure 25 and can be described as follows: (1) Take Off. e UAV takes off by remote controller and flies to the vicinity of the target.   controller designed in this paper starts working automatically until the UAV hits the target.

Flight Test Result.
In this section, the result of the flight test is presented. e final 7 s flight path was selected for presentation purposes during attacking stages of the flight test, and the data are shown in the following figures: Figure 26 shows the thrust, Figure 27 shows the speed, Figure 28 shows the pitch angle, Figure 29 shows the acceleration in z direction in body-coordinate, Figure 30 shows the flight path, Figure 31 shows the first perspective video screenshot, and Figure 32 shows the third perspective video screenshot. Based on the experimental results, following conclusions can be given: First of all, Figures 26 and 27 show that the thrust and speed were dynamic curves, which were very different from the traditional ammunition speed throttle scheme (as can been seen in Figure 3), and their changes were all within reasonable ranges.
Secondly, as can been seen in Table 5, the mean of acceleration in z direction from the scheme flight test was −0.6351 and that from remote control test was −9.323. e standard deviation of acceleration in z direction from the scheme flight test was 1.0596 and that from remote control test was 1.727. Compared with data from the remote control test (which has been shown in Figure 18), acceleration in z direction in body frame has been suppressed clearly (as can been seen in Figure 28).
At the same time, the pitch angle and the flight trajectory curve is smooth (as shown in Figures 29 and 30), which    Acceleration, x Figure 19: Acceleration in x direction in body-coordinate.
strongly verifies that the speed scheme proposed in this paper was helpful to improve the attitude and trajectory stability of the UAV during attacking. Finally, as can been seen in the screenshots of the flight videos (Figures 31 and 32), the UAV successfully hit the target, which further illustrated the effectiveness of the scheme.
3.6.4. Defect Analysis. However, the airspeed during the flight did not follow the expected trajectory perfectly, as shown in Figure 33. Further analysis is given below.  (1) Dynamic Response Characteristics of the rottle. Since the airspeed did not change much during the whole flight process, it can be assumed that the thrust depended only on the working state of the propulsion system and was independent of the atmospheric conditions. According to equation (9), using the thrust as input,     the transfer function of the speed output can be expressed as (43) Assuming that the throttle thrust input signal is stepsignal ΔT1(t), based on final value theorem of Laplace transform, the speed response is (44) is shows that the change of the throttle can only adjust the speed during the transition process but cannot fully achieve the purpose of controlling the speed according to the ideal trajectory.
(2) Stability Analysis of Speed. From equation (24) and the system characteristic matrix equation (25), the differential expression describing the airspeed motion process at any operating point after linearization can be given by When the system is aiming at the target, there is Its characteristic roots are where ρC D SV e m > 0, is means the two eigenvalues of the differential equation have negative real parts, so the change of the airspeed trajectory should satisfy the stability condition, and    the test data also illustrated this feature: the speed curve rises following the expected one ( Figure 33). However, as shown in Figure 23 (the control system block diagram), only an open-loop controller was used to control the speed to achieve the desired value V tc . Obviously, an open-loop controller's performance cannot meet the requirement of a complex testing environment, and the experimental results also verified this conclusion: the speed curve rises but slowly.
It would be nice to design a better controller to improve the response characteristics of the speed, but how to improve the control law requires further analysis, and it is out of range of this research, so it is one of the tasks we will continue to expand in the future.

Conclusions
In this paper, a speed scheme for small loitering munitions has been given for precision strike mission. e pitch angle boundary of the solution during attack was given, and the existence and unique conditions of the system solution were analysed. en, a flight test was taken to verify the scheme, and the result showed that the scheme was effective. At last, based on the analysis of experimental data, this paper summarizes the shortcomings of the program and proposes a clear direction for further research.

L:
Lift force, N D: Drag force, N θ: Pitch angle, rad T: rust, N m: Mass of the loitering munition, kg g: Gravitational acceleration, m/s 2 q: Pitch rate, rad/s η: e angle made by the line of sight with respect to the fixed reference, rad V: e speed of the loitering munition in x direction in body-coordinate, m/s R: e instantaneous separation between LM and target, m V Γ : e speed of the loitering munition in z direction in body-coordinate, m/s I y : e LM's moment of inertia in y direction, kg·m 2 M: Aerodynamic moment component in y direction in fixed plane coordinate, Nm C D : Drag coefficient Q: Dynamic pressure, Pa S: Reference area, m 2 C L : Lift coefficient C m : (Pitching) moment coefficient c A : Mean aerodynamic chord, m C mL : Moment derivative with respect to lift C mδ e : Moment derivative with respect to control-surface deflection δ e : Control-surface deflection, rad C mq : Moment derivative with respect to pitch rate ρ: Density of air, kg/m 3 V tc : Equilibrium speed T tc : Equilibrium thrust T max : Maximum thrust θ lb : Lower boundary of pitch angle θ balance : Balance pitch angle according to speed scheme θ valid : Valid pitch angle according to speed scheme h: Flying height.

Data Availability
e data used to support the findings of this study are available from the corresponding author upon request.

Conflicts of Interest
e authors declare that they have no conflicts of interest regarding the publication of this paper.