This paper presents an optimum control scheme of firing time and firing phase angle by taking impact point deviation as optimum objective function which takes account of the difference of longitudinal and horizontal correction efficiency, firing delay, roll rate, flight stability, and so forth. Simulations indicate that this control scheme can assure lateral impulse thrusters are activated at time and phase angle when the correction efficiency is higher. Further simulations show that the impact point dispersion is mainly influenced by the total impulse deployed, and the impulse, number, and firing interval need to be optimized to reduce the impact point dispersion of rockets. Live firing experiments with two trajectory correction rockets indicate that the firing control scheme works effectively.
1. Introduction
Dispersion characteristics of the trajectory correction projectiles can be dramatically improved by outfitting with a suitable trajectory correction flight control system. The commonly used executive organs are moveable canards [1–3], fixed canards [4, 5], synthetic jets [6], drag brakes [7, 8], and impulse thrusters [9–21].
The work reported here describes the implementation of a flight control system with impulse thrusters on a representative 122 mm artillery rocket. The flight control consists of a finite number of impulse thrusters mounted forward on the rocket body, computes position and velocity errors through comparing the position and velocity measured by GPS with prespecified (reference) trajectory, fires thrusters to change velocity direction, and assists the rocket to follow a prespecified trajectory.
Research and development on the use of impulse thrusters in order to improve the precision of projectiles has been going on for decades. Brandeis and Gill [9] have studied the effect of the impulse thruster on the forces and moments and investigated the impact of size and location of planar surfaces on the amplification of the impulse thruster. Burchett et al. [10] developed a simple set of formulae for the swerve motion of a dual-spin projectile under the action of lateral pulse jets and pointed out that the ability of lateral pulse jets to change the impact point steadily decreases with range to the target. Corriveau et al. [11] described a novel method for correcting the trajectory of spin- and fin-stabilized projectiles using pairs of impulse thrusters located away from the center of mass; they used the linear theory of ballistics to develop an analytical model representing the motion of a projectile subjected to impulse thrusters.
The thrusters’ application on rockets has been originally considered by Harkins and Brown [12]. They have proposed a method using a set of impulse thrusters to marginalize the off-axis angular rates of the rockets just after exiting the launcher and managed to reduce the impact point dispersion by the factor of 4. Jitpraphai and Costello [13] have proposed a simplified control system with impulse thrusters and demonstrated that impact point dispersion of a direct fire rocket could be drastically reduced. Burchett and Costello [14] have investigated the use of a small number of short-duration lateral pulses acting as a control mechanism to reduce the dispersion of a direct fire rocket and studied the effect of the number and amplitude of individual pulse jets. Gupta et al. [15] have explored the feasibility of reducing the impact point dispersion of a routinely used artillery rocket using lateral pulsejets coupled to a trajectory correction flight control system. Pavkovic et al. [16] have presented a new two-phase guidance scheme named trajectory tracking with pulse-frequency modulation which is presented, considering the fact that an artillery rocket flies through different atmospheric environments; the presented guidance scheme achieves excellent accuracy even in the case of a small number of pulse jets. Recently, Pavkovic et al. [17] have presented a simplified control scheme for artillery rockets named the active damping method which performs a correction of disturbances immediately after a rocket exits a launcher tube. It is shown that the application of such a control system achieves a significant dispersion reduction.
Each impulse thruster imparts a single, short-duration, large force to the rocket in the plane normal to the rocket axis of symmetry, the control scheme of impulse thrusters mainly involves two aspects: the firing time and the firing phase angle. Many researchers [13–15] compute firing phase angle by the phase angle of trajectory deviation; theoretical analysis shows that this method may not get the best trajectory correction performance under some conditions, and the induced effects of rocket flight time on trajectory performance are not considered. Yang et al. [18] converts the problem of determining control parameters of the pulsejets into the design of the experiment with multifactor and multilevel. Introducing the firing time as a controllable factor, the relationship model among firing time, number of pulsejets, and total value of trajectory correction was obtained by using regression analysis. This scheme can determine the firing time and number of pulsejets needed quickly and reduce the trajectory error of the rocket effectively, but the difference of longitudinal and horizontal correction efficiency is not considered. Cao et al. [19] present an optimization strategy for the firing phase angle, which makes the total number of thrusters minimized for the residual trajectory deviation after previous correction. However, the foremost goal of trajectory correction is to decrease the trajectory error; the consumption number of thrusters should be taken as the secondary factor.
In this paper, to reduce the impact point dispersion of the trajectory correction rocket using impulse thrusters, a 6-DOF trajectory model with lateral force is established, and then the control algorithm of firing time and firing phase angle is put forward, taking impact point deviation as optimum objective function.
2. 6-DOF Trajectory Model with Lateral Force
The numerical simulation is based on a rigid body six degrees of freedom model typically utilized in flight dynamic analysis of rockets. Figure 1 shows the rocket configuration with an impulse thruster ring mounted on the forward part of rockets. The thrusters are assumed to be located at l from the center of mass. The 6-DOF trajectory model with lateral force can be established by taking the lateral force and its moment into the ballistic motion equations.
Schematic of layout and lateral force of thrusters.
Layout of thrusters
Lateral force
The lateral force in quasibody reference frame can be described as(1)Tx4Ty4Tz4=0TcosθTsinθ.In (1), T is the lateral force of an impulse thruster in quasibody reference frame, and θ is the phase angle of lateral force.
Transforming the lateral force from the quasi body to the aeroballistic reference frame, we have(2)Tx2Ty2Tz2=cosαcosβ-sinαcosβsinβsinαcosα0-cosαsinβsinαsinβcosβ0TcosθTsinθ.In (2), α is the angle of attack of the rocket, and β is its sideslip angle.
The translational kinetic differential equations of the rocket in aero-ballistic reference frame are given in(3)dVdt=Px2+Fx2+Tx2,Vdθ′dt=Py2+Fy2+Ty2,-Vcosθ′dψvdt=Pz2+Fz2+Tz2.The applied loads appearing in (3) consist of main rocket thrust (P), lateral thruster force (T), and other forces (F) components. V, θ′, and ψv are the velocity, trajectory incline angle, and flight path azimuth angle of the rocket, respectively.
Lateral moment of an impulse thruster in quasibody reference frame is given by(4)MTx4MTy4MTz4=0-Tz4lTy4l.The rotational kinetic differential equations of the rocket in quasibody reference frame are given by(5)Jx4dωx4dtJy4dωy4dtJz4dωz4dt=Mx4+MTx4My4+MTy4Mz4+MTz4-0Jx4-Jz4ωx4ωz4Jy4-Jx4ωx4ωy4+0-Jz4ωz4dγdtJy4ωy4dγdt.The applied moments appearing in (5) contain contributions from lateral thruster forces, denoted by MT, and other forces, denoted by M. γ is the Euler roll angle of the rocket. wx4, wy4, wz4 are components of the angular rate vector. Jx4, Jy4, Jz4 are components of the transverse moment of inertia.
Other motion equations of the rocket not involving lateral forces or lateral moments can be obtained in [21].
3. The Establishing of Optimum Objective Function
Figure 2 shows trajectory correction progress of an impulse thruster. Assuming that M is the target and O is the predicted impact point computed by flight control system, the distance between the predicted impact point and the target is OM, and the longitudinal and horizontal deviations are ΔL and ΔH, respectively. Assuming that the longitudinal and horizontal correction distance of an impulse thruster is Px and Pz, respectively, if the impulse thruster is activated, the predicted impact point of the rocket will be changed to A; then the distance between the predicted impact point and the target will be reduced to be AM upon the impact of an impulse thruster.
Schematic of trajectory correction.
Because the aim of trajectory correction is to reduce the impact point deviation, it is foremost to make best use of the energy of impulse thrusters to obtain the minimum impact point deviation. Therefore, the remaining impact point deviation after an impulse thruster was activated, namely, AM, and is taken as the optimization objective function in this paper. It is assumed that the positive direction of ΔL and Px is the same, and the positive direction of ΔH and Pz is the same, and then AM can be described by(6)AM=ΔL+Px2+ΔH+Pz21/2.In (6), ΔL and ΔH are the longitudinal and horizontal predicted deviation, computed by the flight control system and remain constant during a short time, so the change of AM mainly relates to Px and Pz of the impulse thruster activated. Px and Pz depend on roll rate of the rocket, impulse, firing time, and firing angle phase of the thruster; individual thruster impulse and roll rate of the rocket have been determined in the design phase of trajectory correction rockets; therefore, the firing time and firing phase angle are taken as optimization variables of the objective function in the flight phase.
4. Optimization of Firing Phase Angle
The optimization variables firing time and firing phase angle are relatively independent and assumed that the firing time of an impulse thruster has been identified; then the firing phase angle is the only optimization variable of an impulse thruster. To analyze the impact of firing phase angle on the trajectory correction performance, the lateral force is seen as a constant because the thrusters are active over a very short duration of time, and θ denotes the phase angle of lateral force in quasibody reference frame. To get the best trajectory correction performance, the best phase angle of lateral force leads to the minimum value of AM, is denoted by θm, and should be computed firstly. The conventional method to compute θ is by making its value equal to phase angle of trajectory deviation, denoted by α, referencing to Figure 3(a). However, whether or not this method can get the best trajectory correction performance depends on the longitudinal and horizontal correction efficiency (namely, as converting coefficient) of the impulse thrusters.
Trajectory correction performance with equal and unequal converting coefficients.
With equal converting coefficients
With unequal converting coefficients
Given that the correction distance of an impulse thruster is proportional to the thruster impulse, Px and Pz can be described by(7)Px=kx·cosθ,Pz=kz·sinθ.In (7), I is the impulse of a thruster. kx and kz are longitudinal and horizontal converting coefficients of an impulse thruster for the rocket. The converting coefficients reflect the correction efficiency of the impulse thruster. The impact of correction efficiency on the objective function will be analyzed under two conditions.
(1) When kx=kz, longitudinal and horizontal correction efficiency of an impulse thruster is the same and assumed that O is the predicted impact point of the rocket; then possible impact points after correction of an impulse thruster form a circle around O. Assuming that the correction distances are OA1 and OA2, corresponding to θ=α and θ≠α, then OA1=OA2, and OA1 has the same direction with OM, referring to Figure 3(a). The following equation can be obtained according to triangle trilateral theorem:(8)OA1+A1M<OA2+A2M.Because OA1=OA2, A1M<A2M. Therefore, when kx=kz, objective function AM can get the minimum value with θ=α.
(2) When kx≠kz, longitudinal and horizontal correction efficiency of an impulse thruster is different, whether or not OA1=OA2 is not certain and, referring to Figure 3(b), A1M may be bigger than A2M, so objective function cannot always get the minimum value with θ=α. Under this condition, a new method is required to compute θm.
In fact, kx≠kz in most of the flight time, as shown in Figure 6, a new method to compute best phase angle of lateral force is proposed in this paper. Because an analytic solution is not acquired by (6), θm is solved by the binary iteration method given in Figure 4, where [θ1,θ2] is the solution limits and N is the number of iterations.
Solution flow of θm.
To reduce the number of iterations and improve the solution accuracy, solution limits can be acquired by (9).
The solution accuracy can reach 0.1 deg with N=10, so θm can be solved quickly by this method.
The firing phase angle, denoted by Φc, is the phase angle of the impulse thruster in quasibody reference frame when it gets the firing signal from the flight control system. Firing time delay and impulse duration time should be taken into account to make sure the average angle of lateral force is equal to θm:(9)θ1,θ2=0∘,90∘ΔH<0,ΔL<090∘,180∘ΔH<0,ΔL>0180∘,270∘ΔH>0,ΔL>0270∘,360∘ΔH>0,ΔL<0.Schematic of firing phase angle is given in Figure 5. The relation between firing phase angle Φc and θm is given by(10)Φc=θm+π-ωτd+τ.In (10), τd is the firing delay, and τ is a half of the impulse duration. ω is the roll angular rate of rocket.
Schematic of firing phase angle.
Converting coefficients versus flight time.
When the phase angle Φ of the impulse thruster that will be activated is equal to the firing phase angle, the impulse thruster can be activated; the condition to activate an impulse thruster is (11)Φ-Φc≤ε.In (11), ε is the desired activation threshold.
5. Optimization of Firing Time
After the optimization of firing phase angle, the other optimization variable of objective function is the firing time of thrusters. The general firing time control algorithm has two strategies.
(1) Time elapsed from the previous impulse thruster firing must be longer than a specified time duration Δtfire:(12)t-t∗>Δtfire.In (12), t∗ is the firing time of the previous impulse thruster.
(2) If predicted impact point deviation is greater than a specified distance, activate the thruster as soon as possible.
Δtfire is an important design parameter of thrusters. If it is set too low, the rocket does not have sufficient time to respond and too many thrusters will be fired, tending to overcompensate for trajectory deviation. Simultaneously, a lower Δtfire may have an impact on the flight stability of the rocket. On the other hand, if it is set too high, only a small number of thrusters can possibly be fired and the capabilities of the residual impulse thrusters will be wasted.
The impact of flight time on trajectory correction performance of thrusters is not considered by strategy (2), which may lead to the thrusters being activated at the time when the correction efficiency is lower. Gao and Zhang [20] gave the conclusion that the longitudinal correction efficiency of thrusters is lower in ballistic ascending segment of the rocket, and longitudinal and horizontal correction efficiency reduces with the decrease of time to go in ballistic descending segment. Simulation results of converting coefficients of a rocket are shown in Figure 6, the launch angle of the rocket is 45 deg, and the impact angle is −64 deg.
As shown in Figure 6, kx and kz have different change trends, and kx increases gradually in ballistic ascending segment and decreases gradually in ballistic descending segment, while the kz decreases gradually with the decrease of time to go. The behavior of kx and kz has a physical reason. According to the kinetic energy theorem, the rocket will get a velocity increment ΔV when an impulse thruster was activated, the direction of ΔV is nearly perpendicular to the velocity direction of rocket. When an impulse thruster was activated for longitudinal correction, ΔV can be decomposed into ΔVx and ΔVz, as shown in Figure 7; supposing that the lateral force of the impulse thruster is upward, then ΔVy=ΔVcos(θ′), and ΔVx=-ΔVsin(θ′), where θ′ denotes the trajectory incline angle. In ballistic ascending segment, θ′>0, so ΔVx<0, which will shorten the rocket range, while ΔVy>0, which will extend the rocket range, the impact of ΔVx and ΔVy is contradictory, and, therefore, the impact of ΔV is relatively small at the beginning of ballistic descending segment. As the trajectory inclination angle ranges from 45 deg to 0 deg, the absolute value of ΔVy increases gradually, while the absolute value of ΔVx decreases gradually; as a result, the longitudinal correction distance of a thruster increases gradually in ballistic ascending segment. In ballistic descending segment, trajectory inclination angle ranges from 0 deg to −64 deg, so ΔVx>0, which will extend the rocket range, and ΔVy>0, which will also extend the rocket range; the longitudinal correction distance is mainly influenced by the time to go; therefore, longitudinal correction distance of a thruster decreases gradually with the decrease of time to go in ballistic descending segment. In the horizontal plane, ΔVz≈ΔV, the horizontal correction distance is mainly influenced by the time to go, and hence the horizontal correction distance decreases gradually with the decrease of time to go. Therefore, kx and kz have the change trends shown in Figure 6.
Schematic of velocity increment of a thruster.
It can be known from (6) and Figure 6 that the value of the objective function (activated at the optimized firing phase angle) tends to increase in ballistic descending segment. Supposing that tk is the time of ballistic vertex, the value of the objective function at tk can be taken as a reference value, denoted by AM(tk). In ballistic ascending segment, thrusters can be activated if the value of the objective function is less than the reference value AM(tk), otherwise they cannot. This method can guarantee that thrusters are activated at the time when correction efficiency is higher, especially in the ballistic ascending segment.
To sum up, the firing time control algorithm can be improved as the following equation:(13)t-t∗>Δtfire,AMt<AMtkt<tk,AMt<OMtt≥tk.In (13), tk is the time of ballistic vertex.
6. Results and Discussion
To investigate the correction performance of impulse thrusters and verify the effectiveness of the firing control scheme, some simulations of a rocket have been done by numerical integration of the equations described above using a fourth-order Runge-Kutta algorithm. The rocket configuration used in the simulation study is a representative 122 mm artillery rocket, 2.99 m long, fin-stabilized, with four pop-out fins on its rear part. The main rocket motor burns for 2.55 s and imparts an impulse of 54247 N-s to the rocket. During the main rocket motor burns, the forward velocity of the rocket is increased from 46.9 m/s to 935.7 m/s. The rocket weight, mass center location from the nose tip, roll inertia, and pitch inertia before and after burn is 66.1/43.0 kg, 1.43/1.21 m, 0.16/0.12 kg-m2, and 48.42/36.36 kg-m2, respectively. The rocket is launched at sea level toward a target on the ground with altitude and cross range equal to zero at a range of 28000 m. The thruster ring is assumed to be located at 0.869 m from the nose tip of the rocket and contains 50 individual thrusters where each individual thruster imparts an impulse of 15 N-s on the rocket body over a time duration of 0.02 s. The minimum firing interval of thrusters Δtfire is set to 0.2 s. The desired activation threshold ε is set to 3 deg.
The time-varying data of uncontrolled and controlled trajectories with optimum firing control scheme against a nominal command trajectory for the example rocket are compared in Figure 8.
Correction performance of the optimum firing control scheme.
Total number of thrusters tired versus time
Total attack angle versus time
Horizontal impact point deviation versus time
Longitudinal impact point deviation versus time
Cross range versus range
Attitude versus range
Figure 8(a) plots the total number of thrusters fired versus time. Five thrusters are fired in 10~15 s, and 43 thrusters are fired in 45~70 s. No thrusters are fired in 15~45 s which indicates that the correction efficiency is relatively low during this period. No thrusters are fired after 70 s because the predicted impact point deviation is too small to implement any trajectory correction. The activation of thrusters will arouse the increase of attack angle of rocket, as shown in Figure 8(b).
It can be known from Figures 8(c) and 8(d) that the predicted longitudinal and horizontal impact point deviations reduce gradually. The horizontal impact point deviation reduces quickly, while longitudinal impact point deviation has no significant change in 10~20 s, which indicates that the thrusters fired in this period are mainly used for horizontal correction, corresponding to the fact that the horizontal correction efficiency of thrusters in this period is higher (kz is greater, referencing to Figure 6). Both horizontal and longitudinal impact point deviations reduce quickly in 45~70 s, which indicates that thrusters fired in this period are used for both horizontal and longitudinal correction, corresponding to the fact that the horizontal and longitudinal correction efficiency in this period is similar (kx≈kz).
Figures 8(e) and 8(f) show that the trajectory correction rocket flights to target point gradually under the effect of thrusters. The final impact point deviation of the rocket was reduced from horizontal 341.3 m and longitudinal 483.3 m in the uncontrolled case to horizontal 10.9 m and longitudinal 15.9 m in the controlled case, which demonstrates that the firing control scheme works effectively.
Figure 9 shows the impact point distribution using the Monte Carlo method. The cases of the uncontrolled rockets as well as of the rockets with general firing control scheme and optimum firing control scheme are shown. The analysis was performed for a statistical sample of 1024 simulations. All initial conditions were modeled as Gaussian random variables, with mean values and standard deviations given in Table 1. Wind direction is a uniform random variable between 0 and 2π. The nominal launching elevation angle is 42.31 deg, but the mean value of launching elevation angle is set as 42.51 deg, a launching elevation angle error about 0.2 deg and a launching azimuth angle error about 0.1 deg were introduced to simulate the accuracy errors of rockets.
Initial conditions and disturbances.
Parameter
Unit
Mean value
Standard deviation
Launching elevation angle
deg
42.51
0.2
Launching azimuth angle
deg
0.1
0.2
Impulse of the main rocket motor
%
100
0.1
Wind
m/s
0
2
Impact point distribution.
Impact point distribution of uncontrolled rockets
Impact point distribution of rockets with general firing control scheme
Impact point distribution of rockets with the optimum firing control scheme
The result of Monte Carlo simulations is shown in Figure 9, for uncontrolled rockets; the max value of range miss is 589.2 m, while the min value of range miss is −421.9, the mean value of range miss is 95.5, the max value of cross range miss is 1349.8, while the min value of cross range miss is −1063.6, the mean value of cross range miss is 43.8, there are more points in the first quadrant (right upper side) of the plane coordinate system, as shown in Figure 9(a), and the CEP of uncontrolled rockets is 359 m. Impact point distribution of the trajectory correction rockets reduces greatly, as shown in Figures 9(b) and 9(c). The rockets with the general firing control scheme have a CEP of 38 m, while the CEP of rockets with the optimum firing control scheme is 20 m. The average thruster consumption and its standard deviation of the rockets with the general firing control scheme are 31.7 and 15.5, while the average thruster consumption and its standard deviation of the rockets with the optimum firing control scheme are 25.1 and 13.3. The decrease of impact point dispersion of rockets and the reduction of thruster consumption testify the effectiveness of firing control optimization.
Figure 10 shows the impact of thruster parameters on the CEP and flight stability of rockets. Figure 10(a) shows the relation between CEP of rockets, number of thrusters mounted on rocket, and individual thruster impulse. It can be known that the CEP of rockets is highly correlated with the number of thrusters and individual thruster impulse, and CEP reduces steadily as the number of thrusters or individual thruster impulse is increased.
The impact of parameters of thrusters on CEP and flight stability of rockets.
CEP versus number of thrusters and individual thruster impulse
CEP versus total impulse of thrusters and individual thruster impulse
Maximal total attack angle versus individual thruster impulse and minimum firing interval
Figure 10(b) shows the relationship between CEP of rockets, individual thruster impulse, and the total impulse of thrusters. For a value of the total impulse, as the individual thruster impulse is increased, the number of thrusters decreases proportionally for the total impulse to remain as a constant. As shown in Figure 10(b), the evolution of the three curves is similar; when the thruster impulse is small, the number of thrusters is large, and some thrusters will be activated at time when the correction efficiency is lower, which is the unfavorable aspect for trajectory correction; however, large thruster number means more correction opportunity, and the correction precision of a small thruster impulse is higher, which is the favorable aspect for trajectory correction; these two aspects are contradictory, the combined effect of these two aspects leads to the similar evolution of the three curves. The result indicates that CEP is mainly influenced by the total impulse when the individual thruster impulse is reasonable, CEP reduces gradually as the total impulse is increased, and the value of individual thruster impulse has a small impact on CEP.
Figure 10(c) shows the relationship between maximum total attack angle, individual thruster impulse, and the minimum firing interval (Δtfire). The results indicate that individual thruster impulse and minimum firing interval have a direct impact on the maximum total attack angle; as individual thruster impulse increases or minimum firing interval reduces, the maximum total attack angle increases, which may affect the flight stability of trajectory correction rockets.
Based on the above discussion, some suggestions are put forward to determine the thruster configuration parameters of trajectory correction rockets. Firstly, the total impulse should be determined according to CEP needed, because the impact point dispersion is mainly influenced by the total impulse of lateral thrusters deployed. Secondly, individual thruster impulse should be determined according to the total impulse, limits of layout space on the rocket, flight stability of rockets, cost, and so forth. If the individual thruster impulse is too small, there may be too many thrusters needed to be mounted on the rocket. On the other hand, if the individual thruster impulse is too large, the flight stability may deteriorate. Thirdly, the minimum firing interval should be determined according to the individual thruster impulse, flight stability of rockets, correction efficiency of thrusters, and so forth. If the individual thruster impulse is large, a relatively large minimum firing interval should be set to guarantee the flight stability of rockets. If the individual thruster impulse is small, a relatively small minimum firing interval time should be set to insure thrusters can be activated at the segment of trajectory while the correction efficiency is higher.
Live experiments have been done to verify the effectiveness of the firing control scheme. The characters of the trajectory correction rockets (TCR) have been launched same as the rocket described in simulations above. Two TCRs have been launched, toward a target on the ground with altitude equal to 200 m and cross range equal to zero at a range of 28000 m; the launching elevation angle was set to 43.6 deg, according to the weather of the time. Figure 11 shows the trajectory of two TCRs contrasted with the nominal trajectory. As shown in Figure 11(a), at the ballistic vertex, the attitude of two TCRs is 300~400 m higher than the nominal trajectory; then the trajectory of two TCRs gradually close to the nominal trajectory under the effect of flight control system. As shown in Figure 11(b), the cross range of the two TCRs gradually closes to the nominal trajectory under the effect of flight control system. The final impact point deviation of TCR1 is (4.5,21.2) and the final impact point deviation of TCR2 is (-28.3,35.4), which demonstrates the firing control scheme works effectively.
The trajectory of two TCRs.
Attitude versus range
Cross range versus range
7. Conclusion
This paper establishes the 6-DOF trajectory model of a rocket with lateral force and presents an optimum control scheme of firing time and firing phase angle by taking impact point deviation as optimum objective function which takes account of the difference of longitudinal and horizontal correction efficiency, firing delay, roll rate, flight stability, and so forth. The result of Monte Carlo simulations shows that the uncontrolled rockets have a CEP of 359 m; the CEP was improved to 38 m with the general firing control scheme; as a contrast, the CEP was improved to 20 m with the optimum firing control scheme. The average thruster consumption and its standard deviation of the rockets with the general firing control scheme are 31.7 and 15.5, while the average thruster consumption and its standard deviation of the rockets with the optimum firing control scheme are 25.1 and 13.3. The decrease of CEP and the reduction of thruster consumption testify the effectiveness of firing control optimization. The variations of rocket impact point dispersion are analyzed with different impulse and number of impulse thrusters. It is shown that the impact point dispersion is mainly influenced by the total impulse of impulse thrusters deployed and steadily decreases as the total impulse is increased. The impulse, number, and firing interval need to be optimized to insure the flight stability of rockets and impulse thrusters activated at time when the correction efficiency is higher. Two trajectory correction rockets have been launched to verify the effectiveness of the firing control scheme; the final impact point deviations of two trajectory correction rockets are, respectively, (4.5,21.2) and (-28.3,35.4), which demonstrates the firing control scheme works effectively.
Conflict of Interests
The authors declare that there is no conflict of interests regarding the publication of this paper.
CostelloM.Extended range of a gun launched smart projectile using controllable canards200183-420321310.1155/2001/6157482-s2.0-0035191373RogersJ.CostelloM.Design of a roll-stabilized mortar projectile with reciprocating canards20103341026103410.2514/1.478202-s2.0-77956374134CooperG.FresconiF.CostelloM.Flight stability of an asymmetric projectile with activating canards201249113013510.2514/1.A320222-s2.0-84857601124CostelloM.PetersonA.Linear theory of a dual-spin projectile in atmospheric flight200023578979710.2514/2.46392-s2.0-0034273221ElsaadanyA.Wen-JunY.Accuracy improvement capability of advanced projectile based on course correction fuze concept201420141027345010.1155/2014/273450AmitayM.SmithD. R.KibensV.ParekhD. E.GlezerA.Aerodynamic flow control over an unconventional airfoil using synthetic jet actuators200139336137010.2514/2.13232-s2.0-0035280605HollisM. S. L.BrandonF. J.Range correction module for a spin stabilized projectileU.S. Patent Documents US-5816531, 1998HollisM. S. L.BrandonF. J.Design and analysis of a fuze-configurable range correction device for an artillery projectile1999ARL-TR-2074Aberdeen, Md, USAU.S. Army Research Laboratory, Aberdeen Proving GroundBrandeisJ.GillJ.Experimental investigation of super- and hypersonic jet interaction on missile configurations199835329630210.2514/2.33542-s2.0-0032069467BurchettB.PetersonA.CostelloM.Prediction of swerving motion of a dual-spin projectile with lateral pulse jets in atmospheric flight2002357-882183410.1016/S0895-7177(02)00053-5MR1901291ZBL1054.700212-s2.0-0037177271CorriveauD.WeyP.BernerC.Thrusters pairing guidelines for trajectory corrections of projectiles2011344112011282-s2.0-7995997298510.2514/1.51811HarkinsT.BrownT.Using active damping as a precision-enhancing technology for 2.75-inch rockets1999ARL-TR-1772Aberdeen, Md, USAU.S. Army Research LaboratoryJitpraphaiT.CostelloM.Dispersion reduction of a direct fire rocket using lateral pulse jets200138692993610.2514/2.37652-s2.0-0035518661BurchettB.CostelloM.Model predictive lateral pulse jet control of an atmospheric rocket200225586086710.2514/2.49792-s2.0-0036749386GuptaS. K.SaxenaS.SinghalA.GhoshA. K.Trajectory correction flight control system using pulsejet on an artillery rocket2008581153310.14429/dsj.58.16212-s2.0-39749094478PavkovicB.PavicM.CukD.Frequency-modulated pulse-jet control of an artillery rocket201249228629410.2514/1.A321332-s2.0-84861503593PavkovicB.PavicM.CukD.Enhancing the precision of artillery rockets using pulsejet control systems with active damping20126221019YangH.-W.DouL.-H.GanM.-G.Control method of impulse actuator of antiaircraft trajectory correction projectile20113110118411882-s2.0-82055191055CaoY.-J.YangS.-X.LiX.-J.Research on firing phase angle optimization of terminal trajectory correction projectile based on pulse jet control20082988979012-s2.0-52249102013GaoF.ZhangH.Study of 2-D trajectory correction based on geomagnetic detection with impulse force for projectiles20112311231282-s2.0-79951545766QianX.LinR.ZhaoY.2008Beijing, ChinaBeijing Institute of Technology Press