Mathematical Problems in Engineering

Volume 2016 (2016), Article ID 5910928, 8 pages

http://dx.doi.org/10.1155/2016/5910928

## Aerodynamic Parameter Estimation of a Symmetric Projectile Using Adaptive Chaotic Mutation Particle Swarm Optimization

^{1}National Key Laboratory of Transient Physics, Nanjing University of Science & Technology, Nanjing 210094, China^{2}School of Power and Engineering, Nanjing University of Science & Technology, Nanjing 210094, China^{3}Navy Equipment Research Institute, Beijing 100073, China

Received 27 April 2016; Accepted 7 August 2016

Academic Editor: Zhike Peng

Copyright © 2016 Jun Guan 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.

#### Abstract

This article details a new optimizing algorithm called Adaptive Chaotic Mutation Particle Swarm Optimization (ACM-PSO). The new algorithm is used to perform aerodynamic parameter estimation on a spinning symmetric projectile. The main creative ideas of this new algorithm are as follows. First, a self-adaptive weight function is used so that the inertial weight can be adjusted dynamically by itself. Second, the initialized particle is generated by chaos theory. Last, a method that can be used to judge whether the algorithm has fallen into a local optimum is established. The common testing function is used to test the new algorithm, and the result shows that, compared with the basic particle swarm optimization (PSO) algorithm, it is more likely to have a quick convergence and high accuracy and precision, leading to extensive application. Simulated ballistic data are used as testing data, and the data are subjected to the new algorithm to identify the aerodynamic parameters of a spinning symmetric projectile. The result shows that the algorithm proposed in this paper can effectively identify the aerodynamic parameters with high precision and a quick convergence velocity and is therefore suitable for use in actual engineering.

#### 1. Introduction

Artillery forms an important wing of the army in providing firepower during both war and cross-border skirmishes with an enemy. The traditional artillery shell aeroballistic theory has been developed over several years, and it is an important part of the field of aerodynamics. The accuracy of an artillery shell is an important indicator that is used to estimate its effectiveness. The mathematical models of an artillery shell have often been used for predicting its behavior and flight performance from a conventional approach [1].

The estimation of aerodynamic parameters in the preliminary design stage for vehicles, such as airplanes, missiles, and gun-launched weapons by theoretical methods is very useful. Computational fluid dynamics has recently [2] positively influenced the analytical scenario by providing numerical solutions of the total configuration via sophisticated Euler and Navier-Stokes flow solvers. Several methods based on experimental approaches are essential to validating the analytical result. Although wind tunnel experiments have been used to improve the precision of the estimated parameters, it is difficult to simulate the right flight conditions and weather circumstances. Moreover, the vehicle model tested in the tunnel is slightly different from the actual model used in flight due to last-minute changes [3]. Consequently, an analysis guided by flight data appears to be the best option [4]. Aerodynamic parameter identification is the most developed field in conventional aircraft system identification and has been successfully applied to aircraft and missiles [5]. Suk et al. [6] used maximum likelihood estimation and the extended Kalman filter to identify the system of a UAV in 2003. Tang et al. [7] used a numerically robust least-squares estimator in the frequency domain to identify the aircraft flutter modal parameters in 2008. Burchett [8] used an improved gradient-based method to estimate the aerodynamic coefficients of a projectile from flight range data. Singh and Ghosh [9] used the neural network method to identify the aerocoefficients. Wu et al. [10] designed signals to excite the longitudinal motion of a fly-by-wire passenger airliner to identify the aerodynamic parameters in 2013. However, there have been few studies on spinning symmetric projectiles in recent years. Intelligent identification algorithms have been widely used in the field of aerodynamic coefficient identification with the development of optimization theories [11].

This paper proposes a new PSO algorithm to strengthen the properties of the basic PSO. The new algorithm is then used to identify the aerocoefficients of a spinning projectile. The paper is organized as follows. Section 1 introduces the background of gun-launched projectile aerodynamic parameter identification. Section 2 presents the mathematical modeling of a spinning projectile with 6 degrees of freedom. Section 3 gives the algorithm for the basic PSO and ACM-PSO. Section 4 uses some testing functions to verify the performance of the ACM-PSO. Section 5 uses ACM-PSO to identify the aerodynamic parameters using simulated data. Section 6 is the conclusion.

#### 2. Modeling

To acquire a more precise description of the motion of a spinning projectile, six-degree-of-freedom flight dynamic equations are used in this paper. The detailed expressions of these equations are as follows:

The physical meanings of each variable in (1) and other unlisted equations can be found in [12].

#### 3. Algorithm

##### 3.1. Particle Swarm Optimization (PSO)

Particle swarm optimization [13] is an intelligent optimization algorithm that was proposed by Kenny and Eberhart in 1995. PSO is a population-based optimization tool that could be used to solve several function optimization problems or problems that can be transformed to function optimization problems. Compared with other global optimization algorithms, like Genetic Algorithms or Simulated Annealing, the main strength of PSO is its fast convergence. PSO is inspired by the migration and aggregation of bird flocks as they seek foods. PSO is initialized with a population of random solutions [14]. The potential solutions called particles fly through the problem space by following the current optimum particles. Each particle keeps track of its coordinates in the problem space, which are associated with the best solution achieved so far. This value is called . Another value that is tracked by the particle swarm optimizer is the best value obtained so far by any particle in the neighborhood of the particle. This value is called . The particle also takes the entire population as its topological neighbors, and the best value is a global best, which is called . The position of an individual particle is adjusted according to its own previous searching experience. The best solution is determined by its objective function value. The general procedure of PSO is as follows.

*Step 1 (determination of necessary parameters). *The main parameters of a basic PSO include the population size , particle dimensions , inertial weight , personal cognition coefficient , social cognition coefficient , and maximum iterated generations .

*Step 2 (initialization). *The initialized particles are generated by the methods of randomization. For each particle, both its position and velocity must be initialized.

*Step 3 (position and velocity changes). *The position of a particle is influenced by its velocity. Let denote the position of particle in the solution space at time step . Through adding the velocity as shown in formula (2), the position is updated. The velocity-updating law is described in formula (3):where and represent the particle’s previous speed and position, respectively. and are the particle’s current speed and position, respectively. is the inertial weight coefficient; and are random numbers that are in the range of ; and and are the personal and social cognitive coefficients, which have fixed values in the basic PSO. is the previous individual best position of particle , and is the previous best particle position for all of the particles.

*Step 4 (calculation of the fitness function). *Calculate each particle’s fitness value, and update if the current fitness value of the particle is better than the previous value. Compare the fitness value of each particle with its , and if the current value is better, update . can be updated as shown inwhere denotes a fitness function.

*Step 5 (termination). *Stop the algorithm if the stopping condition is satisfied; if not, go to Step .

##### 3.2. Adaptive Chaotic Mutation Particle Swarm Optimization

The classical PSO is widely used in the area of function optimization, parameter identification, and control system design. PSO has advantages including simple computation and quick convergence but also some disadvantages. For example, long computation times, undispersed initial particles, and easing to fall into local optimum. To address these problems, a modified PSO called Adaptive Chaotic Mutation Particle Swarm Optimization (ACM-PSO) is proposed in this paper. The new algorithm has specific powers, such as short computation time; the generated initial particles will be more dispersed in the solution space and more likely to have a global optimization.

###### 3.2.1. Adaptive Inertial Weight

The inertial weight in (3) is employed to manipulate the impact of the previous history of the velocities on the current velocity [15]. The inertial weight resolves the tradeoff between the global and local exploration abilities of the swarm. Many researchers have advocated that the value of should be large in the exploration (global optimization) state and small in the exploitation (local optimization) state [16]. A large inertial weight encourages global exploration, while a small one enhances local optimization, that is, fine-tuning the current search area. A proper value of the inertial weight supplies the expected balance between the global and local optimization abilities of the particles and consequently improves the effectiveness of the PSO algorithm. It has been shown in many experiments that a large value for the initialized gives a quick convergence to the global exploration in the solution space, and gradually decreases to acquire refined solutions. A new adaptive updating law for the inertial weight is proposed in this paper:where denotes the inertial weight of particle in the generation, is the maximum inertial weight, and is the minimum value. Formula (3) can now be rewritten as the following expression:

The physical meaning of formula (5) can be described as follows. The larger the generation number is, the smaller the inertial weight becomes. The norm denotes the distance between the particle and the global best particle. When the value is small and the value of is also small, it means that this particle is near the global solution, and a small value of enables a refined search to obtain more precise solutions.

###### 3.2.2. Chaotic Initialization

In the PSO algorithm, the initialization of particles is usually randomly generated. The convergence speed of the whole search process and the optimization efficiency of the algorithm are directly affected by the initial position of the spread degree and the uniform characteristic of its position in the search space. It is not likely for the initialization of particles to be a uniform distribution if the solution space has a large dimension or only a small number of particles. Chaos is a common nonlinear phenomenon whose behavior is seemingly complex and random but actually has strong internal rules. Chaos has features of randomness, ergodicity, and sensitivity to initial conditions. An initialization strategy utilizing chaotic theory is used in this paper. Cubical mapping is selected to produce the initial sequence. The mathematical form is as follows:

###### 3.2.3. Premature Judge Criterion

PSO can easily fall into local optima. In this paper, a mutate operator is introduced to the classical PSO to enhance the global optimization ability. When the criterion of premature is satisfied, the mutation operator is applied to the global optimization particle to make the algorithm jump out of the local convergence. A fitness variance is used to judge whether the local convergence has occurred. The mathematical form of the fitness variance is shown inwhere is the average fitness value of all of the particles in the generation. denotes the convergence degree of the algorithm. If the value of is not only smaller than a set threshold but also larger than the smallest theoretical fitness variance, we can judge that the algorithm has fallen into a local convergence. Once the algorithm is in a local convergence, a mutation operator in accordance with a certain probability is applied to : where is Gaussian white noise.

###### 3.2.4. Procedure of ACM-PSO

*Step 1 (determination of relevant parameters). *The main parameters have a population size , the dimension of each particle is , the maximum and minimum values of the inertial weight are and , respectively, the personal cognitive coefficient is , the social cognitive coefficient is , and the maximum iteration generation number is .

*Step 2 (initialization of particles). *Use the method detailed in Section 3.2.2 to generate initial chaotic particles.

*Step 3 (update inertial weight). *According to formula (5), update the inertial weight.

*Step 4 (update position and velocity). *According to formulas (2) and (6), update the position and velocity.

*Step 5 (judge whether the algorithm has fallen into a local convergence). *Use the method detailed in Section 3.2.3 to judge whether local convergence has occurred. If premature convergence has occurred, formula (9) can be used to do a mutation operation to help the algorithm jump out. Otherwise, go to Step .

*Step 6 (termination). *Stop the algorithm if the stopping criterion is met; otherwise, go to Step .

#### 4. Algorithm Validation

To verify the effectiveness of the ACM-PSO proposed in this paper, standard testing functions are used. In this paper, the sphere function , Rastrigin function , and Schaffer function are used as testing functions. The mathematical formulas of these functions are as follows: Sphere function: Rastrigin function: Schaffer function:

The necessary parameters are set as follows: , , , , , , , and threshold of 0.05. The testing result is shown in Table 1.