Research Article  Open Access
Improvement of Interior Ballistic Performance Utilizing Particle Swarm Optimization
Abstract
This paper investigates the interior ballistic propelling charge design using the optimization methods to select the optimum charge design and to improve the interior ballistic performance. The propelling charge consists of a mixture propellant of sevenperforated granular propellant and onehole tubular propellant. The genetic algorithms and some other evolutionary algorithms have complex evolution operators such as crossover, mutation, encoding, and decoding. These evolution operators have a bad performance represented in convergence speed and accuracy of the solution. Hence, the particle swarm optimization technique is developed. It is carried out in conjunction with interior ballistic lumpedparameter model with the mixture propellant. This technique is applied to both singleobjective and multiobjective problems. In the singleobjective problem, the optimization results are compared with genetic algorithm and the experimental results. The particle swarm optimization introduces a better performance of solution quality and convergence speed. In the multiobjective problem, the feasible region provides a set of available choices to the charge’s designer. Hence, a linear analysis method is adopted to give an appropriate set of the weight coefficients for the objective functions. The results of particle swarm optimization improved the interior ballistic performance and provided a modern direction for interior ballistic propelling charge design of guided projectile.
1. Introduction
Recently, study of the propelling charge design becomes very essential to achieve the interior ballistic performance and assure the safety firing. This study is known as design of the interior ballistic which is considered as a crucial branch of gun system design. There are a number of computerbased interior ballistic models with different capabilities. These models allow the researchers of interior ballistic to predict the interior ballistic performance of a particular gun, charge, and projectile combination. The classic interior ballistics models including the characteristics of the gun, charge, and projectile are utilized to predict the muzzle velocity and the peak pressure. But these models cannot provide the best design of the propelling charge that gives the optimum solution of the interior ballistic performance.
Optimization techniques allow the designers to evaluate a large number of design alternatives in a systematic and efficient manner to find the best design. The solution of an interior ballistic model is considered improved if there is a net increase in muzzle velocity without violating gun constraints. The optimization techniques are used to improve the gun performance and/or decrease the design cost while meeting the constraints appropriate to the problem. Therefore, some researchers tried to apply the optimization techniques with the interior ballistic models. A numerical optimization method called augmented Lagrange multiplier was carried out and coupled to a classic interior ballistic model. This method was used to design the parameters of sevenperforation propellant and to obtain the best muzzle velocity for the projectile fired from 120 mm tank cannons [1]. Another optimization technique was carried out to optimize the ignition system and propelling charge structure in a high muzzle velocity gun [2]. Some researchers of our group tried to use the genetic algorithm (GA) [3]. This algorithm is considered as one of the most currently popular intelligent optimizations to optimize the interior ballistic problem. GA was introduced as a computational analogy of adaptive systems. It is modeled loosely on the principles of the evolution via natural selection, employing a population of individuals that undergo selection in the presence of variationinducing operators such as mutation and recombination (crossover). A fitness function is used to evaluate individuals, and reproductive success varies with fitness. In GA, each genome (chromosome) contains information of the optimization problem by encoding and then should decode the genome information to the optimization problem after finding the optimum solution. According to these above trials of optimization techniques, it is found that they have some limitations in the accuracy of the solution and speed of the convergence. Moreover, the interior ballistic problem cannot be adjusted in a good manner and sometimes the optimum solution cannot be reached.
Particle swarm optimization (PSO) was firstly developed by Kennedy and Eberhart [4]. It was considered as an expansion of an animal social behavior simulation system that incorporated concepts such as nearestneighbor velocity matching and acceleration by distance. The PSO technique becomes very popular due to its ease of implementation and the quick convergence to the optimum solution [5, 6]. Similarly to evolutionary algorithms, PSO exploits a population (a swarm) of potential solutions (particles) which are modified stochastically at each iteration of the algorithm. However, the manipulation of swarm differs significantly from that of evolutionary algorithms, promoting a cooperative rather than a competitive model. More specifically, the other evolutionary algorithms use explicit mutation and selection operators, while PSO uses an adaptable velocity vector for each particle in order to modify the population and favor the best performing individuals. The velocity vector is used to update the position of the particles at each iteration. The particles are moving towards promising regions of the search space by exploiting information springing from their own experience during the search, as well as the experience of other particles. For this purpose, a separate memory is used where each particle stores the best position it has ever visited in the search space.
In this work, PSO technique is coupled with the interior ballistic model to optimize the propelling charge design. The utilized charge is a mixture charge which consists of two different propellants, granular sevenperforated propellant and tubular propellant. PSO technique is applied to the classic interior ballistic model for 76 mm naval gun with guided projectile to improve the interior ballistic performance. Two types of objective functions are used with PSO techniques, singleobjective function and multiobjective function. Through this work, it is found that the optimization results improved the interior ballistic performance and firing safety with a high quality and quick convergence of the optimum solution.
2. Problem Formulation
The optimization of the charge design in interior ballistic process is considered as nonlinear optimization, multiconstrained problem. It consists of multiobjective function and multiple design variables. The design variables have different effects on the performance of interior ballistic. In the interior ballistic process coupled with the optimization technique, there are restrictions between the objective functions. For example, the designer of the interior ballistic intends to obtain the maximum muzzle velocity as a main objective. Maximum muzzle velocity requires a high muzzle pressure, while one of the objective functions is to minimize the muzzle pressure. Hence, the restriction is clear between the muzzle velocity objective function and muzzle pressure objective function. The objective functions should be weighted according to the advantages and disadvantages. The general optimization problem for the classical interior ballistic process can be described as follows: where is the decision vector (design variables), is the multiobjective vector, is the th objective function, and are inequality and equality constraints, respectively, and and are the lower and upper bound of the design variables, respectively.
The relation between the optimization technique and the interior ballistic model is the objective functions that are calculated from the interior ballistic model and then optimized by the optimization model.
There exist some different tools to simulate the interior ballistic process such as lumpedparameter model and twophase flow model. In this work, lumpedparameter model is utilized to simulate the interior ballistic process and calculate the objective values that will be optimized by using the PSO technique. The lumpedparameter model can be written as follows [7]: where is the relative burnt percentage of the th propellant, is the relative burnt thickness of the th propellant, , , are characteristic parameters of the th propellant, is the burning rate coefficient of the th propellant, is the half web thickness of the th propellant, is the pressure in the chamber, is the burning rate pressure index of the th propellant, is the coefficient accounting for the secondary energy losses, is the projectile mass, is the projectile velocity, is the barrel crosssection area, is the tube length, is the ratio of chamber free volume to the bore area, is the impetus force of the th propellant, is mass of the th propellant, and is propellant specific heat ratio.
According to the form shape function of the sevenperforated propellant, should be changed at the start point of slivers separation. The form shape function can be described in (3) as follows:
3. Formulation of Particle Swarm Optimization Algorithm
In the literature there exist numerous methods of optimization that study and analyze the optimization problems under different conditions. PSO is considered as one of the best computationally efficient optimization techniques. It converges to the optimal solution in many problems where most analytical methods fail to converge. PSO has some advantages over other similar optimization techniques, namely, the following [8].(1)PSO has few and simple parameters. Hence, it is easy to implement.(2)It has a more effective memory capability than the GA as every particle remembers its own previous best value as well as the neighborhood’s best value.(3)PSO is more efficient in maintaining the diversity of the swarm, since all the particles use the information related to the most successful particle in order to improve themselves, whereas, in GA, the worse solutions are discarded and only the good ones are saved; therefore, in GA the populations evolve around a subset of the best individuals.The procedure of PSO algorithm can be described as follows.(1)Initialize the swarm by assigning a random position in the problem hyperspace to each particle (the swarm composed of population of random solutions called particles).(2)Evaluate the fitness function for each particle (fitness function is the objective function obtained from the interior ballistic simulation).(3)For each individual particle, particle’s fitness value is compared with its (personal best position). If the current value is better than the value, then set this value as the and the current particle’s position, , as .(4)Identify the particle that has the best fitness value in the swarm. The value of its fitness function is identified as (global best position) and its position as .(5)Update the velocities and positions of all the particles.(6)Repeat steps 2–5 until the condition of stopping is met (e.g., maximum number of iterations or a sufficiently good fitness value).The updating rules of velocity and position can be described as follows: where is the velocity of particle of the th dimension, is the position of particle of the th dimension, indicates the iteration number, , , is the size of swarm, is the dimension of the search space, and are the acceleration coefficients (generally ), rand_{1} and rand_{2} are random numbers , and is inertia weight which is used as a parameter to control the exploration and exploitation in the search space.
4. Application of PSO Algorithm to Interior Ballistic Model
PSO technique is coupled with the lumpedparameter model to improve the interior ballistic performance via optimizing the propelling charge design. This model is applied to 76 mm naval gun with guided projectile utilizing mixed propellant. Two different approaches will be investigated through this present work:(1)singleobjective function PSO method,(2)multiobjective function PSO method.
4.1. SingleObjective PSO Method
Firstly, the singleobjective PSO method is carried out in order to predict the optimum charge design and to improve the performance of interior ballistic process. The singleobjective optimization problem contains objective function, design variables, and constraints. Through this problem the muzzle velocity is considered as the objective function.
4.1.1. Objective Function
The muzzle velocity is considered as the most important parameter in the interior ballistic process. Hence, the only objective function in this method is maximizing the muzzle velocity.
4.1.2. Design Variables
Selection of the design variables is considered very crucial to obtain the optimum solution. Half web thickness and loading density of the propellant are considered the main characteristic parameters of the propellant [3]. In this work, two types of propellants are used. Hence, the design variables are considered as , where and are half web thickness and loading density of the granular propellant and and are half web thickness and loading density of the tubular propellant.
4.1.3. Constraints
The penalty method is utilized to treat the constraints of the optimization problem. The original constrained problem is replaced by a sequence of unconstrained problems [9, 10]. The unconstrained problem consists of the objective function and the penalty terms according to each constraint. Penalty terms penalize the objective function when its corresponding constraint is violated.
The selection of the constraints limits is considered according to the launch safety, the technical requirements, the experimental work for charge design, and the data of the typical gun.
The details of the constraints can be described as follows.
(1) Projectile Muzzle Velocity. The minimum value of muzzle velocity is determined according to the sufficient performance required to satisfy the tactical requirements. The muzzle velocity should be higher than or at least equal to the :
(2) Maximum Chamber Pressure. According to the launch safety requirements, the strength of the gun tube has a limit. The peak pressure must be lower than the maximum pressure value ; this value is considered a design constant for each gun caliber. Hence, the peak pressure should be lower than this value:
(3) Charge Loading Density. The chamber volume is a fixed value. Hence the propelling charge density has a limit. Consider where is considered the maximum charge density limit calculated according to the fixed chamber volume.
The loading density limits depend on the launch safety and the required interior ballistic performance. The limits of this constraint can be considered as follows:
(4) Relative Charge Burnout Point. Relative burnout point is considered as an indicator that explains whether the propellant is completely burnt or not. can be calculated using the following equation: where is the bore length and is the projectile travel inside the bore associated with the complete burning of the propellant.
Low relative burnout point means that the propellant was burnt out early while the projectile was still inside the bore. Hence, some of the kinetic energy of the projectile will be lost to overcome the engraving force. Large relative burnout point means that the projectile exited from the muzzle and the propellant was not completely burnt. Hence, the remaining amount of the propellant is considered as a waste propellant. Due to charge design experience, the limits of this constraint can be considered as follows:
(5) Energy Efficiency of the Charge. Energy efficiency of the charge indicates whether the charge is fully utilized. It explains the ratio between the energy delivered to the projectile and the energy of the charge burning. Low value of means small kinetic energy delivered to the projectile. is calculated using the following equation:
The limits of can be considered as follows:
(6) Half Web Thickness. Values of the half web thicknesses and of the propellant depend on the gun caliber. The limits of and should be selected carefully according to the charge manufacturing, launch safety, and technical requirements. These limits can be considered as follows:
The nonoptimized parameters of the interior ballistic process for the 76 mm gun are tabulated in Table 1.

In this section, the general optimization model can be written in (1) with the design variables :
In the singleobjective problem, two different optimization techniques, GA and PSO, are carried out in conjunction with the lumpedparameter model. The values of the main parameters utilized in GA and PSO algorithm are tabulated in Tables 2 and 3, respectively.


The optimization results for GA and PSO algorithms compared with the experimental data are tabulated in Table 4. The convergence of the two different algorithms for the objective function is illustrated in Figure 1.

From the optimization results illustrated in Table 4 and Figure 1, it is clear that PSO technique is the winner. PSO provides a better performance of interior ballistic propelling charge design. The better performance is represented in the speed of convergence (PSO converges after 170 iterations, while GA converges after 307 iterations) and the solution quality m/s in PSO, while m/s in GA).
Figure 2 represents the convergence of different parameters of the interior ballistic obtained by using PSO technique. These parameters include the design variables and some interior ballistic performance indexes.
4.2. Multiobjective PSO Method
In the actual design process, the designers have different interests in the objective function. Hence, the multiobjective model should be considered to satisfy the different requirements. This section deals with multiobjective optimization model for interior ballistic propelling charge design. Three objective functions, the muzzle velocity , the muzzle pressure , and the energy efficiency of propelling charge are considered to improve the interior ballistic performance and to satisfy the different technical requirements.
The design variables and the constraints are the same as in the singleobjective problem, but the objective function is not the same. The objective function contains the following objectives.
Muzzle Velocity . In gun design, the muzzle velocity is required to be maximized as much as possible to satisfy the tactical requirements. Hence, muzzle velocity will be the key objective function in interior ballistic design.
Muzzle Pressure . The muzzle pressure affects the muzzle blast, the firing precision, and the initial dispersion. Hence, the muzzle pressure should be minimized as possible.
Energy Efficiency of the Charge . The energy efficiency of the charge is considered as one of the main goals in interior ballistic charge design. It should be maximized as possible.
The multiobjective model for the specific gun can be described as follows: where , , and are the weight coefficients of the objective functions determined by the gun designer according to the design requirements. For example, if the designer considers the muzzle velocity to be the main target, the weight coefficient of the muzzle velocity should be increased at the expense of the other coefficients and .
Figure 3 shows the convergence graph for the three objective functions, muzzle velocity, muzzle pressure, and energy efficiency of the charge. It shows a good and rapid convergence to certain values determined by PSO technique.
Figure 4 represents the feasible region of web thickness and propelling charge density with the corresponding objective functions. Three design points marked in Figures 4(a), 4(b), and 4(c) are selected to stand for the maximum muzzle velocity, minimum muzzle pressure, and maximum energy efficiency, respectively. The values of the objective functions according to the three design points are listed in Table 5.

(a) Muzzle velocity
(b) Muzzle pressure
(c) Energy efficiency
Based on the optimization results at point A, Figures 5 and 6 show the pressuretime and velocitytime curves for different schemes, original scheme and optimized scheme. The original scheme refers to the lumpedparameter model described in (2) without optimization. In Figure 5, the optimized scheme decreases the gradient of the pressuretime curve and the muzzle pressure. From Figure 6, the optimized scheme increases the muzzle velocity better than the original scheme. Therefore, the optimized scheme improves the interior ballistic performance and provides a better firing safety.
The feasible region shown in Figure 4 provides a large number of appropriate solutions. Obviously, the optimum solution cannot be found due to the contradictory tradeoffs among the selection conditions. Hence, an appropriate method for adapting the weight coefficients is required for helping the designers of propelling charge to find the optimum solution [11].
4.2.1. Linear Analysis Method
According to the feasible region obtained by the multiobjective function, it is crucial to investigate the effect of changing the weight coefficients on the optimal solution. Hence, the linear analysis method is utilized to obtain a suitable combination of the weight coefficients. Three cases are tabulated in Table 6. One coefficient is changed linearly in each case, and the other coefficients are selected randomly in such a way that , where one has the following.
Case 1. It indicates that is changed linearly, and represents the weight coefficient of muzzle velocity.
Case 2. It indicates that is changed linearly, and represents the weight coefficient of energy efficiency of the charge.
Case 3. It indicates that is changed linearly, and represents the weight coefficient of muzzle pressure.
According to the multiobjective PSO method, the three objective values, , , and are illustrated in Figure 7 to Figure 9 with the different weight coefficients in each case.
Through the observations of Figures 7, 8, and 9, it is obvious to conclude that the weight coefficients have a significant effect on the objective values. The highest muzzle velocity is attained at the highest value of . The change of the energy efficiency has enormous effect on the other objective values. Accordingly, the energy efficiency can be considered as crucial index of interior ballistic performance. It can be seen from the data in Figures 7, 8, and 9 that is strongly recommended to take a value greater than 0.6, to take a value greater than 0.1, and to take a value greater than 0.1. The designers should weigh gains and losses according to the change of the weight coefficients on the objectives to obtain the optimum performance and achieve the tactical requirements and firing safety.
5. Conclusion
This paper has given an account of and the reasons for the widespread use of the particle swarm optimization. PSO algorithm was developed in conjunction with interior ballistic lumpedparameter model utilizing mixed propellant. The mixed propellant is composed of granular sevenperforated propellant and tubular onehole propellant. PSO was applied to optimize the interior ballistic propelling charge design in order to improve the interior ballistic performance. The following conclusions can be drawn from the present study.(1)Based on the characteristics of interior ballistic propelling charge design, particle swarm optimization technique was developed with two different approaches: singleobjective problem and multiobjective problem.(2)GA and PSO techniques were applied to the single objective problem of interior ballistic propelling charge design. The optimization results were compared with the experimental data. The optimized solution showed that PSO technique has a better performance than GA. The better performance was represented in the speed of convergence and the quality of the solution.(3)PSO with multiobjective problem provided a good opportunity for the charge’s designers. The best parameters can be selected through the feasible region to attain the optimum solution according to the tactical requirements. Utilizing the multiobjective method improved the interior ballistic performance and assured the launch safety of the guided projectile.(4)The linear analysis method was developed to present the most appropriate set of the weight coefficients of objectives.
Conflict of Interests
The authors declare that there is no conflict of interests regarding the publication of this paper.
Acknowledgments
The research was supported by the Research Fund for the Natural Science Foundation of Jiangsu province (BK20131348) and the Key Laboratory Fund (Grant no. 9140C300103140C30001), China.
References
 J. R. Gonzalez, Interior ballistics optimization [Ph.D. thesis], Kansas State University, 1990.
 H. S. Zhang, The optimization studies of a new ignition system and charge structure in highvelocity guns [Ph.D. thesis], Nanjing University of Science and Technology, 1998.
 K. Li and X. Zhang, “Multiobjective optimization of interior ballistic performance using NSGAII,” Propellants, Explosives, Pyrotechnics, vol. 36, no. 3, pp. 282–290, 2011. View at: Publisher Site  Google Scholar
 J. Kennedy and R. Eberhart, “Particle swarm optimization,” in Proceedings of the IEEE International Conference on Neural Networks, pp. 1942–1948, December 1995. View at: Google Scholar
 V. Kalivarapu, J.L. Foo, and E. Winer, “Improving solution characteristics of particle swarm optimization using digital pheromones,” Structural and Multidisciplinary Optimization, vol. 37, no. 4, pp. 415–427, 2009. View at: Publisher Site  Google Scholar
 M. Marinaki, Y. Marinakis, and G. E. Stavroulakis, “Fuzzy control optimized by a MultiObjective Particle Swarm Optimization algorithm for vibration suppression of smart structures,” Structural and Multidisciplinary Optimization, vol. 43, no. 1, pp. 29–42, 2011. View at: Publisher Site  Google Scholar
 Y. X. Yuan and X. B. Zhang, Multiphase Hydrokinetic Foundation of High Temperature and High Pressure, Publishing Company of Harbin Institute of Technology, Harbin, China, 2005.
 Y. del Valle, G. K. Venayagamoorthy, S. Mohagheghi, J.C. Hernandez, and R. G. Harley, “Particle swarm optimization: basic concepts, variants and applications in power systems,” IEEE Transactions on Evolutionary Computation, vol. 12, no. 2, pp. 171–195, 2008. View at: Publisher Site  Google Scholar
 W. Y. Sun and Y. Yuan, Optimization Theory and Methods Nonlinear Programming, Springer Optimization and Its Applications, Springer, 2006.
 A. Ravindran, K. M. Ragsdell, and G. V. Reklaitis, Engineering Optimization: Methods and Applications, John Wiley & Sons, Hoboken, NJ, USA, 2nd edition, 2006.
 K. Li and X. Zhang, “Using NSGAII and TOPSIS methods for interior ballistic optimization based on onedimensional twophase flow model,” Propellants, Explosives, Pyrotechnics, vol. 37, no. 4, pp. 468–475, 2012. View at: Google Scholar
Copyright
Copyright © 2014 Hazem El Sadek et al. 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.