Research Article  Open Access
Qiushi Zheng, Zhiming Zhou, "Flight Stability of CanardGuided DualSpin Projectiles with Angular Rate Loops", International Journal of Aerospace Engineering, vol. 2020, Article ID 2705175, 9 pages, 2020. https://doi.org/10.1155/2020/2705175
Flight Stability of CanardGuided DualSpin Projectiles with Angular Rate Loops
Abstract
Generally, as a precisionguided weapon, the missile has many disadvantage such as high price, difficult maintenance, and low yield. Modern war requires more and more new guided ammunition with high precision, low cost, and low collateral damage. Therefore, as a simple guided conventional ammunition technology, the dualspin projectile has attracted the attention of ammunition experts recently. This paper proposes a dualspin projectile scheme based on the rotation control method. Firstly, the concept of the dualspin projectile is introduced. Secondly, the mathematical model of the dualspin projectile is established, and the angular motion equation is finally obtained by using some linearized assumptions. Finally, the sufficient and necessary conditions of coning motion stability for dualspin projectile with angular rate loops are analytically derived and further verified by numerical simulations. It is noticed that the upper bound of the control gain is affected by the delay angle of the control system and the spinning rate of the projectile.
1. Introduction
The trajectory correction projectile is modified from the conventional ammunition, equipped with the ballistic detection device and the correction module. The ballistic detection device measures the trajectory elements and calculates the trajectory deviation in real time during the whole flight. Then, the projectile is guided by canards to reduce the impact miss distance and improve the damage probability. For the twodimensional rajectory correction, the spin rate of the projectile is a major constraint to the correction effect. Conventional ammunitions always keep a very high spinning rate after launching to maintain the fly gyroscope stability. Thus, it is difficult to use canards to guide the projectile because general actuators do not have such a large operating frequency. In order to solve this problem, a dualspin configuration principle has been proposed.
The concept of lowcost trajectory correction ammunition first proposed by the United States in the 1970s adopted the dualspin structure in [1]. The socalled dualspin structure means that the projectile can be divided into two decoupled parts along the axial direction, which are connected by roller bearings. In the following years, many scholars have studied the dualspin projectile. Costello and Peterson [2] developed the equations of motion for a dualspin projectile in atmospheric flight. Guan and Yi [3] established a sevendegreeoffreedom projectile trajectory model for the dualspin projectile and obtained a more accurate and reliable measurement data by filtering the trajectory of the dualspin projectile using the UKF algorithm. Burchett et al. [4] obtained the closed form expressions for swerving motion under the action of lateral pulse jets for a dualspin projectile in atmospheric flight by using the linear theory. Chang et al. [5, 6] studied the spinrate characteristic of the dualspin projectile and the dynamic response of the dualspin projectile under the canard control effect. Liu et al. [7] deduced the analytical expressions without attitude information, which can predict the increment orientation for trajectory angular rate and swerving motion induced by canard control force. Wernert et al. [8] carried out the wind tunnel test of a canards guided 155 mm artillery and got the result of openloop trajectory simulations.
There has been a considerable literature on the flight stability of conventional spinstabilized projectiles, such as [9–12]. But for this kind of dualspin projectiles, there are some differences. Costello and Peterson [2] analysed the stability characteristic of the dualspin projectile using modified projectile linear theory and obtained the expression of the gyroscopic stability factor and the dynamic stability factor which take on the same form as the rigid projectile. Wernert [13] addressed the problem of stability of such canardguided dualspinstabilized projectiles and gave the stability criterion considering the influence of canards. He pointed out that the stability for this kind of projectiles is reduced due to the following two factors: the addition of canards and the reduced angular momentum of the forward part. Theodoulis and Wernert [14] present the quasilinear parametervarying model of the dualspinstabilized projectile and analysed the stability of the projectile on different characteristic points of the whole flight envelope. Zhu [15] took account of the control effects of canards and then established a revised stability criterion according to the Hurwitz stability criterion which can be reduced to that of conventional spinstabilized projectiles.
However, few of the existing literatures have considered the coning motion of dualspin projectile with the control loop. Thus, this paper focuses on the stability of coning motion for a dualspin projectile with angular rate loops. The mathematical model of the dualspin projectile is established. The sufficient and necessary condition of coning motion stability for dualspin projectile with angular rate loops is analytically derived and further verified by numerical simulations. The stability boundary of control gains is obtained, and moreover, the influence of the total delay angle for the control system and spinning rate of the projectile on stability is analysed.
2. System Configuration
The socalled dualspin projectile in this paper is improved from a conventional 155 mm spinstabilized projectile by replacing the original fuse for a novel course correction fuse (CCF), which is illustrated in Figure 1. The projectile can be divided into two parts: the forward part contains the necessary guidance and control hardware and software and the aft part contains the payload. These two parts are connected by rotating bearings and are completely roll decoupled from each other.
As can be seen from Figure 1, the forward part is equipped with two pairs of rotating fins which are used to generate normal force to guide the projectile. After launching, the aft part keeps a high spinning rate of several hundreds of radians per second in order to maintain the fly gyroscopic stability, while the spinning rate of the forward part can be reduced to keep a constant low value in a few seconds under the control of an axial motor. Therefore, this dualspin structure design can not only solve the technical problems faced by the actuator when the projectile rotates at a high speed but also facilitate the design and normal work of the missileloading measuring device and can also maintain the gyroscopic stability of the projectile.
3. Mathematical Model
According to the previous studies [2], the 7DOF equations of motion of the dualspin projectile are generally established in the body fixed plane (BFP) reference frame, as shown in Figure 2.
The translational and attitude kinematics in the BFP reference frame can be derived as
The translational and attitude airframe dynamics can be derived in the BFP reference frame as where , , and are the linear velocities, , , and are the linear positions, the roll angle of the forward part, the roll angle of the aft part, the pitch angle, the yaw angle, the roll angular rate of the forward part, the roll angular rate of the aft part, and and the pitch and yaw angular rate. , , and are the external forces containing aerodynamic and gravity components and can be derived as
, , and are external moments containing aerodynamic and friction components:
The term is the drag coefficient, the lift force coefficient, the rollinduced moment coefficient, the rolldamping moment coefficient of the forward part, the rolldamping moment coefficient of the aft part, the static moment coefficient, the lateral damping moment coefficient, the Magnus moment coefficient, and the canard control moment coefficient. and are the equivalent variables of the canard in the pitch and yaw channels.
The definitions of and are illustrated in Figure 1 and can be expressed as
Other variables are the mass of the projectile , the dynamic pressure (where is the air density), the airframe velocity , the reference area and length , the moments of inertia , , and , and the static and viscous friction coefficients and .
It can be seen obviously from Equations (2) and (4) that the longitudinal and lateral dynamics of the dualspin projectile contain two main couplings. On the one hand, the inertial coupling term greatly depends on the high spin rate of the aft part. On the other hand, the aerodynamic coupling caused by Magnus moments depends on the angle of attack , angle of sideslip , and the high spinning rate . Furthermore, the aerodynamic coefficients appearing in Equations (5) and (6) are heavily nonlinear.
4. Angular Motion Equations of the DualSpin Projectile
Even though the mathematical model described in Equations (3) and (4) is more accurate and close to the real case, due to the highly nonlinear equations of motion, it is difficult to get the analytical solution and the obvious relationship between the flight characteristics of the missile and control parameters. To facilitate theoretical analysis, the general method is to apply the linearization theory of projectile. This theory has been regarded as an effective tool to analyse the flight stability of projectiles and applied [8–13]. Therefore, in order to linearize these two equations, the following assumptions are introduced firstly: (1)Relative to the longitudinal velocity and the angular velocity , variables , , , and can be regarded as small quantities, so we have (2)During a short flight, the varieties of the projectile parameters are neglected, and the velocity and the spinning rates are considered to be constant(3)The Magnus force and canard control force are negligible quantities compared to aerodynamic forces and the gravity effect is negligible(4)Angle of attack and angle of sideslip are small, so we have and .
Under these assumptions, the lateral equations for the dualspin projectile can be formulated as
According to assumptions (1) and (4), we can get
Thus, we can obtain the lateral equations of the dualspin projectile in a simplified form as where
By defining the complex variables , , and , we can further obtain the angular motion equation of the dualspin projectile in the complex form as
5. Stability of the Moving Mass Spinning Missile with the Angular Rate Loop
5.1. DualSpin Projectile Control System with Angular Rate Loops
The control system with angular rate loops is shown in Figure 3, in which and are control commands, and are feedback signals which can be measured by IMU mounted on the forward part, is the control gain, and and are equivalent canard deflections.
(a)
(b)
When the control commands for the pitch and yaw channels are and , the ideal equivalent output of the actuator system should also be and . However, due to the signal delay, when the actuator system receives the control command, the angular position of the forward part has shifted by . Furthermore, due to the dynamic response delay of the actuator, there will be a certain phase lag in the response output of the actuator system, which causes the actual equivalent output to generate an angular position deviation of . Therefore, the total delay angle of the control system are .
The actuator system is assumed to be a secondorder system, only considering the steadystate output coupling, we can get where is the time constant of the system, is the damping ratio, is the gain of the servo system, and is the command transmission delay time. Obviously, the control action in the pitch channel will cause the movement in the yaw channel, and the control action in the yaw channel will also affect the movement in the pitch channel.
It can be seen from Figure 3 that the input commands to the actuators can be described as
According to the definition of coordinate system and angle, positive equivalent canard deflections will generate positive angle of attack and negative angle of sideslip. And the negative angle of attack will generate negative pitching acceleration, and the positive angle of sideslip will generate positive yaw acceleration. Therefore, the displacement instruction of the canard is obtained as
Meanwhile, based on the assumption that the missile is in horizontal flight, there exists an approximation relationship: , . Thus, Equation (17) can be expressed as
Finally, substituting Equation (18) into Equation (13), we can get the actual equivalent canard outputs with all the delays as
Converting Equation (19) into the complex form, one has
5.2. Stability Conditions for the DualSpin Projectile Control System
Substituting Equation (20) into the angular motion Equation (12), we can get
Rewrite the above equation into the simplified form as where
Apparently, the characteristic equation of Equation (22) has the form
Then, we can obtain the characteristic roots as
Assuming that , one has
The real part of is where ; it is obvious that the control system is stable only if the real parts of the characteristic roots are negative, that is
The stability condition (28) can be formulated as
By substituting the expressions of these variables, the necessary and sufficient conditions for the stability of the dualspin projectile control system can be obtained as
5.3. Stability Analysis
Firstly, by analysing the first inequality in the stability condition (30), which is the pitching damping moment coefficient , it can be seen obviously that when the delay angle , the first inequality is satisfied. When , the control parameter has a upper limit as follows:
Then, analysing the second inequality in the stability condition (30), since the dualspin projectile is static unstable, we have . The first and the second terms of the second inequality are negative; thus, whether the second inequality is true depends on the third term. The magnitude and sign of the third term are mainly affected by the spinning rate of the aft part , the delay angle , and the control gain .
The second inequality can be rewritten as a polynomial in terms of : where
When , one gets , , and . The two zero points of are obtained as
Only when one gets . The sufficient and necessary condition for the coning motion stability can be derived as
When , one gets , , and ; the two zero points of can be obtained as
Only when or one gets . The sufficient and necessary condition for the coning motion stability can be derived as
6. Numerical Simulations
To demonstrate the proposed stability condition in the last section, numerical simulations are run for a 155 mm dualspin projectile. All simulation results are based on a fixedstep (0.001 s), fourthorder Runge–Kutta numerical program by using the MATLAB simulation software. The parameters of the dualspin projectile are listed in Table 1.

The delay angle of the control system is set as ; then, the calculated upper bound of the control loop gain can be obtained as 0.31 according to Equation (37). When the control loop gain , which satisfies the stability condition, the simulation results are obtained as shown in Figure 4. It can be seen obviously that the coning motion of the projectile converges to zero quickly. When , which is the critical value, the simulation results are illustrated as shown in Figure 3. It is observed that the coning motion of the missile neither converges nor diverges but presents a critical stable state. When , the simulation results are obtained as shown in Figure 5. It can be seen that the coning motion is divergent.
(a)
(b)
(a)
(b)
Furthermore, the influence of the total delay angle of the control system and the spinning rate of the aft part to the stability criterion is demonstrated. The stable regions corresponding to different total delay angles for the dualspin projectile control system are illustrated in Figure 6. It can be seen that the stable region of the control system shrinks as the delay angle increases. Therefore, in the design of the control system for the dualspin projectile, the delay angle should be reduced as far as possible to ensure that the control system is stable.
The stable regions corresponding to different spinning rates of the aft part are illustrated in Figure 7. It is noticed that the stable region of the control system shrinks as the spinning rate decreases. Furthermore, when the spinning rate of the aft part drops to a certain value, any control parameter cannot guarantee the stability of the dualspin projectile system. In the actual flight, the spinning rate of the aft part gradually decreases with flight time, which will reduce the stability region of the projectile system. For this reason, the flight stability of the dualspin projectile is more easily destroyed, which needs to be paid enough attention in the design process.
7. Conclusions
In this paper, the mathematical equation of a dualspin projectile is established. The sufficient and necessary condition of the coning motion stability for dualspin projectiles with angular rate loops is analytically derived and further verified by numerical simulations. Simulation results show that there exists a stability boundary value for the control gain. If the control gain exceeds it, the coning motion of the dualspin projectile will diverge, and the system will become unstable. It is also noticed that as the delay angle increases, the stability region of the system decreases, while, as the spinning rate of the aft part decreases, the stability region of the system decreases greatly. This paper is mainly based on the linearization theory of projectiles, so the stability condition obtained in this paper is applicable to the linearized missile model. In the future, we will focus on the stability analysis of nonlinear model of the dualspin projectile.
Data Availability
All data, models, and code generated or used during the study appear in the submitted article.
Conflicts of Interest
The authors declare that there is no conflict of interest regarding the publication of this paper.
References
 F. J. Reg and J. Smith, “Aeroballistics of a Terminally Corrected Spinning Projectile (TCSP),” Journal of Spacecraft and Rockets, vol. 12, no. 12, pp. 733–738, 1975. View at: Publisher Site  Google Scholar
 M. Costello and A. Peterson, “Linear theory of a dualspin projectile in atmospheric flight,” Journal of Guidance, Control, and Dynamics, vol. 23, no. 5, pp. 789–797, 2000. View at: Publisher Site  Google Scholar
 J. Guan and W. Yi, “Modeling of dualspinning projectile with canard and trajectory filtering,” International Journal of Aerospace Engineering, vol. 2018, Article ID 1795158, 7 pages, 2018. View at: Publisher Site  Google Scholar
 B. Burchett, A. Peterson, and M. Costello, “Prediction of swerving motion of a dualspin projectile with lateral pulse jets in atmospheric flight,” Mathematical and Computer Modelling, vol. 35, no. 78, pp. 821–834, 2002. View at: Publisher Site  Google Scholar
 S. Chang, Z. Wang, and T. Liu, “Analysis of spinrate property for dualspinstabilized projectiles with canards,” Journal of Spacecraft and Rockets, vol. 51, no. 3, pp. 958–966, 2014. View at: Publisher Site  Google Scholar
 S. Chang, “Dynamic response to canard control and gravity for a dualspin projectile,” Journal of Spacecraft and Rockets, vol. 53, no. 3, pp. 558–566, 2016. View at: Publisher Site  Google Scholar
 X.d. Liu, D.g. Li, and Q. Shen, “Swerving orientation of spinstabilized projectile for fixedcant canard control input,” Mathematical Problems in Engineering, vol. 2015, Article ID 173571, 11 pages, 2015. View at: Publisher Site  Google Scholar
 P. Wernert, F. Leopold, D. Bidino, and J. Juncker, “Wind tunnel tests and openloop trajectory simulations for a 155 mm canards guided spin stabilized projectile,” in AIAA Atmospheric Flight Mechanics Conference and Exhibit, Honolulu, Hawaii, August 2008. View at: Google Scholar
 C. H. Murphy, Free Flight Motion of Symmetric Missiles, Ballistic Research Laboratories, Aberdeen Proving Ground, 1963.
 C. H. Murphy, “Angular motion of spinning almost symmetric missiles,” Journal of Guidance, Control and Dynamics, vol. 2, no. 6, pp. 504–510, 1979. View at: Publisher Site  Google Scholar
 K. L. Nielsen and J. L. Synge, “On the motion of a spinning shell,” Quarterly of Applied Mathematics, vol. 4, no. 3, pp. 201–226, 1946. View at: Publisher Site  Google Scholar
 A. E. Hodapp, “Effect of Mass Asymmetry on Ballistic Match of Projectiles,” Journal of Spacecraft and Rockets, vol. 13, no. 12, pp. 757–760, 1976. View at: Publisher Site  Google Scholar
 P. Wernert, “Stability analysis for canard guided dualspin stabilized projectiles,” in AIAA Atmosphere Flight Mechanics Conference, Chicago, IIIinois, August 2009. View at: Publisher Site  Google Scholar
 S. Theodoulis and P. Wernert, “Flight control for a class of 155mm spinstabilized projectiles with course correction fuse (CCF),” in AIAA Guidance, Navigation, and Control Conference, Portland, Oregon, August 2011. View at: Publisher Site  Google Scholar
 D. Zhu, S. Tang, J. Guo, and R. Chen, “Flight stability of a dualspin projectile with canards,” Proceedings of the Institution of Mechanical Engineers, Part G: Journal of Aerospace Engineering, vol. 229, no. 4, pp. 703–716, 2014. View at: Publisher Site  Google Scholar
Copyright
Copyright © 2020 Qiushi Zheng and Zhiming Zhou. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.