Uncertain optimization of in-bore launching performance of artillery based on interval method

Fengjie Xu1 , Guolai Yang2

1, 2School of Mechanical Engineering, Nanjing University of Science and Technology, Nanjing, China

2Corresponding author

Vibroengineering PROCEDIA, Vol. 28, 2019, p. 189-194. https://doi.org/10.21595/vp.2019.21030
Received 18 September 2019; accepted 26 September 2019; published 19 October 2019

Copyright © 2019 Fengjie Xu, 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.
Creative Commons License

In order to study the uncertainty optimization of the in-bore launching performance of larger-caliber artillery, the interior ballistic program was compiled and embedded in ABAQUS finite element software for secondary development, and the dynamic model of in-bore launching was built. The structural parameters of bore, the structural parameters of projectile and parameters of launching propellant were considered, and the uncertainty was described by interval number. BP neural network was used to construct a surrogate model for the dynamic model of in-bore launching. The deterministic transformation of uncertain objective function and uncertain constraints was carried out by using interval order model and interval probability model respectively. The interval uncertainty optimization model of the in-bore launching performance artillery was established by taking the interval radius and midpoint of the projectile muzzle velocity as objective functions and the maximum chamber pressure as constraint. Multi-objective genetic algorithm was used to solve the problem, and the optimal solution and reasonable interval of uncertain parameters were obtained.

Uncertain optimization of in-bore launching performance of artillery based on interval method

  • The analysis model of artillery in-bore launching performance is established
  • The interval method is applied to the uncertain optimization of artillery system
  • Finite element model and internal ballistic subroutine

Keywords: larger-caliber artillery, uncertain optimization, interior ballistic, interval optimization.

1. Introduction

The in-bore launching performance is the key performance of artillery. The motion of projectile in the chamber is affected by many uncertain factors which is the fundamental reason of the difference of ballistic performance. These uncertain factors include: the structural parameters of the gun barrel, the structural parameters of the projectile and parameters of launching propellant.

At present to research in the field of artillery systems rarely involve uncertainty optimization research. Most research focuses on artillery firing dynamics. The research of artillery optimization rarely involves parameter uncertainty. Wang Liqun [1] used the stochastic robust design theory to optimize the firing density, but the research neglected the stage of projectile motion in the chamber. The design variable optimization results obtained by deterministic optimization of the in-bore launching performance are all deterministic parameter combinations, which are not suitable for solving the uncertain problems in actual production. Therefore, this paper, from the point of view of the integrated design of projectile, gun and propellant, carried out the uncertain optimization of the in-bore launching performance. Firstly, the interior ballistic model and finite model element model were combined to construct the gun in-bore firing dynamics model. Then, based on the interval uncertainty optimization theory, the uncertain optimization model of gun bore firing performance is constructed, and the parameter interval of uncertain design variables is obtained by combining BP neural network surrogate model technology and multi-objective genetic algorithm for optimization solution.

2. Artillery in-bore launching model

2.1. Finite element mesh model

The gun barrel, projectile and rotating band were modeled discretized. Except for a few wedge-shaped elements at the front of the rotating band, the hexahedral elements with eight nodes were adopted. The deformation of the gun barrel and the projectile is small, so the mesh of this part is large. Influenced by the driving side force, the contact force between the rotating band and rifling is very large, so the mesh size of rotating band should be small. The mesh model is shown in Fig. 1.

Fig. 1. Finite element mesh model

 Finite element mesh model

2.2. Interior ballistic model

In the process of projectile motion in the chamber, the energy required for motion is provided by the combustion of the propellant. In order to describe the impact of propellant, an interior ballistic amplitude subroutine (VUAMP) is written to provide the driving force for the projectile. The interior ballistic model adopted the classical interior ballistic model based on thermodynamics, and its mathematical model is shown in Eq. (1):

ψ i = χ i Z i 1 + λ i Z i + μ i Z i 2 , d Z i d t = u 1 i p n i e 1 i ,                         i = 1,2 , . . . , n , φ m d v d t = S p , S p l ψ + l = i = 1 n f i ω i ψ i - θ 2 φ m v 2 , d l d t = v ,

where, lψ=l01-i=1nΔiρpi(1-ψi)-i=1nαiΔiψi. i is the number of composite charge, χi, λi and μi are charge form coefficient, m is the mass of projectile, ωi is the mass of charge, S is the equivalent cross-sectional area, φi is the coefficient of second work, l0 is the chamber volume-to-bore area ratio, lψ is the free chamber volume-to-bore area ratio, p is the propellant gas pressure, αi is the co-volume and l and v are the displacement and velocity of projectile respectively.

In the above formula, the burned thickness Z can be solved by the fourth-order Runge-Kutta method to solve the differential equation. The velocity and displacement of the projectile can be directly obtained by the sensor in ABAQUS, and the combination of the two can solve the chamber pressure. The specific execution process is as follows: the initial chamber pressure is obtained through the initial parameters of the interior ballistic, which pushes the projectile forward at the bottom; the velocity and displacement of the projectile can be read directly by the sensor in ABAQUS as the initial conditions for the next calculation to calculate the next chamber pressure, and so on until the projectile is out of the muzzle. At this point, the artillery in-bore launching model was established. This model can be used to calculate the parameters needed for optimization.

3. Uncertain optimization method for artillery in-bore launching performance based on interval optimization

There are generally three uncertain optimization methods: probability model [2], fuzzy model [3] and interval model [4]. The interval model is based on interval mathematics and interval programming. The uncertain variables are regarded as the upper and lower intervals of the uniform distribution. It is not necessary to know the probability distribution of the parameters or the fuzzy membership function. Therefore, this paper chooses the interval method to optimize the uncertainty.

The interval uncertain optimization problem can be expressed as:

m i n f X , s . t .         g i X = , b i I = b L , b R ,           i = 1,2 , . . . , k , X X I = X L , X R = X c , X w ,           X i X i I = X L , X R ,           i = 1,2 , . . . , n ,

where X denotes an uncertainty interval vector, which is composed of an n-dimensional interval vector XI. the superscripts L, R, c and w represent the lower, upper, center, and radius of the interval, respectively. f and g denote the uncertain objective function and uncertain constraints of the interval optimization, respectively. Due to the transitivity of uncertain factors, objective function and constraints also interval numbers rather than real numbers.

3.1. Uncertain design variables

The uncertain design variables considered in this paper include the following three categories (totaling 10). The initial values and ranges of the design variables are listed in Table 1.

a) Projectile parameters: rotating band location ld, rotating band width H, mass eccentricity e0, CM (center of mass) axial location from shell base lR, projectile radius d.

b) Barrel structure parameters: rifling depth t, groove engravings width b, chamber volume W0.

c) Propellant charge parameters: web thickness 2e1, propellant mass ω.

Table 1. Initial values and ranges of design variables

l d (mm)
H (mm)
e 0 (mm)
l R (mm)
2e1 (mm)
d (mm)
t (mm)
b (mm)
W 0 L
ω (kg)

3.2. Treatment of the uncertain objective function

In interval optimization, interval order relations are often used to qualitatively determine whether an interval is better or worse than another interval. In this paper, the interval order relation cw proposed by Jiang [5] was adopted to deal with uncertain objective functions. AIcwBI means that interval B is better than interval A only when the radius and midpoint of interval B are both smaller than interval A. Therefore, the uncertain objective function of Eq. (2) can be transformed into the following multi-objective optimization problem:

m i n X f c X , f w X , f c X = f R X + f L X 2 ,         f w X = f R X - f L X 2 ,

where, fRX and fLX are upper and lower bounds of the objective function respectively.

3.3. Treatment of the uncertain constraints

For the uncertain constraints in Eq. (2), Jiang C. [5] proposed an improved interval probability degree model to compare the two intervals:

P A I B I = 0 ,           A L B R , B R - A L 2 ( A R - A L ) B R - A L B R - B L ,           B L A L < B R A R , B L - A L A R - A L + B R - A L 2 ( A R - A L ) ,           A L < B L < B R A R , B L - A L A R - A L + A R - B L A R - A L B R - A R B R - B L + A R - B L 2 ( A R - A L ) A R - B L B R - B L ,   A L < B L A R < B R , B L - A L A R - A L + B R - A L 2 ( A R - A L ) ,           B L A L < A R < B R , 1 ,           A R < B ,

where PAIBI or PAIBI is the probability of a random variable of interval A being smaller or larger than interval B. The uncertain constraints are transformed into deterministic constraints by the interval probability degree model in Eq. (4).

The inequality constraint of type , such as giIXbiI, can be transformed into the following deterministic constraint:

P g i I X b i I λ i ,

where, λi denotes the interval possibility degree value set by the decision maker. Larger value of λi indicates stricter constraints.

The inequality constraint of type =, such as giIX=biI, can be transformed into the following two deterministic constraint:

P g i I X b i L λ i ,           P g i I X b i R λ i .

Through the above treatment, the interval uncertain optimization problem expressed in Eq. (2) can be converted into the following deterministic optimization problem:

m i n f c X , f w X , s . t .         P g i I X = , b i I λ i ,           i = 1,2 , . . . , k , X X I = X L , X R = X c , X w ,             X i X i I = X L , X R ,             i = 1,2 , . . . , n ,


f c X = f R X + f L X 2 ,           f w X = f R X - f L X 2 .

From the perspective of launch safety and gun performance, the maximum chamber pressure Pmax and muzzle velocity Vel are selected as constraint and objective function respectively. The maximum chamber pressure and projectile velocity can be obtained by finite element model. Thus, the following deterministic multi-objective optimization model is obtained:

m i n - V e l c X , V e l w X , s . t .           P 310 P m a x i I X 340 λ = 0.99 , X X I = X L , X R ,  
V e l c X = V e l R X + V e l L X 2 ,       V e l w X = V e l R X - V e l L X 2 .

The above problem is a two-layer nested optimization problem. The inner optimizer is used to solve the upper and lower bounds of uncertain objective functions and constraints. The outer optimizer is used to search for optimal design variable. Considering that there are two objective functions in Eq. (8), NSGA-II is used to solve the optimal solution. The inner optimizer is the GA. Although the GA has a strong global search capability, its needs too much computation.

In addition, for the artillery in-bore launching model established in this paper, using actual simulation model for nested optimization will bring unacceptable computation cost. In order to solve the contradiction between computational accuracy and computational efficiency, BP neural networks are adopted to replace the initial model. The above method is illustrated in Fig. 2.

Fig. 2. Optimization flowchart

 Optimization flowchart

4. Optimization results

A set of non-dominated solutions organized in a Pareto-optimal front is illustrated in Fig. 3. As can be seen from the Fig. 3, as the midpoint of the velocity interval increases, the radius of the interval also increases, which indicates that blindly pursuing high muzzle velocity of the projectile will lead to large velocity fluctuation. Therefore, in order to give consideration to the high muzzle velocity and robustness, point A in the figure X is selected as the final design scheme. The interval of each design variable is obtained and the results are listed in Table 2.

Fig. 3. Pareto front

 Pareto front

Table 2. Intervals of the design variables

l d (mm)
H (mm)
e 0 (mm)
l R (mm)
2e1 (mm)
X c
X w
d (mm)
t (mm)
b (mm)
W 0 L
ω (kg)
X c
X w

5. Conclusions

In this paper, the interval uncertain optimization theory is utilized to research the artillery in-bore launching performance uncertain optimization problem. Reasonable ranges of projectile, gun barrel and propellant charge parameters are obtained. The midpoint of the interval in the optimization result can be regarded as nominal size, and the interval radius can be regarded as tolerance range. It can be used as a reference in the design of artillery weapon system.


  1. Wang L., et al. An uncertain optimization method for overall ballistics based on stochastic programming and a neural network surrogate model. Engineering Optimization, Vol. 51, Issue 4, 2018, p. 663-679. [Publisher]
  2. Ke H., Su T., Ni Y. Uncertain random multilevel programming with application to production control problem. Soft Computing, Vol. 19, Issue 6, 2014, p. 1739-1746. [Publisher]
  3. Luhandjula M. K. Fuzzy optimization: Milestones and perspectives. Fuzzy Sets and Systems, Vol. 274, 2015, p. 4-11. [Publisher]
  4. Jiang C., Han X., Liu G. P. A sequential nonlinear interval number programming method for uncertain structures. Computer Methods in Applied Mechanics and Engineering, Vol. 197, Issues 49-50, 2008, p. 4250-4265. [Publisher]
  5. Jiang C., et al. A nonlinear interval number programming method for uncertain optimization problems. European Journal of Operational Research, Vol. 188, Issue 1, 2008, p. 1-13. [Publisher]