Parameters identification of a cannon counter-recoil mechanism based on PSO and interval analysis theory

Dan Bao1 , Baolin Hou2

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

2Corresponding author

Vibroengineering PROCEDIA, Vol. 20, 2018, p. 248-253. https://doi.org/10.21595/vp.2018.20243
Received 18 September 2018; accepted 1 October 2018; published 19 October 2018

Copyright © 2018 Dan Bao, 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
Table of Contents Download PDF References
Cite this article
Views 41
Reads 17
Downloads 928
CrossRef Citations 0
Abstract.

In this paper, two methods based on PSO algorithm and interval sequence conversion model or interval analysis theory are proposed to identify two kinds of uncertain parameters of a cannon counter-recoil mechanism during the manual recoil and forward moving process. Before identifying, some test data were obtained during the manual recoil process. Then, the uncertain parameters were described by interval number and a mathematical model about recoil process of a cannon was built. Taking the similarity degree of time-series data as an optimization objective function, Particle Swarm Optimization (PSO) algorithm was used to solve the deterministic optimization problem which was transformed by interval sequence conversion model, and the parameter identification of recuperator in the manual recoil process of a cannon was achieved. On this basis, combining PSO algorithm and Krawczyk algorithm of interval analysis theory, the uncertain parameters of recoil brake were identified. Finally, the results of identification proved that the above-mentioned two methods had a relatively high identification accuracy.

Keywords: counter-recoil mechanism, parameter identification, interval analysis theory.

1. Introduction

The motions of counter-recoil mechanism have a great influence on firing accuracy of a cannon. Therefore, it is significant to identify the relevant parameters of counter-recoil mechanism [1].

The performance parameters of counter-recoil mechanism can be roughly classified into three types: (a) parameters that can be calculated or tested to get specific values; (b) parameters that can be identified to get specific values; (c) parameters that can only be determined the change interval. Generally, decision making is based on the specific values of the objective function and constraint function with respect to each design variable in traditional optimization algorithms. However, when the uncertain parameters are described by interval number, the identification problem cannot be solved by traditional certain identification algorithms because the possible values of the objective function constitute an interval depending on each design variable [2]. This kind of optimization problem often can be solved by converting the optimization problem with uncertain interval into a certain optimization problem by a suitable method.

To solve the parameters identification problem of the counter-recoil mechanism during the manual recoil and forward moving process, a mathematical model was built, and the performance parameters were classified. For parameters that can be identified to get specific values, the parameter identification problem with uncertain interval was converted into a certain single objective optimization problem via an interval number transformation model. The similarity degree between the simulation curve and identification curve was set as the objective function. Then the unknown parameters are identified by particle swarm optimization (PSO). For the parameters that can only be determined their intervals, an interval analysis theory and particle swarm optimization were used to optimize their ranges.

2. Modeling of the cannon counter-recoil process

2.1. Mathematical model

An manual recoil motion was carried out on the counter-recoil mechanism of the cannon. After the manual recoil motion, the recoiling parts of the cannon return to the in-battery position under the action of the force from the recuperator. According to working principle of counter-recoil mechanism, the forces from the recuperator and recoil brake can be expressed as follows:

(1)
F f = A f P f 0 V 0 V 0 - A f x n ,
(2)
F z = k 2 2 ρ A f i A f i + a f a f 2 u 2 + k 1 2 ρ A o f A o f a x 2 u 2 ,           0 < x < λ d T 2 D T 2 - d T 2 , F z = k 2 2 ρ A f i A f i + a f a f 2 u 2 ,         λ d T 2 D T 2 - d T 2 < x < λ .

The motion differential equation of cannon recoil process is:

(3)
M d x d t = F f - F z .

In Eqs. (1-3), Ff is the recuperator force, Fz is recoil brake force, x is the recoil displacement of the cannon, u is the recoil velocity, λ is the maximum recoil length, Pf0 is the initial pressure of the recuperator, n is the gas nonlinear exponential of the recuperator, k1 is the main stream hydraulic resistance coefficient of recoil brake, k2 is the tributary hydraulic resistance coefficient of recoil brake, ax is the area of flow hole, Af is the piston working area of the recuperator, Afi is the gauge working area of the recuperator, af is the groove area, DT is the inner diameter of the drawing cylinder, dT is the outer diameter of the drawing cylinder.

2.2. Parameter classification of the counter-recoil mechanism

According to the mathematical model established above, there are many parameters in the counter-recoil mechanism which can be divided into three types. (a) The first type contains the parameters that are certain and can be determined by calculating or measuring. (b) The second type contains parameters that are unknown or hard to be measured. However, they are certain during a certain action cycle due to the short time of the cycle. This type of parameters can be obtained by identification methods. The two unknown parameters need to be identified in this paper are the initial pressure Pf0 and the gas nonlinear exponential n of the recuperator. (c) The third type contains parameters that are uncertain, and the true values of this type parameters are changing even in a certain action cycle. The present cognitive means can only determine the range of the third type of parameters instead of an exact value or the definite distribution rules. The two interval uncertain parameters are the mainstream hydraulic resistance coefficient of k1 and the tributary hydraulic resistance coefficient k2 of recoil brake.

In order to identify the second and third type parameters above, the entire identification process is divided into two parts. First, identifying Pf0 and n when the uncertain parameters k1 and k2 are changing in their interval ranges. Then, identifying the uncertain parameters k1 and k2 after determining the values of Pf0 and n.

3. Parameters identification of cannon recuperator with PSO algorithm

3.1. Description of the identification problem

The parameters identification problem of the recuperator can be solved by regarding it as an optimization problem: seeking out an optimum Pf0 and n that make the outputs of the mathematic model agree best with the objective function when the uncertain parameters k1 and k2 are changing in their interval ranges. The mathematical model is expected as follows:

(4)
max X f ( X , D ) ,

where, X is the unknown vector that need to be identified, X=(Pf0,  n). D is the uncertain vector, D=(k1,k2). f is the objective function. Because the possible values of D constitute an interval, the results of the objective function is an interval number. In this paper, the interval uncertain identification was turned to a certain optimization problem by interval sequence conversion model, and PSO algorithm was proposed to identify the parameters.

Interval sequence relation is used to evaluate the interval numbers quality. According to the interval number theory, an interval number can be denoted by midpoint Ac and radius Aw:

(5)
A I = A c , A w = A c - A w x A c + A w , x R .

We expect to seek the best unknown parameters that can make the objective function has the largest midpoint and radius. First, the optimization problem with uncertain interval expressed by Eq. (4) is converted into a certain optimization problem. Then, with the linear weighting method, transforming it into a certain single-objective optimization problem [3]:

(6)
max X f d X = 1 - β f c X + ξ ϕ + β f w X + ξ ψ , X = P f 0 , n ,

where, fd is the single-objective evaluation function, β is the weight value with range [0, 1], ξ is a parameter to ensure that fcX+ξ and fwX+ξ are nonnegative, here ξ=0 due to the nonnegative similarity degree. ϕ and ψ are regularization factors of the multi-objective function, they can be obtained by Eq. (7):

(7)
ϕ = min f c X + ξ ,         ψ = m i n f w X + ξ .

Fig. 1. Sensor installation figure

 Sensor installation figure

Fig. 2. Recoil brake structure and pressure test chamber figure

 Recoil brake structure and pressure test chamber figure

In this paper, the similarity degree of the time-series data is used to evaluate the identification results. The objective function is the similarity degree S between the simulated curve or test curve and identification curve and it is comprised of the numerical similarity S1 and the shape similarity S2. The range of S is [0, 1], S= 0 denotes that the two systems are completely alien, S= 1 denotes that the two systems are identical, and 0<S<1 denotes the level of similarity of the two systems [4]. According to multi-attribute and multi-index synthesis theory, when calculating the similarity degree, it is better to use multiplication relation than addition relation between numerical similarity and shape similarity [5]. The numerical similarity is represented by the weighted average of the relative errors and the shape similarity is evaluated by landmark dividing method.

The simulated curve is the pressure of chamber P3 in the recoil brake and it is obtained by virtual model simulation while the test curve is obtained during the manual recoil process using a piezoresistive pressure sensor PN3000.

3.2. Parameters identification results

The unknown parameters of the counter-recoil mechanism are identified by PSO in this paper [6]. During the identification process, the parameters of PSO are as follows: the population size is 50, the maximum iteration step is 80, c1=1.2, c2=1.7, and the inertia weight δ is linearly decreasing, δk=δmin+(δmax-δmin)(Tmax-k)/Tmax, where k is the current iteration step, Tmax is the maximum iteration step, δmax=0.9, δmin=0.4. The range of the parameters needed to be identified are as follows: the range of the initial pressure of counter-recoil mechanism Pf0 is [2.08, 3.12], and the gas nonlinear exponential of counter-recoil mechanism n is [1.04, 1.56]. The range of the uncertain parameters k1 and k2 are [1.2, 1.76] and [1.8, 2.64].

Table 1 and Table 2 give the identification results of the simulation data and test data. Comparing with the theoretical values, the mean error of the identification results of two simulations are 0.360 % and 0.865 %. we can see that the identification accuracy is relatively high, and the model can be used for parameters identification of the real system.

Table 1. Identification results of simulation curves

Parameter
Theoretical value
Identification value 1
Identification value 2
Similarity interval
P f 0 / MPa
2.6
2.6108
2.6079
[0.9136, 0.9927]
n
1.3
1.3136
1.3089
[0.9023, 0.9891]

Table 2. Identification results of test curves

Parameter
Theoretical value
Test value 1
Test value 2
Similarity interval
P f 0 / MPa
2.6
2.5735
2.5591
[0.6080, 0.6536]
n
1.3
1.2878
1.2857
[0.6171, 0.6522]

Fig. 3. Pressure curves identified by mathematical model

Pressure curves identified by mathematical model

Fig. 4. Pressure curves identified by test data

Pressure curves identified by test data

4. Parameters identification of cannon recoil mechanism with interval analysis theory

4.1. Description of the identification problem

The mathematical model used for identification is still the mathematical model of the manual recoil process of the cannon established above. In this section, the two parameters to be optimized are respectively the mainstream hydraulic resistance coefficient of recoil brake k1 and the tributary hydraulic resistance coefficient of recoil brake k2. The basic idea of interval analysis theory is to represent the real value by an interval including real value, and then to narrow the range of the interval by a series of interval operations. For the third type of parameters k1 and k2, a method combining PSO algorithm and interval analysis theory Krawczyk algorithm was proposed for identification. The identification process is as follows:

1) Giving the initial interval of uncertain parameters and sampling within the interval. Then the sampling values are substituted into the model to obtain the corresponding output, and the relative error between it and the target value is calculated.

2) Selecting the corresponding sampling value with the minimum relative error in the first step and taking it as the midpoint of the new interval. The middle point is subtracted from the upper and lower bounds of the original interval, and the result with small absolute value is taken as the radius of the new interval to generate the new interval of uncertain parameters.

3) If the midpoint calculated in the second step happens to be the midpoint of the original interval, reducing the original interval by 10 % and returning to the first step.

4) Repeat steps 1-3 until the relative errors of all the output values meet the accuracy requirement.

On the basis of identification in the above section, after determining the values of Pf0 and n, identifying the uncertain parameters k1 and k2 according to the above steps. The initial value of k1 and k2 are set to [1.2, 1.76] and [1.8, 2.64]. After setting the value of each parameter in the model, the iterative calculation can be conducted according to the Krawczyk algorithm process. After each iteration, choose new sample within the new range generated. The simulation is carried out according to the sample values and calculating the similarity degree of chamber P3 pressure between simulated data and target data. Finally, the iterative process is terminated according to whether the similarity meets the requirements.

4.2. Parameters identification results

Table 3 and Table 4 give the identification results of the simulation data and test data. The results of identification are evaluated by the uncertainty level of interval γAI=Aw/Ac and relative error.

Table 3. Identification results of simulation curves

Parameter
Theoretical value
Identification value
Relative error / %
γ
k 1
1.4
[1.3468, 1.5617]
3.875
0.0739
k 2
2.2
[2.0693, 2.4821]
3.441
0.0907

Table 4. Identification results of test curves

Parameter
Theoretical value
Test value
Relative error / %
γ
k 1
[1.3659, 1.6908]
0.1063
k 2
[1.9856, 2.4127]
0.9711

5. Conclusion

1) In this paper, a mathematic simulation model of a cannon counter-recoil mechanism during the manual recoil and forward moving process was built. Because the cannon weapon is a complicated electro-mechanical and hydraulic integrated system with multi-parameters, in order to use different identification algorithms for different types of parameters, the performance parameters are classified into three types.

2) The parameters of a cannon recuperator may be considered as the second type parameters that can be identified their specific values. The parameter identification problem with uncertain interval is converted into a certain optimization problem via an interval number transformation model. The similarity degree of time-series data between the simulation curve and identification curve is set as the objective function. Then the unknown parameters are identified by particle swarm optimization (PSO).

3) The parameters of the cannon recoil brake are the third type parameters that are changing during a certain firing process. On this basis of the first part identification, combining PSO algorithm and Krawczyk algorithm of interval analysis theory, the uncertain parameters of recoil brake were identified. Finally, the results of identification prove that the above-mentioned two methods have a relatively high identification accuracy.

References

  1. Hou B. L., Qiao J. M., Liu C. M. Ammunition Automatic Loading System of Howitzer. Publishing of Ordnance Industry, Beijing, 2010. [Search CrossRef]
  2. Jiang C. Uncertainty Optimization Theory and Algorithm Based on Interval Number. Hunan University, Changsha, 2008. [Search CrossRef]
  3. Sun X. H., Ding X. X. Research on multi-objective topology optimization design methods for structure. Machine Design and Research, Vol. 28, Issue 4, 2012, p. 1-4. [Search CrossRef]
  4. He J. Research on Validity Evaluation and Parameter Calibration of Generator Model Based on Hybrid Dynamic Simulation. Chongqing University, Chongqing, 2011. [Search CrossRef]
  5. Xu Z. Uncertain Multi-Attribute Making: Methods and Applications. Springer, 2015. [Search CrossRef]
  6. Kennedy J., Eberhart R. Particle swarm optimization. IEEE International Conference on Neural Network, Perth, Australia, 1995, p. 1942-1948. [Publisher]