1 en 12 Chapter Author
1 en 12 Chapter Author
1 en 12 Chapter Author
Author Proof
Abstract. One of the current tendencies is to equip the unguided munition AQ1
designed for ballistic shooting with low cost, solid propellant lateral thrusters
based, actuators to achieve the controlled flight functionality and reduce the
collateral damage. The basic technical challenge connected with this type of
pulsed control system is that each from the thrusters could be used only once
which results in projectile low control authority. The thruster parameters have
crucial impact on the achieved hit accuracy. The main goal of this article was to
investigate and understand the influence of control force magnitude on the
guidance process using six-degree-of-freedom numerical simulation. 122 mm
artillery rocket controlled with the aim of 30 solid propellant thrusters mounted
before center of mass was chosen as a test object. Single channel control was
considered. The impact point prediction algorithm based on point mass model
was developed and implemented into a Matlab software. Using Monte-Carlo
simulations the optimum lateral thruster force amplitude was obtained for shots
at low elevation angle. The numerical experiments showed that with the pro-
posed method the circular error probable of the projectile might be reduced 10
times when compared to unguided case.
1 Introduction
One of the modern requirements is to improve the delivery accuracy of the existing
unguided munitions. The object of interest of this work is concentrated on rocket
artillery. The main advantages of this kind of weapon are: long range, high firepower
and high mobility of the launcher [1]. The main disadvantage is high dispersion of the
impact points which makes its area weapon. There exists a tendency to equip the
unguided projectiles with a various kind of actuators like body-fixed canards [2–4], roll
decoupled fins [5], lateral thrusters [6, 7], foldable drag plates and movable internal
parts to achieve the controlled flight functionality. Using this approach, the collateral
damage and the number of projectiles required to eliminate the target could be reduced.
On the other hand, surface to surface rockets have low control authority which
limits the amount of possible control strategies. The guidance methods designed for
guided munition fall into three main categories [5]: reference trajectory tracking,
impact point prediction and trajectory shaping. The main requirements for control
algorithms intended for smart weapons are the guidance accuracy and low computa-
Author Proof
tional burden due to limited resources available onboard of the rocket. The trajectory
tracking method [8] relies on flight along prespecified ballistic curve. The path defi-
nition is preloaded onto memory of the onboard computer before the projectile launch.
The advantage of this method is its low cost. The predictive guidance [9, 10] uses a
model of the projectile implemented into a control system to obtain the impact point
location coordinates. In this method the control effort rely on driving the impact point
to the target. The trajectory shaping algorithm is based on the following non-ballistics
trajectories [11]. This method allows to for e.g. extend the range and control the impact
angle of the projectile but the main drawback of this approach is that generally requires
more control authority than predictive guidance.
The main goal of the article was to further extend the study on pulse thrusters to
minimize the dispersion of the rocket and obtain the most suitable thrust force
amplitude of the single lateral thruster when the predictive guidance is applied.
The parameters before/after main motor burnout was assumed to be: center of mass
location xcg0 = 1.446/xcgk = 1.645 m from the rocket base, mass m0 = 64.75/
mk = 38.5 kg, axial moments of inertia Ix0 = 0.137/Ixk = 0.092 kgm2 and lateral
moments of inertia Iy0 = Iz0 = 38.80/Iyk = Izk = 29.16 kgm2. In the basic unguided
configuration, the rocket parameters were designed for ballistic flight.
The functionality of controlled flight was obtained with equipping the rocket with a
guidance unit composed from power supply, inertial navigation system, onboard
computer and a set of 30 small solid propellant lateral thrusters mounted before the
center of mass and located 1.9 m from the rocket base. The nozzles of thrusters are
spaced equally around the fuselage. In this solution the thrusters generate not only force
perpendicular to the rocket longitudinal axis but also the moment around its center of
mass. The main difficulty connected with this type of actuator is the fact that each
thruster could be fired only once. The limited energy stored in the grain of the thrusters
and limited number of firings provides very low control authority of the rocket. It
means that the achievable lateral accelerations are much lower than in case of high
maneuver air-to-air missiles and the resulting the impact point could be translated only
Monte-Carlo Based Lateral Thruster Parameters Optimization for 122 mm Rocket 3
where U, V, W are rocket linear velocities in body-fixed coordinate frame, P, ::, R are
rocket roll, pitch and yaw rates, U, H, W are Euler angles (roll, pitch and yaw), xn , yn ,
zn are position coordinates of the rocket center of mass in the north-east-down frame.
The m is projectile mass. Ix is the longitudinal and Iy , Iz lateral moments of inertia. The
dot symbol above variables at left side of abovementioned equations is a first derivative
with respect to time. Xb , Yb , Zb are sum of gravity, propulsive, aerodynamic and control
forces and moments in body fixed frame [12]. In similar way, Lb , Nb , Zb are moments
acting on the rocket body.
The main motor thrust force was obtained from the ground tests. The operation time
of the motor was 3 s and the total impulse was 60000 Ns. Thrust force misalignment
angles in pitch HT and yaw planes WT were included into the model.
The aerodynamic characteristics of the rocket were obtained with the aim of the-
oretical and engineering, semi-empirical methods [13, 14]. The nondimensional forces
and moments coefficients were implemented as a function of Mach number, angles of
attack and sideslip into the software using lookup table methodology. Base drag
variations due to main motor state (powered/unpowered phase of flight) were included.
Air thermodynamic properties were used according to U.S. Standard Atmosphere [15].
4 M. Jacewicz et al.
At first the firing logic of the thrusters was developed. The thrusters were assumed to be
fired according to scheme presented in Fig. 2a. Using the single channel control the
actuation system generates forces only over a part of single roll cycle (Fig. 2b). In
result, the lateral thruster must be fired at appropriate roll angle.
Fig. 2. (a) Lateral thrusters firing sequence (b) control force direction as a function of roll angle
The thrusters firings conditions were formulated as follows [12, 16, 17]:
• the lateral thruster was not already consumed
• the time t, from the last thruster firing which has occurred at time tlast is longer than
s 2 ð0; 1Þ
t tlast [ s ð5Þ
• lateral thruster must be fired when is at opposite side than desired flight direction
where c is angle defined in Fig. 2b, is the lateral thruster location angle Ui . sd is
ignition delay (0.001 s) and ssk is half time of the lateral pulse operation time. e is
angular tolerance of firing error and it was chosen to be 2°.
• the pitch angle H of the rocket is smaller than Hg or the time from launch t is
greater than tg
^
H Hg t tg ð7Þ
where xt , yt are target coordinates, xp , yt are predicted impact point location and r2tt
is a small threshold (e.g. 1 m) responsible for guidance accuracy. It prevents from
the undesired consumption of thrusters when the hit error is relatively small (below
several meters) from practical point of view.
Various models could be used to determine the impact point location [9]. Linear
theory based predictive guidance was investigated in [18]. It uses analytical closed form
solution to evaluate the future projectile states. To extend the capabilities of this approach,
modified projectile linear theory intended for fast calculation of the impact point when
projectile is launched at high elevation angles was presented in [19]. The basic constraint
connected with the predictive control is the limited computational power of the onboard
computer. The simple point mass model was used to calculate the impact point of the
rocket. The equations of motion in the navigation On xn yn zn frame were as follows [20]:
dxn
¼ Un ð9Þ
dt
dyn
¼ Vn ð10Þ
dt
dzn
¼ Wn ð11Þ
dt
qffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi
dUn pqd 2 CX0
¼ Un2 þ Vn2 þ Wn2 Un ð12Þ
dt 8m
qffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi
dVn pqd 2 CX0
¼ Un2 þ Vn2 þ Wn2 Vn ð13Þ
dt 8m
qffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi
dWn pqd 2 CX0
¼ Un2 þ Vn2 þ Wn2 Wn þ g ð14Þ
dt 8m
where Un , Vn , Wn are velocities with respect to navigation frame. The air density q was
assumed to be the same as in 6 DOF model and g is gravity acceleration. Drag
coefficient CX0 as a function of Mach number for unpowered phase of flight was used.
The equations of motion are numerically integrated during the rocket flight.
When the target location xt , yt and estimated impact point coordinates xp , yp are
known the range and crossrange errors were computed [5, 21, 22]:
Dx ¼ xp xt Dy ¼ yp yt ð15Þ
No guidance process takes place in the ascending portion of trajectory because the
roll rate of the rocket is too high. The controlled flight starts after vertex of the
6 M. Jacewicz et al.
projectile trajectory and last to the rocket touchdown. The developed mathematical
Author Proof
model of the rocket with the guidance algorithm was implemented in the MATLAB
2018a software. The equations of motion were integrated numerically using ode4
algorithm with fixed step size 1e-4 s. In the simulation the impact points were predicted
with the frequency 10 Hz.
Monte-Carlo methodology and the developed mathematical model of the rocket were
used to obtain the dispersion patterns. The rocket dispersion is a result of combination
of a various errors such as initial aiming inaccuracies, motor thrust misalignments,
launcher vibrations and wind disturbances, which influence the location of impact
point. The thrust force of the single thruster was modelled as a rectangular pulse. The
operation time of the single thruster was fixed and assumed to be 0.02 s. Other pre-
dictive algorithm parameters were: s = 0.5 s, Hg = 0°, tg = 7 s and r2tt = 1 m. To
investigate the influence of the lateral thruster parameters on rocket delivery accuracy a
series of simulations for various force magnitudes was evaluated. The simulation inputs
(Table 1) were assumed to be an Gaussian random independent variables pseudo-
generated with the aim of Marsenne-Twister algorithm. The standard deviations were
assumed using data from [6, 10, 17].
In the each case the 400 simulations was evaluated using desktop PC with Intel i7 and
16 GB RAM. Circular error probable (CEP) was used as a measure of the rocket impact
Monte-Carlo Based Lateral Thruster Parameters Optimization for 122 mm Rocket 7
Author Proof
(a) (b)
(c) (d)
(e) (f)
Fig. 3. The impact point patterns (a) unguided flight and guided flights for thrust force
magnitudes: (b) 200 N (c) 400 N (d) 600 N (e) 800 N (f) 1000 N.
points dispersion. CEP is a radius of circle centered on the target which include 50% of the
impact points. At first, the unguided case was considered. Next, the simulations for
guided flights were evaluated. The obtained hitting patterns were presented in Fig. 3.
8 M. Jacewicz et al.
For unguided flight the lateral dispersion is bigger than longitudinal because of low
Author Proof
launch elevation angle and the CEP is 184.66 m. When the impact point guidance was
applied the CEP values were reduced rapidly. The smallest dispersion was obtained for
the thrust amplitude 1000 N. No bias error in mean impact point location was intro-
duced when the control algorithm was used.
The obtained hitting errors were presented in the form of histograms in Fig. 4.
(a) (b)
(c) (d)
(e) (f)
Fig. 4. The rocket miss distance distributions (a) unguided flight and guided flights for thrust
force magnitudes: (b) 200 N (c) 400 N (d) 600 N (e) 800 N (f) 1000 N.
For the unguided case the miss distances forms a shape similar to Gauss curve.
Several rockets achieved error more than 500 m. For subsequent guided cases a large
number of hit errors below several meters was obtained.
Monte-Carlo Based Lateral Thruster Parameters Optimization for 122 mm Rocket 9
In the Fig. 5 the CEP (vertical axis in logarithmic scale) as a function of lateral
Author Proof
Fig. 5. Circular error probable of the rocket for various thrust force magnitudes
The obtained relationship shows that for values larger than 600 N the profit from
increasing the amplitude is negligible. The obtained results proved that with the aim of
carefully selected lateral thruster parameters and a suitable guidance law, the accuracy
of the projectile might be improved significantly.
5 Conclusion
In this article the influence of solid propellant lateral motor side force magnitude on the
dispersion of 122 mm artillery projectile was presented. 6DoF mathematical model in
Matlab was used to investigate behavior of the rocket. The characteristic feature of
lateral thrusters controlled projectile is its low control authority. Predictive algorithm
based on point mass model was used to steer the rocket precisely to the target. Using
the calculated impact error, the control is applied to reduce the predicted miss distance.
Even relatively small 200 N side force allows to reduce the dispersion significantly.
Numerical Monte-Carlo simulations proved, that with the presented approach the cir-
cular error probable could be reduced approximately ten times when compared to
unguided rocket.
Acknowledgements. This work was supported by The National Centre for Research and
Development (NCBiR) under project DOB-BIO8/10/01/2016 “Projectiles control system tech-
nology development”.
References
1. Pavković, B., Pavić, M., Ćuk, D.: Enhancing the precision of artillery rockets using pulsejet
control systems with active damping. Sci. Tech. Rev. 62(2), 10–19 (2012)
2. Gamble, A., Jenkins, P.: Low cost guidance for the multiple launch rocket system (MLRS)
artillery rocket. IEEE Aerosp. Electron. Syst. Mag. 16(1), 33–39 (2001)
3. Mandić, S.: Dispersion reduction of artillery rockets guided by flight path steering method.
Aeronaut. J. 120(1225), 435–456 (2016)
10 M. Jacewicz et al.
4. Guo, Q.-W., Song, W.-D., Gao, M., Fang, D.: Advanced guidance law design for trajectory-
Author Proof
corrected rockets with canards under single channel control. Eng. Lett. 24(4), 469–477
(2016)
5. Gagnon, E., Vachon, A.: Efficiency analysis of canards-based course correction fuze for a
155-mm spin-stabilized projectile. J. Aerosp. Eng. 29(6) (2016) AQ2
6. Pavkovic, B., Pavic, M., Cuk, D.: Frequency-modulated pulse-jet control of an artillery
rocket. J. Spacecr. Rockets 49(2), 286–294 (2012)
7. Drescher, T., Nielson, J.: Rocket trajectory correction using strap-on GPS guided thrusters.
In: IEEE 1998 Position Location and Navigation Symposium (Cat. No. 98CH36153), Palm
Springs (1998)
8. Jitpraphai, T., Burchett, B., Costello, M.: A comparison of different guidance schemes for a
direct fire rocket with a pulse jet control mechanism. In: AIAA Atmospheric Flight
Mechanics Conference and Exhibit, Guidance, Navigation, and Control and Co-located
Conferences (2001)
9. Fresconi, F., Cooper, G., Costello, M.: Practical assessment of real-time impact point
estimators for smart weapons. J. Aerosp. Eng. 24(1), 1–11 (2011)
10. Gross, M., Costello, M., Fresconi, F.: Impact point model predictive control of a spin-
stabilized projectile with instability protection. In: AIAA Atmospheric Flight Mechanics
(AFM) Conference, Boston (2013)
11. Gagnon, E., Lauzon, M.: Course correction fuze concept analysis for in-service 155 mm
spin-stabilized gunnery projectiles. In: AIAA Guidance, Navigation and Control Conference
and Exhibit, Honolulu (2008)
12. Jitpraphai, T., Costello, M.: Dispersion reduction of a direct-fire rocket using lateral pulse
jets. J. Spacecr. Rockets 38(6), 929–936 (2001)
13. Krasnov, N.F.: Rocket Aerodynamics (NASA Technical Translation). National Aeronautics
and Space Administration (1971)
14. Department of Defense: MIL-HDBK-762 Design of Aerodynamically Stabilized Free
Rockets. United States Department of Defense (1990)
15. National Aeronautics and Space Administration: U.S. Standard Atmosphere 1976,
Washington, D.C. (1976)
16. Gupta, K.: Trajectory correction flight control system using pulsejet on an artillery rocket.
Def. Sci. J. 58(1), 15–33 (2008)
17. Gao, M., Zhang, Y., Yang, S.: Firing control optimization of impulse thrusters for trajectory
correction projectiles. Int. J. Aerosp. Eng. 2015 (2015)
18. Slegers, N.: Predictive control of a munition using low-speed linear theory. J. Guid. Control
Dyn. 31(3), 768–775 (2008)
19. Hainz, L., Costello, M.: Modified projectile linear theory for rapid trajectory prediction.
J. Guid. Control Dyn. 28(5), 1006–1014 (2005)
20. Wang, Z., Chang, S.: Impact point prediction and analysis of lateral correction analysis of
two-dimensional trajectory correction projectiles. Def. Technol. 9(1), 48–52 (2013)
21. Zhang, Y., Gao, M., Yang, S., Fang, D.: Optimization of trajectory correction scheme for
guided mortar projectiles. Int. J. Aerosp. Eng. 2015 (2015)
22. Park, W., Yun, J., Ryoo, C.-K., Kim, Y.: Guidance law for a modern munition. In:
International Conference on Control, Automation and Systems 2010, Gyeonggi-do (2010)
23. Fresconi, F.: Guidance and control of a projectile with reduced sensor and actuator
requirements. J. Guid. Control Dyn. 34(6), 1757–1766 (2011)
24. Burchett, B., Costello, M.: Model predictive lateral pulse jet control of an atmospheric AQ3
rocket. J. Guid. Control Dyn. 25(5), 860–867 (2002)
25. Guo, Q., Song, W., Wang, Y., Lu, Z.: Guidance law design for a class of dual-spin mortars.
Int. J. Aerosp. Eng. 2015, 1–12 (2015)
Author Proof
Book ID : 492180_1_En
Chapter No : 12
Please ensure you fill out your response to the queries raised below
and return this form along with your corrections.
Dear Author,
During the process of typesetting your chapter, the following queries have
arisen. Please check your typeset proof carefully against the queries listed below
and mark the necessary changes either directly on the proof/online grid or in the
‘Author’s response’ area provided below
or and/or
Insert double quotation marks (As above)
or
Insert hyphen (As above)
Start new paragraph
No new paragraph
Transpose
Close up linking characters