Next Article in Journal
Polymer Electrolyte Fuel Cell Degradation Mechanisms and Their Diagnosis by Frequency Response Analysis Methods: A Review
Next Article in Special Issue
Passive Heating and Cooling of Photovoltaic Greenhouses Including Thermochromic Materials
Previous Article in Journal
Frequency-Domain Modeling of Harmonic Interactions in Voltage-Source Inverters with Closed-Loop Control
Previous Article in Special Issue
Hydrothermal Investigation of a Microchannel Heat Sink Using Secondary Flows in Trapezoidal and Parallel Orientations
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Numerical Calculation and Uncertain Optimization of Energy Conversion in Interior Ballistics Stage

1
School of Mechanical Engineering, Nanjing University of Science & Technology, Nanjing 210094, China
2
Inner Mongolia North Heavy Industries Group Corp. Ltd., Baotou 014033, China
*
Author to whom correspondence should be addressed.
Energies 2020, 13(21), 5824; https://doi.org/10.3390/en13215824
Submission received: 23 August 2020 / Revised: 26 October 2020 / Accepted: 2 November 2020 / Published: 7 November 2020
(This article belongs to the Special Issue Modelling of Thermal and Energy Systems)

Abstract

:
Gun firing is a process that converts propellant chemical energy to projectile kinetic energy and other kinds of energies. In order to explore the energy conversion process, firstly, the interior ballistics mathematical model and the barrel-projectile finite element model are built and solved. Then, the related variable values and energy values are obtained and discussed. Finally, for improving energy efficiency, the interval uncertainty optimization problem is modeled, and then solved using the two-layer nested optimization strategy and back-propagation (BP) neural network surrogate model. Calculation results show that, after optimization, the heat efficiency raises from 31.13% to 33.05% and the max rifling stress decreases from 893.68 to 859.76 Mpa, which would improve the firing performance and prolong the lifetime of the gun barrel.

1. Introduction

A gun is a kind of mechanical device which is used to fire projectiles through the gun firing system. The firing process of a gun can be briefly described as follows: firstly, the propellant combusts with generation of propellant gas, and then the propellant gas propels the projectile until it moves out of the gun barrel. This integral process is also known as the interior ballistics process. Besides movement of the projectile and combustion of the propellant, another three physicochemical phenomena also happen in the interior ballistics stage, namely friction and plastic deformation of the rotating band, recoil of the gun barrel and heat transfer of the gun barrel.
Many researchers have studied these physicochemical phenomena by numerical and experimental approaches. For the research on propellant combustion and the interior ballistics model, Liao et al. [1] established an innovative interior ballistics model which adopts actual combustion characteristics between propellants and pressure impulse. Monreal-González et al. [2,3] developed a program named UXGun to solve the one-dimensional interior ballistics model. They also compared numerical results with test date and the results showed a good agreement. An interior ballistics model based on conserved variables was built by Otón-Martínez et al. [4] and calculation results were discussed by comparing with the model of Gough. A two-phase flow model of interior ballistics which contains igniters was given and four experiments were carried out by Jaramaz et al. [5] for investigating the interior ballistics process. To investigate detonation of granulated solid propellants, López et al. [6] used Rusanov and MacCormack-TVD (short for total variation diminishingschemes) numerical schemes to solve combustion equations and compared with calculation results by these two schemes. For the research on the movement of a projectile inside the gun barrel, Wu et al. [7,8,9,10,11] performed a series of engraving tests by employing a gas gun and various kinds of rotating bands. He discussed the effect of shape, diameter, material and processing technology of rotating bands on engraving characteristics. Alexander et al. [12] described a finite element model of the firing system of a 155 mm gun. Solving this model by Abaqus software, he analyzed velocity, acceleration and displacement of the projectile. Sudarsan et al. [13] built a dynamics model of a 155 mm gun which contains a projectile gun barrel and rotating band, and then he conducted sensitivity analysis of the dependent parameters of shot start pressure. The influence of gun–projectile interactions as well as projectile- and gun-dependent factors on the dispersion were investigated by Dursun et al. [14]. Carlucci et al. [15] employed a simplified model to study the dynamic response of the projectile. His results indicated that accelerometer data should be used with caution. Ding et al. [16] proposed a mesh generation method for building a finite element model of the gun barrel. Moreover, he explored the relationship between a worn gun barrel and projectile velocity. For the research on heat transfer of the gun barrel, Zhang et al. [17] employed a method which combines the finite difference method with the sequential function specification method. His method can predict the heat flux of the gun barrel and has high accuracy by making comparisons of numerical tests. Considering material property to vary with temperature, Mathematica software was used by Hill et al. [18] for solving the one-dimensional transient heat transfer numerical model. In order to study the self-ignition phenomenon, Işık et al. [19] and Değirmenci et al. [20] used the finite element model of a 7.62 mm rifle to get temperature distribution of the gun barrel and discussed the effect of time interval, grain sizes and initial temperatures parameters. Talaee et al. [21] discussed the effect of different shooting rates in consecutive shooting on temperature distribution of the gun barrel. For the research on recoil of the gun barrel, a new mechanism is proposed by Sefiane et al. [22] to explain the occurrence of boiling crisis and it could improve recoil instability. Kathe and other researchers [23,24] suggested the concept of the soft recoil mechanism accelerating first in the forward direction in the armament technology. Mitchell et al. [25] studied the dynamics of the recoil system when firing projectiles with Mach 4.4 muzzle velocity. Hassaan et al. [26] investigated the dynamics of the barrel assembly-recoil mechanism of military cannons when using a hydraulic damper and a constant stiffness helical spring in their recoil mechanisms and evaluated the performance through the minimum and maximum displacements and the settling time of its response upon firing.
Heat machine is defined as a mechanical device which converts chemical energy to kinetic energy. From the energy point of view, guns also belong to the category of heat machine. Heat efficiency refers to the ratio of output kinetic energy to chemical energy of the heat machine. In engineering, heat efficiency is widely studied, such as heat efficiency of the gas turbine system [27], gasoline and ethanol mixture engine [28] and diesel engine [29]. Moreover, many methods have been employed to improve heat efficiency, such as optimizing structure and improving reliability [30,31]. For the gun, heat efficiency is also an important performance indicator, because higher heat efficiency means longer range and more penetration capability of the projectile. However, many researchers have built a simplified model and studied one or two physicochemical phenomena for getting several variable values, such as stress and velocity. Little attention has been paid to energy conversion in the interior ballistics stage.
In this paper, an energy calculation model of a fixed gun firing system is established and solved. Subsequently, energy values are obtained, and heat efficiency is calculated. Finally, the uncertain optimization is carried out for higher heat efficiency. This paper would provide an approach to investigate energy conversion of a fixed gun and would also be beneficial for investigating other kinds of thermal propulsion systems.

2. Energy Calculation Model of Firing System of Fixed Gun

Concerning the gun firing system, the propellant, projectile and gun barrel are three basic components. According to the function of each component, an integrated energy calculation model of the firing system of a fixed gun is divided into two sub-models, namely the interior ballistics mathematical model and the barrel-projectile finite element model. Additionally, the fixed gun ignores movement of the gun bracket during the firing process. Table 1 shows some parameters of the propellant and mechanical structure of these two sub-models.

2.1. Interior Ballistics Mathematical Model

Before the propellant starts to combust, it is placed inside the gun barrel from breech to projectile bottom, as shown in Figure 1. When the propellant combusts, its products are gas-phase propellant and solid-phase propellant particles.
The interior ballistics mathematical model describes the law of propellant combustion and product’s state from the moment the propellant is ignited to the moment the projectile leaves the muzzle. Based on multiphase fluid mechanics theory and combustion theory, a modified one-dimensional two-phase flow of the interior ballistics model is established. In addition, for establishing this model, some hypotheses are stated below:
(1)
The gas-phase propellant and solid-phase propellant have independent element volume, and they are considered as homogeneous continuum.
(2)
The appropriately defined average of flow properties is used to describe the heterogeneous compound flow, which consists of two interacting continua.
(3)
The solid-phase propellant is incompressible and its density value is fixed.
(4)
The Nobel-Abel state equation is used to describe the gas-phase state.
(5)
The geometrical combustion theory and exponential combustion theory are used to describe the law of propellant combustion.
This model takes the characteristic that the cross-section area of the barrel bore varies with its distance from the breech into consideration. Firstly, five main equations are given.
Mass conservation equation of gas-phase:
( φ ρ g A ) t + ( φ ρ g u g A ) x = m ˙ c A
where φ is the porosity, ρ g is the gas-phase propellant density, u g is the gas-phase propellant velocity, A is the cross-section area of gun bore and m ˙ c is the generation mass of propellant gas per unit time.
Momentum conservation equation of gas-phase:
( φ ρ g u g A ) t + ( φ ρ g u g 2 A ) x + ( A φ ) p x = f s A + m ˙ c u p A
where p is the propellant gas pressure, u p is the solid-phase propellant velocity and f s is the interphase resistance.
Energy conservation equation of gas-phase:
( φ ρ g u g A ( e g + u g 2 2 ) )   t + ( φ ρ g u g A ( e g + p ρ g + u g 2 2 ) )   x + p ( A φ )   t = Q p A f s u p A + m ˙ c A ( e p + p ρ p + u g 2 2 )
where Q p is the interphase heat transfer, e g is the internal energy of gas-phase propellant, e p is the potential energy of solid-phase propellant and ρ p is the solid-phase propellant density.
Mass conservation equation of solid-phase:
[ ( 1 φ ) ρ p A ]   t + [ ( 1 φ ) ρ p   u p A ]   x = A   m ˙ c
Momentum conservation equation of solid-phase:
[ ( 1 φ ) ρ p u p A ]   t + [ ( 1 φ ) ρ p u p 2 A ]   x + A ( 1 φ )   p   x + [ ( 1 φ ) R p A ]   x = f s A m ˙ c u p A
where R p is the particle stress.
To solve these equations, it is necessary to add another seven auxiliary equations, which are given in the following.
Gas-phase state equation:
p ( 1 ρ g β ) = R T
where R is the gas constant, T is the propellant gas temperature and β is the co-volume.
Particle stress equation:
R p = { ρ p c 2 φ φ   0 1 φ φ φ   0 , φ φ   0 ρ p c 2 2 k ( 1 φ ) [ 1 e 2 k ( 1 φ ) ] φ   0 < φ < φ * 0 φ φ *
c ( φ ) = { c 1 φ 1 φ φ φ 0 c 1 e k ( φ φ 0 ) φ 0 < φ < φ * 0 φ φ *
where c 1 is the sound velocity of the solid-phase particle when φ   =   φ 0 , k is the stress attenuation factor, φ 0 is the initial porosity constant by natural packing of the propellant and its value is 0.45. Additionally, φ *   =   φ 0   +   0.151 .
Interphase resistance equation:
f s = 1 φ d p | u g u p | ( u g u p ) ρ g × { 1.75 φ φ 0 1.75 ( 1 φ 1 φ 0 φ 0 φ ) 0.45 φ 0 < φ < φ 1 0.3 φ φ 1
φ 1 = [ 1 + 0.02 1 φ 0 φ 0 ] 1
where d p is the equivalent diameter of the propellant particle.
Interphase heat transfer equation:
Q p = ρ p ( 1 φ ) S p q / M p
where S p is the superficial area of the propellant particle, M p is the mass of the propellant particle and q is the heat flux coefficient between two phases.
Temperature equation of particle surface:
d T p s d t = 2 q λ p a p π t + δ t t δ t
where T p s is the temperature of the particle surface, a p is the temperature conductivity coefficient and λ   is the thermal conductivity coefficient. Here, the temperature principle is adopted as combustion physical condition. Specifically, the initial temperature of the particle is 293 K. If the temperature of the particle surface is higher than 615 K, which is the combustion temperature, then the propellant begins to combust.
Combustion rate equation:
d r p d t = d ˙
d ˙ = b   p n
m ˙ c = d ˙ ( 1 φ ) ρ p 2   S p / M p
where r p is the equivalent radius of the propellant particle, d ˙ is the linear combustion rate in porous bed, b is the combustion rate coefficient and n is the combustion rate exponent.
Projectile motion equation:
d ν d t = A d p d f m
where ν   is the projectile velocity, A d is the projectile bottom area, p d is propellant gas pressure at projectile bottom, f is projectile motion resistance and m is projectile mass.
In order to solve the equations above, the MacCormack difference scheme which has second order accuracy is adopted, as shown in Equations (13a)–(13c). Fortran language is used to develop a solving program called TFIBS (short for two-phase flow interior ballistics solver).
U ¯ j n + 1 = U j n Δ t Δ x ( F j + 1 n F j n ) + Δ t H j n
U ¯ ¯ j n + 1 = U ¯ j n + 1 Δ t Δ x ( F ¯ j n + 1 F ¯ j 1 n + 1 ) + Δ t H ¯ j n + 1
U j n + 1 = 1 2 ( U j n + U ¯ ¯ j n + 1 )

2.2. Barrel-Projectile Finite Element Model

The barrel-projectile finite element model contains four parts: gun barrel, rifling, projectile and rotating band. This model, which is established by Hypermesh software, is shown in Figure 2, with partially hidden gun barrel and riflings for better visualization. Related material parameters are listed in Table 2. Totally, there are 1,239,092 elements and 1,421,076 notes in the model. In addition, element sizes of the gun barrel, rifling and projectile are approximately 5 mm. As for the rotating band, considering the accuracy of the engraving process, its meshing size is approximately 0.5 mm. The contact algorithm between rotating band and bore surface is set as general contact in Abaqus software. The dynamics and partial energy values are calculated by Abaqus solver.

2.3. Calculation Process of Energies

Figure 3 depicts the calculation process of energies. Firstly, we solve the interior ballistic model by TFIBS program and get the related variable values of propellant gas. Then, we divide the calculation process into two routes. One is to integrate the values of pressure and porosity to calculate energy Q1 and E4. Another one is adopting Abaqus subroutine Vuamp and Film to set the values of pressure and temperature as initial load and boundary condition in the finite element model, and then to use Abaqus to calculate energy E1, E2, E3 and E5.

3. Calculation Results and Analysis of Energy Conversion

As mentioned above, several physicochemical phenomena happen in the interior ballistics stage, and appearance of these phenomena means occurrence of energy conversion. From the perspective of energy, the firing process can be described as: firstly, chemical energy of propellant converts into propellant gas energy through combustion; then, the propellant gas energy, as intermediate energy, converts into projectile kinetic energy, thermal energy of gun barrel and recoil energy, etc. According to the law of energy conservation, these energies can be given as:
Q = Q1 + E1 + E2 + E3 + E4 + E5 + E6
where Q is the chemical energy of total propellant, Q1 is the chemical energy of the uncombusted propellant, E1 is the translational kinetic energy of the projectile, E2 is the rotational kinetic energy of the projectile, E3 is the sum of frictional energy and plastic deformation energy, E4 is the recoil energy, E5 is the thermal energy of the gun barrel and E6 is the propellant gas energy. Among them, Q is a constant while Q1, E1, E2, E3, E4, E5 and E6 are time-varying. The energy conversion process is displayed in Figure 4 and will be explained in this section.

3.1. Numerical Calculation of E1 and E2

As is exhibited in Figure 4, combustion of the propellant produces propellant gas, which is filled in the space from breech to projectile bottom and pushes the projectile to do translational and rotational movement inside the gun barrel. At the same time, the projectile gets translational kinetic energy and rotational kinetic energy. E1 and E2 can be expressed by:
E 1 = 1 2 m v 2
E 2 = 1 2 j w 2
where m is the projectile mass, v is the translational velocity of the projectile, j is the rotational inertia of the projectile and w is the rotational velocity of the projectile. m and j are constants and their values are listed in Table 1. Solving the interior ballistics model by the TFIBS program, propellant gas pressure is obtained and its data at some time points are shown in Figure 5a. In addition, Figure 5b shows propellant gas pressures at the projectile bottom.
By loading the pressure shown in Figure 5b to projectile bottom in the finite element model with adoption of Abaqus subroutine Vuamp, translational velocity and rotational velocity of the projectile can be obtained by Abaqus, as shown in Figure 6a,b. Obviously, curves of the projectile translational velocity and rotational velocity are highly similar. This is because the entanglement angle of the rifling is fixed, so that the translational velocity is proportional to the rotational velocity. E1 and E2 are calculated by Equations (15) and (16). Their values are shown in Figure 6a,b, respectively. From Figure 6a, the muzzle translational velocity is 909 m∙s−1. Table 3 gives muzzle translational velocity data of six firing tests, which are obtained by high-speed photography. Comparing E1 obtained by numerical simulation with E1 measured by firing tests, the error is −2.58%.
As shown in Figure 6a, duration from the moment the propellant starts combusting to the moment the projectile leaves the muzzle is 18.9 ms, and this process is the interior ballistics stage. According to projectile bottom pressure, the generation period of E1 can be divided into three periods. The time between 0 and 3.6 ms belongs to the first period. In this period, the propellant combusts apart, and the propellant at the projectile bottom has not yet combusted. Hence, the projectile bottom pressure is 0 MPa and the projectile is static. Taking the pressure data at 3.2 ms in Figure 5a as an example, the total propellant length is 1.08 m. From 0 to 0.7 m, the propellant gas pressure is greater than 0 Mpa, and the propellant gas pressure from 0.7 to 1.08 m is equal to 0 Mpa. It means that the propellant only combusts to 0.7 m and the propellant from 0.7 to 1.08 m has not yet started to combust. The time of 3.6–4.7 ms belongs to the second period. During this period, the projectile bottom pressure gradually increases from 0 to 25.6 Mpa. Meanwhile, the static motion resistance of the projectile also increases and is equal to it. So, the projectile is still stationary. The third period is from 4.7 to 18.9 ms. During this period, as the projectile bottom pressure continues to increase, the projectile begins to move, and the translational velocity, rotational velocity, translational kinetic energy and rotational energy of the projectile increase with time.

3.2. Numerical Calculation of E3

The rotating band, which is ring-shaped and surrounds the projectile, moves along with the projectile inside the gun barrel. By partially hiding the gun barrel and riflings for better visualization, Figure 7 presents the size relationship between rotating band diameter (D1) and gun bore diameter (D2) before the projectile moves. Apparently, contact will occur between the rotating band, riflings and inner surface of the gun barrel after the projectile moves. The purpose of this structural design is to make the projectile rotate and prevent leakage of propellant gas. In addition, the type of projectile in this paper has two rotating bands.
Again, the finite element model and the projectile bottom pressure in Section 3.1 are adopted for calculation, and Abaqus is used to output the Mises stress and morphology of the rotating band and riflings at 6.2 and 9.5 ms, as shown in Figure 8 and Figure 9. According to the characteristics of energy, the generation period of E3 can be divided into two distinct parts. The engraving period is from 4.7 to 8.2 ms. During this period, the rotating band is squeezing into the riflings and is subjected to the cutting and squeeze of the riflings. When the Mises stress experienced by the rotating band is greater than the yield limit of its material, the rotating band will undergo plastic deformation and form rotating band grooves, which are shown in Figure 8. Meanwhile, the rotating band also produces friction with riflings and the inner surface of the gun barrel. The translation-rotation period is from 8.2 to 18.9 ms. During this period, rotating band grooves have been completely formed, and the rotating band no longer undergoes plastic deformation. The shape of the rotating band is shown in Figure 9. However, the rotating band still produces friction with riflings and the inner surface of the gun barrel. In conclusion, both frictional energy and plastic deformation energy generate in the engraving period, while only frictional energy generates in the translation-rotation period.
The frictional energy and the plastic deformation energy are calculated by Abaqus energy output functions ALLFD (short for all frictional dissipated energy) and ALLPD (short for all plastic dissipated energy), respectively. The value of E3 is shown in Figure 10. From 8.2 to about 15 ms, the increase rate of E3 is less than its increase rate during the engraving period, which is due to the lack of plastic deformation energy during this period. From about 15 to 18.9 ms, the average increase rate of E3 is greater than its average increase rate during the engraving period, because after 15 ms, the value of projectile velocity is increasing and relatively high, which makes the projectile displacement higher in unit time.

3.3. Numerical Calculation of Q and Q1

Q, Q1 and b can be calculated by the following Equations ((17)–(20)):
Q = k p m p
m 1 = 0 x A ( 1 φ ) ρ p d x
Q 1 = k p m 1
b = Q 1 Q
where k p is the coefficient of heat value, m p is the propellant mass and ρ p is the propellant density, and the values of these three constants are listed in Table 1. The calculated Q is 6.05 × 107 J. φ is the porosity and is obtained by the TFIBS program. The data of porosity at some time points are exhibited in Figure 11a. m 1 is the uncombusted propellant mass. b is the percentage of uncombusted propellant energy, and its value is shown in Figure 11b.
As can be seen from Figure 11a, the porosity tends to increase roughly before 16.4 ms. This indicates that the propellant is combusting and producing propellant gas. At 16.4 ms, the porosity value of all spaces inside the gun barrel is 1, which means that the propellant has completely combusted and there is no solid-phase propellant. In other words, from 16.4 to 18.9 ms, no new propellant gas is produced, and the movement of the projectile depends on the propelling effect of the existing propellant gas. This propelling force is greater than the frictional resistance experienced by the projectile, so the projectile kinetic energy is still increasing. The absolute value of the slope of the curve in Figure 11b represents the combustion rate of the propellant, and the combustion rate generally shows a trend of increasing first and then decreasing. Referring to Equation (11b), this is because the combustion rate of the propellant is a monotonous increasing function of the propellant gas pressure, and the propellant gas pressure first increases and then decreases in the interior ballistics stage.

3.4. Numerical Calculation of E4 and E5

As is exhibited in Figure 4, the propellant gas also exerts force on the breech, which propels the gun barrel to move with a direction opposite to the projectile movement. p b is the propellant gas pressure at the breech and is obtained by the TFIBS program, as shown in Figure 12a. Recoil force, F z , generates from the gun recoil mechanism and can be expressed by Equation (21). By loading p b and F z on the breech of the finite element model in Section 3.1, L is the gun barrel displacement and is obtained by Abaqus, as shown in Figure 12a. Finally, E 4 can be calculated by Equation (22) and its value is shown in Figure 12b.
F z = K 1 ρ Z 2 ( A 0 A p ) 3 α x 2 v l 2 + K 2 ρ Z 2 A f j 3 Ω 1 2 v l 2
E 4 = 0 t A b p b L d t
where K 1 is the hydraulic resistance coefficient, Ω 1 is the minimum cross-sectional area of the hydraulic branch, α x is the area of fluid orifice, A f j is the cross-sectional area of piston, v l is the recoil velocity and A b is the breech area.
As is exhibited in Figure 4, due to the temperature difference between the propellant gas and the gun barrel, heat transfer occurs at the contact surface between them. This process causes the gun barrel temperature to increase and the gun barrel to obtain thermal energy. The average temperature of propellant gas calculated by the TFIBS program is shown in Figure 13a. By using Abaqus subroutine Film to apply the average temperature value of propellant gas as the boundary condition to the inner surface of the gun barrel, the temperature distribution is calculated by Abaqus, as shown in Figure 14. Furthermore, E5 is calculated by Abaqus energy output function ALLIE (short for all internal energy), and the result is shown in Figure 13b.
At 18.9 ms, the ratio between E4 and Q is 0.56%, and the ratio between E4 and Q is 1.18%. The values of both E4 and E5 are relatively low. This is the reason that E4 and E5 are ignored sometimes. In engineering, for convenience, E4 and E5 can be calculated roughly in the way that Q is multiplied by the ratio. This ratio is determined by the type of gun and it can be found in the interior ballistics manual.

3.5. Numerical Calculation of E6

The propellant gas energy, which contains thermal energy of propellant gas and kinetic energy of propellant gas, is intermediate energy and could convert to other kinds of energies. After the projectile leaves the muzzle, the propellant gas will emit out of the gun barrel as exhaust gas. In Equation (14), the other kinds of energies have been obtained, except E6. So, E6 can be given by Equation (23) and its value is shown in Figure 15.
E6 = QQ1E1E2E3E4E5
In the early stage of interior ballistics, since the combustion rate of the propellant is increasing, the values of E1, E2, E3, E4 and E5 are relatively low. So, E6 increases with time. In the later stage of interior ballistics, conversely, because the combustion rate of the propellant continuously decreases until zero, the values of E1, E2, E3, E4 and E5 are relatively high. So, E6 decreases with time.

3.6. Energy Distribution Proportion

According to application characteristics, E1 and E2 belong to the effective energy, which designers aim to get, E3, E4 and E5 belong to the secondary energy, which appears inevitably, and E6 belongs to the potential energy, which is able to convert to other kinds of energies. Based on the calculation results above, Table 4 summarizes energy values at 18.9 ms. Additionally, Figure 16 illustrates the proportionate relationship between each kind of energy and total propellant chemical energy. For instance, the percentage of effective energy, η 1 , can be given as:
η 1 = E 1 + E 2 Q
For this type of fixed gun, its heat efficiency, namely percentage of effective energy, is 31.13%. That is to say, the projectile gets 31.13% of Q as kinetic energy, while the rest (68.87%) of Q converts to other kinds of energies. More specifically, 4.09% of Q dissipates directly and it is a relatively low percentage. E6 accounts for 64.78% of Q and its percentage is relatively high. So, it seems that there is a possibility to decrease the percentage of E6 further. Generally, two methods can be employed to pursue higher heat efficiency. One is recycling the propellant exhaust gas, and the other one is optimizing the parameters of the gun firing system. Considering the difficulty of recycling propellant exhaust gas, an optimization work will be carried out and discussed in the next section.

4. Uncertain Optimization of Energy Conversion

A gun is a firing device, and as a result, the projectile kinetic energy is an important performance indicator. Moreover, propellant gas pressure and rifling stress should also be considered in design since they have a significant effect on the lifetime of the gun barrel. Due to the uncertainty in the interior ballistics stage, an uncertain optimization is performed to obtain higher heat efficiency in this section. This optimization takes projectile kinetic energy as the objective function and takes maximum propellant gas pressure and rifling stress as constraint functions.

4.1. Algorithm Model of Uncertain Optimization

Without loss of generality, the interval method is employed to model the uncertainty, hence the objective function and the constraints with respect to the interval uncertain vector can be written in a form of nonlinear continuous function. By this, an interval optimization problem can be formulated as follows:
min f ( U ) subject   to g i ( U ) ( = , ) b i I = [ b i I , b i R ] ,   i = 1 , 2 , , k , U U I = [ U L , U R ] = U c , U w , U i U i I = [ U i L , U i R ] , i = 1 , 2 , , n ,
where U is a n-dimensional interval uncertainty vector, and the superscripts I , L , R , c and w are the interval, the lower bounds, the upper bounds, the center and the radius of an interval number, respectively. f and g are the interval objective function and the interval constraints, respectively. b i I is the permissible interval of the i-th constraint, and it will degenerate into a real number when b i I =   b i R . Obviously, the traditional optimization algorithms cannot directly solve Equation (25), and some special treatments are required.
The interval order relations can judge that an interval number is better than another, and thus they are usually used to treat the uncertainty objective function. The interval order relation c w , proposed by Ishibuchi et al. [32], is herein adopted to rank the interval numbers. Assuming that there are two interval numbers, A and B, if and only if the center and the radius of B are both smaller than that of A, then interval B is better than interval A, i.e., A c w B (Jiang et al. [33]). Using c w , we can transform the interval objective function in Equation (25) into a multi-objective optimization problem:
min X ( f c ( U ) , f w ( U ) ) f c ( U ) = [ f R ( U ) + f L ( U ) ] / 2 f w ( U ) = [ f R ( U ) f L ( U ) ] / 2
The lower and upper bounds of the interval objective function, i.e., f L ( U ) and f R ( U ) , can be computed through two deterministic optimization processes, which can be formulated as follows:
f L ( U ) = min f ( U ) f R ( U ) = max f ( U ) U U I = [ U L , U R ]
The above equations transform an uncertain objective function into two deterministic ones. Similar to the stochastic programming, f c ( U ) reflects the average performance of the objective function with respect to the uncertainty, and f w ( U ) is the fluctuation range of the objective function caused by the uncertainties. Therefore, we can consider the robustness in this optimization model by minimizing f w ( U ) .
Besides, the interval probability degree model in Reference [34] is adopted to complete the deterministic transformation of uncertain constraints. By this, the inequality constraint, g i ( U ) b i I , can be rewritten as a deterministic one:
P ( g i I ( U ) b i I ) λ i
where λ i is the presupposed threshold value of interval possibility degree. The larger λ i means the more rigorous constraint.
The equality constraint g i I ( U ) = b i I can be transformed as follows:
b i L g i I ( U ) b i R
It can be further transformed through two deterministic inequality constraints:
{ P ( g i I ( U ) b i L ) λ i P ( g i I ( U ) b i R ) λ i
The interval of the constraint, g i I ( U ) , can also be computed through two deterministic optimization processes:
g i I ( U ) = [ g i L ( U ) , g i R ( U ) ] g i L ( U ) = min g i ( U ) g i R ( U ) = max g i ( U )
Through Equations (26)–(31), the interval uncertain optimization problem as in Equation (25) can be rewritten as a multi-objective optimization model:
min [ f c ( U ) , f w ( U ) ] subject   to P ( g i I ( U ) ( = , )   b i I ) λ i , i = 1 , 2 , , k , U U I = [ U L , U R ] where f c ( U ) = [ f R ( U ) + f L ( U ) ] / 2 = ( max f ( U ) + min f ( U ) ) / 2 f w ( U ) = [ f R ( U ) f L ( U ) ] / 2 = ( max f ( U ) min f ( U ) ) / 2 g i I ( U ) = [ g i L ( U ) , g i R ( U ) ] = [ min g i ( U ) , max g i ( U ) ]
Through the above modeling, the obtained Equation (32) is generally solved by the two-layer nested optimization strategy. The inner optimizer is used to compute the interval bounds of the objective function and constraints with respect to the uncertainty, and the outer optimizer is adopted to search for the optimal solution of the design variables. The multi-objective optimization algorithm, namely the non-dominated sorting genetic algorithm-II (NSGA-II), is employed as the outer optimizer since Equation (32) is a typical multi-objective optimization problem. Moreover, the inner optimization is completed by the genetic algorithm (GA), which has fine global search ability so it avoids the local optimal solution. However, it takes more computational costs than the gradient-based algorithms. For more details of this solving strategy, the reader is referred to Reference [34]. Finally, the solving process of the uncertain optimization is depicted in Figure 17.

4.2. Optimization Model and Result Analysis

In design work, the maximum gas pressure that the gun barrel can bear depends on the material of the gun barrel. Once the maximum gas pressure exceeds the yield strength of the gun barrel material, the plastic expansion of the gun barrel would happen. For a large-caliber gun, the lifetime of the gun barrel usually allows several hundred uses due to the wear of the riflings. So, the rifling stress should be in a proper range. By maximizing the interval center of projectile kinetic energy and minimizing its interval radius, and meanwhile considering the constraints of the gas pressure and rifling stress, the optimization model is formulated as follows:
min X [ E A c ( U ) , E A w ( U ) , ζ ] = [ E A c ( U ) , E A w ( U ) , 1 n i = 1 n | U i w A X i w w | ] P ( P b m a x I ( U ) 350 ) λ 1 = 0.9 P ( S r m a x I ( U ) 1000 ) λ 2 = 0.9 U U I = [ U L , U R ]
where E A is the projectile kinetic energy and is equal to the sum of E1 and E2, P b m a x is the max propellant gas pressure at the breech and S r m a x is the max rifling stress at the initial position. Moreover, in actual engineering practice, smaller interval radius of the uncertain design variables usually represents higher manufacturing precision and higher production cost. Therefore, the best design scheme is to achieve the required performances under the maximum parameter errors, which means the low manufacturing costs. Thus, ζ in Equation (33) is a non-negative and dimensionless evaluation coefficient which is designed herein to evaluate the overall uncertainty of all design variables:
ζ = 1 n i = 1 n | U i w A X i w w | ,
where A X i w is the preset maximum possible error of the design variable, and n is the number of design variables. The bigger ζ has the higher overall error level and the lower cost. The initial values of design variables are listed in Table 5.
In Table 5, x 1 is the external diameter of the rotating band, x 2 is the width of the rotating band, x 3 is the rifling width, x 4 is the rifling height, x 5 is the projectile mass, x 6 is the forcing cone taper, x 7 is the propellant thickness and x 8 is the combustion rate coefficient.
Figure 18 depicts optimization results, i.e., pareto front. Here, E A c and E A w denote the center and the radius of the projectile kinetic energy interval, respectively. All solutions in pareto front are feasible and satisfy the constraints. The choice of a final solution depends on the preference of the decision maker, and the maximum projectile kinetic energy is herein taken as the decision preference. Subsequently, Table 6 lists the optimal combination and intervals of design variables. Unlike deterministic optimization, the specific tolerance range of uncertain parameters are obtained, which provides a theoretical reference for the error formulation of design variables. By substituting the data in Table 6 to the original model, Figure 19 shows the specific, optimal ranges of objective function and constraint functions. Detailed optimization results are listed in Table 7.
Figure 19 shows that the performance indicator from the uncertain optimization method is also an interval value. This fully reflects the fluctuation of these dynamic responses which is caused by uncertainty. Hence, the upper and lower bounds of these performance indicators can provide a good reference for designers.
From Table 7, the projectile kinetic energy increases significantly, and the increment is 6.78% with constraint conditions satisfied, which would improve the heat efficiency of the gun firing system. Furthermore, the max rifling stress decreases from 893 to 859.76 Mpa, which conduces to prolong the lifetime of the gun barrel.
Table 8 shows a comparison of energies before and after optimization. As can be seen, after optimization, the percentage of E6 decreased from 64.78% to 61.91%, while the other kinds of energy increased. It means that the optimization improves utilization of E6 and the reduced part of E6 distributes to other kinds of energies. The ideal optimization result is that of increasing E1 and E2, as well as decreasing the other energies. But, it is hard to realize this result. For example, the higher propellant gas caused the increase in not only E1 and E2 but also E4.

5. Conclusions

In this paper, an interior ballistics mathematical model and barrel-projectile finite element model were established and solved. Then, the values of related variables such as pressure and porosity were obtained, the generation process of energies was explained, and especially, the energy versus time data were obtained and discussed. Furthermore, the uncertainty optimization was carried out for higher heat efficiency. From the calculation results, the following conclusions can be drawn:
(1)
For this type of fixed gun, interior ballistics duration was 18.9 ms when using initial parameters. At 18.9 ms, percentage of effective energy was 31.13%, percentage of secondary energy was 4.09% and percentage of potential energy was 64.78%. It indicates that a while small proportion of propellant chemical energy dissipates directly in the interior ballistics stage, most of it converts to propellant gas energy.
(2)
After the uncertain optimization, the percentage of the projectile kinetic energy rose from 31.13% to 33.05%, and the percentage of the improvement was 6.78%. It means that heat efficiency and firing performance improved. In addition, the max rifling stress decreased by 3.79%, which is a benefit for prolonging the lifetime of the gun barrel.
The gun firing model and optimization model built in this paper could calculate and optimize not only energies but also other dynamical variables such as stress and strain. They would be a benefit for the improvement of the comprehensive performance of the gun firing system.

Author Contributions

T.X. completed the calculations and the framework of paper. G.Y. reviewed the paper and validated the calculation results. L.W. supplied the optimization model. Q.S. supplied parameter values and experiment data. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the National Natural Science Foundation of China (Grant number: 11572158 and 51705253).

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Liao, X.; Zhao, J.; Xu, B.; Wang, Z.S. Modeling and simulation of interior ballistics based on actual combustion characteristics of propellants. J. Energetic Mater. 2015, 33, 155–166. [Google Scholar] [CrossRef]
  2. Monreal-González, G.; Otón-Martínez, R.A.; Velasco, F.J.S.; García-Cascales, J.R.; Moratilla-Fernández, D. A one-dimensional code for the analysis of internal ballistic problems. In Proceedings of the 28th International Symposium on Ballistics, Atlanta, GA, USA, 22–26 September 2014. [Google Scholar]
  3. Monreal-González, G.; Otón-Martínez, R.A.; Velasco, F.J.S.; García-Cascáles, J.R.; Ramírez-Fernández, F.J. One-Dimensional Modelling of Internal Ballistics. J. Energetic Mater. 2017, 35, 397–420. [Google Scholar] [CrossRef]
  4. Otón-Martínez, R.A.; Monreal-González, G.; García-Cascales, J.R.; Vera-García, F.; Velasco, F.J.S.; Ramírez-Fernández, F.J. An approach formulated in terms of conserved variables for the characterisation of propellant combustion in internal ballistics. Int. J. Numer. Methods Fluids 2015, 79, 394–415. [Google Scholar] [CrossRef]
  5. Jaramaz, S.; Mickovic, D.; Elek, P. Two-phase flows in gun barrel: Theoretical and experimental studies. Int. J. Multiph. Flow 2011, 37, 475–487. [Google Scholar] [CrossRef]
  6. López, M.; García, C.; Velasco; Otón, M. An Energetic Model for Detonation of Granulated Solid Propellants. Energies 2019, 12, 4459. [Google Scholar] [CrossRef] [Green Version]
  7. Wu, B.; Zheng, J.; Tian, Q.-T.; Zou, Z.-Q.; Chen, X.-L.; Zhang, K.-S. Friction and wear between rotating band and gun barrel during engraving process. Wear 2014, 318, 106–113. [Google Scholar] [CrossRef]
  8. Wu, B.; Fang, L.H.; Zheng, J.; Yu, X.H.; Jiang, K.; Wang, T.; Hu, L.M.; Chen, Y.L. Strain Hardening and Strain-Rate Effect in Friction Between Projectile and Barrel During Engraving Process. Tribol. Lett. 2019, 67. [Google Scholar] [CrossRef]
  9. Wu, B.; Zheng, J.; Tian, Q.-T.; Zou, Z.-Q.; Yu, X.-H.; Zhang, K.-S. Tribology of rotating band and gun barrel during engraving process under quasi-static and dynamic loading. Friction 2014, 2, 330–342. [Google Scholar] [CrossRef] [Green Version]
  10. Wu, B.; Fang, L.-H.; Chen, X.-L.; Zou, Z.-Q.; Yu, X.-H.; Chen, G. Fabricating Aluminum Bronze Rotating Band for Large-Caliber Projectiles by High Velocity Arc Spraying. J. Therm. Spray Technol. 2013, 23, 447–455. [Google Scholar] [CrossRef]
  11. Wu, B.; Zheng, J.; Qiu, J.; Zou, Z.-Q.; Yan, K.-B.; Chen, X.-L.; Hu, L.-M.; Zhang, K.-S.; Chen, R.-G. Preparation of the projectile rotating band and its performance evaluation. J. Adhes. Sci. Technol. 2016, 30, 1143–1164. [Google Scholar] [CrossRef]
  12. Edward Alexander, J. Advanced Gun System Gun and Projectile Dynamic Model Results and Correlation to Test Data. J. Adhes. Sci. Technol. 2012, 134. [Google Scholar] [CrossRef]
  13. Sudarsan, N.V.; Sarsar, R.B.; Das, S.K.; Naik, S.D. Prediction of Shot Start Pressure for Rifled Gun System. Def. Sci. J. 2018, 68, 144–149. [Google Scholar] [CrossRef]
  14. Dursun, T. Effects of Projectile and Gun Parameters on the Dispersion. Def. Sci. J. 2020, 70, 166–174. [Google Scholar] [CrossRef]
  15. Carlucci, D.E.; Frydman, A.M.; Cordes, J.A. Mathematical Description of Projectile Shot Exit Dynamics (Set-Forward). J. Appl. Mech. 2013, 80, 979–985. [Google Scholar] [CrossRef]
  16. Ding, C.; Liu, N.; Zhang, X. A mesh generation method for worn gun barrel and its application in projectile-barrel interaction analysis. Finite Elem. Anal. Des. 2017, 124, 22–32. [Google Scholar] [CrossRef]
  17. Zhang, L.; Chen, Z.; Wen, D.; Wang, X.; Zhang, D.; Liang, J. Estimation of the Time-Varying High-Intensity Heat Flux for a Two-Layer Hollow Cylinder. Energies 2018, 11, 3332. [Google Scholar] [CrossRef] [Green Version]
  18. Hill, R.D.; Conner, J.M. Transient Heat Transfer Model of Machine Gun Barrels. Mater. Manuf. Process. 2012, 27, 840–845. [Google Scholar] [CrossRef]
  19. Işık, H.; Göktaş, F. Cook-off analysis of a propellant in a 7.62 mm barrel by experimental and numerical methods. Appl. Therm. Eng. 2017, 112, 484–496. [Google Scholar] [CrossRef]
  20. Değirmenci, E.; Evci, C.; Işık, H.; Macar, M.; Yılmaz, N.; Dirikolu, M.H.; Çelik, V. Thermo-mechanical analysis of double base propellant combustion in a barrel. Appl. Therm. Eng. 2016, 102, 1287–1299. [Google Scholar] [CrossRef]
  21. Talaee, M.R.; Hosseinli, S.A. Theoretical simulation of temperature distribution in a gun barrel based on the DPL model. J. Theor. Appl. Mech. 2019, 57, 685–696. [Google Scholar] [CrossRef]
  22. Sefiane, K.; Benielli, D.; Steinchen, A. A new mechanism for pool boiling crisis, recoil instability and contact angle influence. Colloids Surf. A Physicochem. Eng. Asp. 1998, 142, 361–373. [Google Scholar] [CrossRef]
  23. Kathe, J.K. Design on the recoil control orifice for high energy recoil force. In Proceedings of the ADD 11th Land Vehicle Conference, Daejon, Korea, 7–11 July 2003. [Google Scholar]
  24. Kathe, J.K. Analysis on the forward motion for the soft recoil mechanism. In Proceedings of the Agency for Defense Development (ADD) 40th Foundation Conference, Seoul, Korea, 23–28 August 2010. [Google Scholar]
  25. Mitchell, M.R.; Link, R.E.; Wu, Y.C.; Chang, H.; Tsung, T.T. Dynamic Characteristics of a Recoil System When Firing Projectiles with Mach 4.4 Muzzle Velocity from a 105 mm Cannon. J. Test. Eval. 2011, 39, 704–714. [Google Scholar] [CrossRef]
  26. Hassaan, G.A. Dynamics of a cannon barrel-recoil mechanism with a nonlinear hydraulic damper. Int. J. Sci. Eng. Technol. 2014, 5, 82–91. [Google Scholar]
  27. Bartnik, R.; Buryn, Z.; Hnydiuk-Stefan, A.; Skomudek, W.; Otawa, A. Thermodynamic and Economic Analysis of Trigeneration System Comprising a Hierarchical Gas-Gas Engine for Production of Electricity, Heat and Cold. Energies. 2020, 13, 1006. [Google Scholar] [CrossRef] [Green Version]
  28. Jang, J.; Choi, J.; Yi, H.; Park, S. Effects of the Bore to Stroke Ratio on Combustion, Gaseous and Particulate Emissions in a Small Port Fuel Injection Engine Fueled with Ethanol Blended Gasoline. Energies 2020, 13, 321. [Google Scholar] [CrossRef] [Green Version]
  29. Nonthakarn, P.; Ekpanyapong, M.; Nontakaew, U.; Bohez, E. Design and Optimization of an Integrated Turbo-Generator and Thermoelectric Generator for Vehicle Exhaust Electrical Energy Recovery. Energies 2019, 12, 3134. [Google Scholar] [CrossRef] [Green Version]
  30. Pagán Rubio, J.A.; Vera-García, F.; Hernandez Grau, J.; Muñoz Cámara, J.; Albaladejo Hernandez, D. Marine diesel engine failure simulator based on thermodynamic model. Appl. Therm. Eng. 2018, 144, 982–995. [Google Scholar] [CrossRef]
  31. Vera-García, F.; Pagán Rubio, J.A.; Hernández Grau, J.; Albaladejo Hernández, D. Improvements of a Failure Database for Marine Diesel Engines Using the RCM and Simulations. Energies 2019, 13, 104. [Google Scholar] [CrossRef] [Green Version]
  32. Ishibuchi, H.; Tanaka, H. Multiobjective programming in optimization of the interval objective function. Eur. J. Oper. Res. 1990, 48, 219–225. [Google Scholar] [CrossRef]
  33. Jiang, C.; Han, X.; Liu, G.P. A sequential nonlinear interval number programming method for uncertain structures. Comput. Methods Appl. Mech. Eng. 2008, 197, 4250–4265. [Google Scholar] [CrossRef]
  34. Jiang, C.; Han, X.; Liu, G.R.; Liu, G.P. A nonlinear interval number programming method for uncertain optimization problems. Eur. J. Oper. Res. 2008, 188, 1–13. [Google Scholar] [CrossRef]
Figure 1. Schematic drawing of propellant location.
Figure 1. Schematic drawing of propellant location.
Energies 13 05824 g001
Figure 2. Barrel-projectile finite element model.
Figure 2. Barrel-projectile finite element model.
Energies 13 05824 g002
Figure 3. Flow chart of energies calculation.
Figure 3. Flow chart of energies calculation.
Energies 13 05824 g003
Figure 4. Schematic representation of the energy conversion process.
Figure 4. Schematic representation of the energy conversion process.
Energies 13 05824 g004
Figure 5. (a) Distribution of propellant gas pressure; (b) projectile bottom pressure.
Figure 5. (a) Distribution of propellant gas pressure; (b) projectile bottom pressure.
Energies 13 05824 g005
Figure 6. (a) Translational velocity and translational kinetic energy of the projectile, (b) rotational velocity and rotational kinetic energy of the projectile.
Figure 6. (a) Translational velocity and translational kinetic energy of the projectile, (b) rotational velocity and rotational kinetic energy of the projectile.
Energies 13 05824 g006
Figure 7. Representation of size relationship of rotating band diameter and barrel bore diameter.
Figure 7. Representation of size relationship of rotating band diameter and barrel bore diameter.
Energies 13 05824 g007
Figure 8. Mises stress and morphology at 6.2 ms in the engraving period.
Figure 8. Mises stress and morphology at 6.2 ms in the engraving period.
Energies 13 05824 g008
Figure 9. Mises stress and morphology at 9.5 ms in the translation-rotation period.
Figure 9. Mises stress and morphology at 9.5 ms in the translation-rotation period.
Energies 13 05824 g009
Figure 10. The sum of frictional energy and plastic deformation energy.
Figure 10. The sum of frictional energy and plastic deformation energy.
Energies 13 05824 g010
Figure 11. (a) Porosity distribution, (b) percentage of uncombusted propellant energy.
Figure 11. (a) Porosity distribution, (b) percentage of uncombusted propellant energy.
Energies 13 05824 g011
Figure 12. (a) Pressure at breech and gun barrel displacement, and (b) recoil energy.
Figure 12. (a) Pressure at breech and gun barrel displacement, and (b) recoil energy.
Energies 13 05824 g012
Figure 13. (a) Average temperature of propellant gas, (b) thermal energy of the gun barrel.
Figure 13. (a) Average temperature of propellant gas, (b) thermal energy of the gun barrel.
Energies 13 05824 g013
Figure 14. Temperature distribution of inner surface of the gun barrel at 10.5, 14.8 and 18.9 ms.
Figure 14. Temperature distribution of inner surface of the gun barrel at 10.5, 14.8 and 18.9 ms.
Energies 13 05824 g014
Figure 15. Propellant gas energy.
Figure 15. Propellant gas energy.
Energies 13 05824 g015
Figure 16. Sector diagram of energy distribution proportion.
Figure 16. Sector diagram of energy distribution proportion.
Energies 13 05824 g016
Figure 17. The solving process of the uncertain optimization.
Figure 17. The solving process of the uncertain optimization.
Energies 13 05824 g017
Figure 18. Pareto front.
Figure 18. Pareto front.
Energies 13 05824 g018
Figure 19. (a) Optimized value and intervals of propellant gas pressure at the breech, (b) optimized value and intervals of rifling stress and (c) optimized value and intervals of projectile kinetic energy.
Figure 19. (a) Optimized value and intervals of propellant gas pressure at the breech, (b) optimized value and intervals of rifling stress and (c) optimized value and intervals of projectile kinetic energy.
Energies 13 05824 g019
Table 1. Parameter values of propellant and mechanical structure.
Table 1. Parameter values of propellant and mechanical structure.
Parameter NameParameter ValueParameter NameParameter Value
gun barrel length8.06 mcaliber155 mm
propellant mass15.6 kgprojectile mass45.5 kg
propellant length1.08 mprojectile rotational inertia0.11 kg∙m2
heat value coefficient3.87 × 106 J∙kg−1
Table 2. Material parameters of the barrel-projectile finite element model.
Table 2. Material parameters of the barrel-projectile finite element model.
PartMaterialYoung’s ModulusPoisson’s RatioDensitySpecific HeatThermal Conductivity
projectilesteel190 GPa0.35420 kg ·   m   3 460 J ·   kg 1 ·   K 1 40 W ·   m 1 ·   K 1
gun barrelsteel210 GPa0.37830 kg ·   m   3 473 J ·   kg 1 ·   K 1 30 W ·   m 1 ·   K 1
riflingssteel210 GPa0.37830 kg ·   m   3 473 J ·   kg 1 ·   K 1 30 W ·   m 1 ·   K 1
rotating bandcopper124 GPa0.348960 kg ·   m   3 420 J ·   kg 1 ·   K 1 380 W ·   m 1 ·   K 1
Table 3. Projectile muzzle translational velocity measured by firing tests.
Table 3. Projectile muzzle translational velocity measured by firing tests.
123456Average
919 m∙s−1913 m∙s−1925 m∙s−1933 m∙s−1919 m∙s−1917 m∙s−1921 m∙s−1
Table 4. Energy values at 18.9 ms.
Table 4. Energy values at 18.9 ms.
t/msQ/JQ1 /JE1 /JE2 /JE3 /JE4/JE5/JE6/J
18.96.05 × 10701.87 × 1071.34 × 1051.42 × 1063.39 × 1057.19 × 1053.92 × 107
Table 5. The nominal values and the searching ranges of the design variables.
Table 5. The nominal values and the searching ranges of the design variables.
Lower BoundInitiationUpper Bound
x 1 / ( mm ) 79.880.081.0
x 2 / ( mm ) 5458.5262.0
x 3 / ( mm ) 5.906.786.9
x 4 / ( mm ) 2.102.953.10
x 5 / ( kg ) 45.245.5046.0
x 6 / ( ° ) 1.782.055.71
x 7 / ( mm ) 2.102.302.40
x 8 / ( m / ( s . Pa n ) ) 1.57 × 10−81.67 × 10−81.93 × 10−8
Table 6. The optimized intervals of design variables.
Table 6. The optimized intervals of design variables.
Lower BoundOptimizationUpper Bound
x 1 / ( mm ) 80.1280.2680.49
x 2 / ( mm ) 54.9555.2156.88
x 3 / ( mm ) 6.126.316.55
x 4 / ( mm ) 3.013.063.09
x 5 / ( kg ) 45.2845.3145.47
x 6 / ( ° ) 5.215.405.66
x 7 / ( mm ) 2.212.282.33
x 8 / ( m / ( s . Pa n ) ) 1.69 × 10−81.72 × 10−81.78 × 10−8
Table 7. Detailed optimization results and comparison of objective function and constraints.
Table 7. Detailed optimization results and comparison of objective function and constraints.
InitiationLower BoundOptimizationUpper Bound
E A 1.881.932.012.08
P b m a x 337.61344.54347.48349.91
S r m a x 893.68839.38859.76872.31
Table 8. Comparison of energies before and after optimization.
Table 8. Comparison of energies before and after optimization.
QE1E2E3E4E5E6
initial energies6.05 × 107 J1.87 × 107 J1.34 × 105 J1.42 × 106 J3.39 × 105 J7.19 × 105 J3.92 × 107 J
initial percent100%30.91%0.22%2.35%0.56%1.18%64.78%
optimized energies6.05 × 107 J1.98 × 107 J2.11 × 105 J1.61 × 106 J4.63 × 105 J9.62 × 105 J3.74 × 107 J
optimized percent100%32.71%0.34%2.66%0.77%1.56%61.91%
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Xin, T.; Yang, G.; Wang, L.; Sun, Q. Numerical Calculation and Uncertain Optimization of Energy Conversion in Interior Ballistics Stage. Energies 2020, 13, 5824. https://doi.org/10.3390/en13215824

AMA Style

Xin T, Yang G, Wang L, Sun Q. Numerical Calculation and Uncertain Optimization of Energy Conversion in Interior Ballistics Stage. Energies. 2020; 13(21):5824. https://doi.org/10.3390/en13215824

Chicago/Turabian Style

Xin, Tong, Guolai Yang, Liqun Wang, and Quanzhao Sun. 2020. "Numerical Calculation and Uncertain Optimization of Energy Conversion in Interior Ballistics Stage" Energies 13, no. 21: 5824. https://doi.org/10.3390/en13215824

APA Style

Xin, T., Yang, G., Wang, L., & Sun, Q. (2020). Numerical Calculation and Uncertain Optimization of Energy Conversion in Interior Ballistics Stage. Energies, 13(21), 5824. https://doi.org/10.3390/en13215824

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop