Stress Analysis of Gun Barrel Subjected to Dynamic Pressure
H. Babaei^{*}, M. Malakzadeh, H. Asgari
Department of Mechanical Engineering, Engineering Faculty, University of Guilan, Rasht, Iran
Email address:
To cite this article:
H. Babaei, M. Malakzadeh, H. Asgari. Stress Analysis of Gun Barrel Subjected to Dynamic Pressure. International Journal of Mechanical Engineering and Applications. Vol. 3, No. 4, 2015, pp. 71-80. doi: 10.11648/j.ijmea.20150304.14
Abstract: In the optimal design of a modern gun barrel, there are some aspects to be considered. One of the main factor is internal ballistic which consist of pressure-time, pressure-distance, velocity-time and distance-time curves. In this paper, a simple analytical solution for the plastic stress of an internally pressurized open-ended thick-walled cylinder made of hardening steel which is the closest model to gun barrel is obtained in perfectly plastic and plane stress condition by using energy method and the yield criterion of Von Mises and adding rifle grooves and choosing stress components as basic unknowns and ballistic pressure equation as known. Then results of analytical solution are compared to a numerical model and verified a very well and reliable accuracy. So the resultant can be used easily in calculation of radial expansion velocity and compressive pressure.
Keywords: Gun Barrel, Stress, Dynamic Pressure
1. Introduction
Many papers are published around the matter of gun design and especially about gun tube design from the beginning of improved guns design generation until now, but reaching to variation sources around this subject is impossible because of preventing of publishing of this technology and security and military problems of some countries. Therefore many engineers and designers who work on this field assume a reasonable simplification that considers gun barrel as an open-ended thick-walled cylinder in which an explosive causes stress and deformation on the wall by creation an immediate pressure. According to this assumption investigation must be done under plane stress condition (. Few satisfactory theoretical solution based on Von Mises yield criterion that are reliable and convenient for engineering use have been obtained for an open-ended thick-walled cylinder. A closed-form solution for the stress components was given by Nadai using an auxiliary-variable method and the deformation theory of Hencky [1]. A set of analytical expressions for the elasto-plastic stress and displacement components were obtained by Davidson et al. on the basis of an empirical relationship resulting from tests [2]. Two mathematically consistent analytical solutions to the strains and displacements were obtained according to the deformation theory of Hencky and to the flow theory of Prandtl-Reuss, respectively by P.C.T. Chen using a modified Nadai's auxiliary variable method [3]. D.R. Bland et al. derived equations of stress distribution in thick-walled cylinder under internal and external pressure using numerical solution in 1956 [4]. They studied terms of equivalent stress and strain for the two conditions of open-ended () and close-ended () for this analysis but the way of calculating of stresses in the open-ended condition which is the appropriate model for gun barrel had not been mentioned. V. A. Adintsov et al. achieved an equation for radial and hoop stresses and radial expansion velocity by analytical investigating on perfectly plastic thick-walled cylinder behavior under internal explosive pressure using energy method [5]. One of the weak points of this study is analyzing the material as a perfectly plastic one but not considering the strain rate. An analytical pattern and some numerical results for an internal-pressurized open-ended thick-walled cylinder made of linear-hardening material was given by Xu Hong and Chen Shuning on the basis of the Prandtl-Reuss flow theory [6]. A theoretical and numerical analysis for a similar problem was performed by LiuYong using the elasto-plastic mixed-boundary-element method [7]. Gao Xin-lin et al. obtained a closed form analytical solution for stress, strain and displacement components for an internal-pressurized elasto-plastic open-ended thick-walled cylinder in 1991 [8]. In this paper it is shown that the solution is a general one with on the one hand Nadai's known solution for stress components and on the other P.C.T. Chen's solution for strain components of an open-ended thick-walled cylinder made of elastic-perfectly-plastic material as its two specific cases. In this analysis it is assumed an auxiliary function for radial stress in terms of current radius after performing pressure for obtaining stress components at plane stress condition in plastic zone which is very complicated and Time-consuming. X.-L Gao et al. exhibited a technique for elastic-plastic internally pressurized cylinder analyzing in that both Tresca and Von Mises criterion is used for three conditions of plane stress, plane strain and close-ended in order to calculate the stresses and strains in 2003 [9]. Also the deformation energy of cylinder wall is obtained after the pressure is performed. Lack of considering the term of wall expansion acceleration is one of it’s weaknesses. Li Mao-lin et al. investigated plastic limit load of viscoplastic thick-walled cylinder and spherical shell subjected to internal pressure analytically using a strain gradient plasticity theory in the paper which is presented in 2008 [10]. Results show that the size eﬀect is more evident with increasing strain or strain rate sensitivity index, but the weak point of this approach is that the viscoplastic model analysis doesn't include material behavior in all part of the wall. Bagheri et al. presented a paper with the purpose of deriving a mathematical model for expansion isotropic thick-walled aluminum cylinder containing TNT in which, JWL equation of state is considered for explosive products [11]. As a result the equations of radial and hoop stress and radial expansion velocity is obtained. However, it must be noted that, each of these known analytical solutions is either over-simplified in the material model but time-consuming in calculation [12,13] or too complicated in the expressions proposed but incomplete in content [14,15].
In this analysis after-explosion phenomena in a thick walled steel cylinder, such as radial pressure on the wall of cylinder, the radial drift, the rate of expansion of the radial drift and radial and circumferential stresses, are considered. This study starts with a simple analytical model using energy-based methods for mentioned parameters. Subsequently, by choosing a cylinder which has specific material and specific geometric dimensions as well as an explosive substance, the mentioned analysis is carried out. One of the advantages of this study is simplicity in calculating the rate of the radial expansion and hence, changes in the radius of the cylinder and its stresses and also the acceleration of the wall are measurable.
2. Theoretical Analyses
2.1. Basic Equations of Thick-Walled Cylinders
Considering Figure 1, the following equations can be obtained [16]
(1)
Where and are radial and hoop stresses respectively. The equation above shows the equation of the thick-walled cylinder in static state [17] and considering the open-ended condition (), and in a dynamic state it turns to the following equation considering the movement of the element shown in figure 1. It should be mentioned that the term related to density in acceleration is added to the equation [18].
(2)
Where and are radial expansion velocity and acceleration respectively, and is the density of cylinder material.
2.2. Conservation Law
Conservation law in physics expresses that the whole energy of an isolated system is consistent with time [19]. Energy neither produces nor wastes, but transforms from a state to another. In this study the chemical energy of explosives turns to kinetic energy. this law is explained in the form of the equation below.
where, and represent initial internal energy and current energy of the substances of explosion, respectively. is the kinetic energy of the wall of the cylinder, is the work of plastic deformation and is the work of friction force of the rifle. If we represent the equation above per unit of length, will have the following equation.
(3)
The terms of the equation above are determined during the trend of analysis in this study.
2.3. Adiabatic Expansion Law
In this research the expansion of the wall due to explosion is considered adiabatic and obeys the equation below [20].
(4)
Where p and V is pressure and special volume respectively and 𝛾 is adiabatic expansion coefficient which is specified for any material.
2.4. Gunnery Internal Ballistic Equations
Internal ballistic equations consist of formulas and diagrams related to time-pressure, distance-pressure and distance-time. It should be mentioned that, because the selected weapon in this study is M24 which has a 500 mm barrel and the firing time is 1 ms, we consider the pressure curve yielded from Russel Model within these ranges and subsequently, we curve fit this diagram in software. The results are as follows.
The equation of pressure based on time:
(5)
where P represents internal ballistic pressure and t represents the time interval between triggering and firing bullet through the barrel.
, , ,
, ,
,
Hence, figure 2 shows the diagram of equation (5)
The equation of distance (the amount of space passed by the bullet into the barrel) based on time is presented as follows.
(6)
In which x is the length of the barrel
, , ,
The yielded diagram from the equation above is demonstrated in the Figure 3.
As a result, the diagram of pressure versus distance is demonstrated as Figure 4.
3. Research Method
3.1. Assumptions
• In this analysis the cylinder containing explosives is assumed perfectly plastic, homogeny and incompressible.
• The used criterion is Von Mises. As a result of recent studies and comparing with the analysis a more appropriate criterion is introduced which is Von Mises criterion. Due to utilizing plastic work, the simulation and the analysis are independent of stiffness strains and thermal variations.
• The volume of the barrel is assumed constant before and after explosion and there is no body force affecting the barrel but the pressure of explosion.
• In this analysis the shock wave loading and swinging effects on the diameter of the barrel are ignored.
• The barrel analysis is carried out under the conditions of plane stresses ().
• It is assumed that axial stresses of the barrel caused by the friction force of the projectile do not affect the analysis results.
• In the equation of conservation law, the term related to the energy of explosion is neglected.
3.2. Analysis
The stresses and are principal. The internal and external initial cylinder radii before the deformation due to explosion are denoted by and respectively and radii after deformation by a and b, and r is the current radii. We assume the deformation of the tube wall is perfectly plastic, so there is no changing in volume, therefore from equality of current and initial volume we have:
(7)
By substituting b in (5) and differentiating, (internal radius expansion velocity) and (current radius expansion velocity) is obtained, so [21]:
(8)
By differentiating form (6) and substituting in (2) we have
(9)
The yield criterion adopted here is the Von Mises criterion according to recent studies which say it is more accurate than others like Tresca, therefore
(10)
Where for Von Mises criterion and is yield stress. By integrating (7) and applying the boundary condition in :
(11)
Where is the pressure dependent to distance and time that is considered as p from now on. We have in , the equation of internal radii acceleration is as follow:
(12)
Now we want to calculate the amount of and . Accordingly, in order to preventing of using complicated calculation of differential equation, we present a simple solution using energy method (conservation energy law) and thermodynamic fundamental. Hence we calculate each term of equation (3).
3.2.1. Calculation of Detonation Products Internal Energy Per Unit of Length ()
is instantaneous pressure and obtained from
(13)
is the density of the explosive and D is detonation velocity which are specified for any explosive and obtained empirically and accessible in engineering handbooks.
Internal energy for an ideal gas is in form in which is volume per unit of length, so we have
(14)
3.2.2. Calculation of Wall Kinetic Energy Per Unit of Length ()
The equation of wall kinetic energy is in form , substituting (8) and in the equation, it turns to , in which is barrel material density, so
(15)
3.2.3. Calculation of Work of Plastic Deformation Per Unit of Length ()
Work of plastic deformation equation is , where is equivalent strain, We assume , for solving the integral we need to obtain an equation for the equivalent strain in plane stress condition [4].
(16)
Where , in which is the boundary of elastic and plastic domain and is poisson coefficient. Since the workplace of this research is perfectly plastic, we have , as a result the equation of equivalent is
(17)
In order to obtain the amount of we differentiate from (8), therefore
(18)
Substituting (12) and (18) into (17)
(19)
As a result
Finally work of plastic deformation equation is in form
(20)
3.2.4. Calculation of Work of Rifiling Friction Force Per Unit of Length ()
According to number of grooves the equation can be in form , in which is gas pressure in the place x and time t, is gas force on the chamber basis in the place x and time t, D is the caliber, n is the number of grooves, h is the grooves depth, and is land width.
On the other hand work of friction force is
Where is average friction coefficient, is of physical properties of bullet that is 0.74 on average, α is groove angle and equal to 8 degree according to compare variation of guns. So . we choose the barrel length 500 mm. finally rifle work of friction force per unit length is
(21)
Substituting equations
(22)
For simplification of recent equation, we assume , , , , , , and
So is obtained as
(23)
Now by substituting into (12), (11) and (10) we obtain , and respectively.
3.3. Numerical Simulation
In this research ABAQUS 6.12 and dynamic explicit solver is used to simulate and solve the problem.
CATIA software is used to model the M24 barrel. The barrel have four grooves with the pitch of a round per twelve inches in accordance with Figure 5.
This model imported to ABAQUS to numerical analysis. Other assumptions and geometrical and physical properties of material and explosive properties is shown in table 1 and table2 respectively.
Dimension and size | |
internal radius (mm) | 3.81 |
External radius (mm) | 14.6 |
Final expansion radius (mm) | 14 |
Barrel length (mm) | 600 |
Groove width (mm) | 3 |
Groove depth (mm) | 2 |
Number of rifles | 4 |
Physical properties | |
Density (kg/) | 7800 |
Yield stress (Mpa) | 330 |
Poisson coefficient | 0.3 |
Density (kg/) | 1717 |
Detonation velocity (m/s) | 7980 |
Gas coefficient 𝛾 | 2.7 |
The equation of pressure-distance imported in analytical field section to investigate pressure effect. The initial part of barrel is bounded to create boundary condition.
According to the analytical solution, period of 0.001 sec is chosen to numerical analysis in step section in ABAQUS. Mesh section is used for modeling and reticulation of problem geometry. The model have free reticulation and the elements are quadratic. Form of meshing used in barrel is shown as Figure 6.
Solving the problem is done by job section. Other setting related to problem solution in parallel processing condition, output numbers accuracy, etc can be done in this section. Stress contour among barrel length is shown in Figure 7.
Graph of stress distribution on grooves among of barrel length is presented in Figure 8.
4. Conclusion
Diagram of radial expansion velocity variation per time is shown in Figureure 9 as one of the important results of this research. Maximum tube expansion velocity is 722.04 m/s in 0.39 sec for 4 rifiled barrel and 709.02 m/s in same time for 3 rifiled barrel on internal surface. The tube has maximum deformation and minimum thickness in this time. So the more number of grooves, the more expansion velocity.
Figureure 10 shows comparison of stress diagram per barrel length between two condition of analytical and numerical solution on internal surface. Maximum compressive pressure is about 455 Mpa due to analytical solution and about 465 Mpa due to numerical one. As a result the equations obtained from theoretical analysis are in well accordance on numerical solution for calculating the amounts of compressive pressure in order to acceptable difference between the two graphs.
• In this paper, equilibrium equations of cylinder is considered as the closest model to gun barrel and a simple equation for radial expansion velocity is obtained considering boundary conditions. Expansion velocity can be calculated in every times entering boundary condition. The method used in this research is unique according to comparison analytical and numerical output, therefore it can be used in real situation, lab environment and testing experimental related to gunnery in order to reducing costs and time-consuming.
• Basically choosing the sort of equation of state governing the explosive after explosion and material model is not affecting on the analytical technique but it is necessary to be chosen in such a way to conform to the terms and actual results. Hence, internal ballistic pressure equation is considered in general form, but it should be mentioned the equation has been extracted with very high precision software. But from the convergence of analysis and simulation results can be deduced the equation of state describes explosive products expansion behavior very well with the related coefficients. So the resultant can be used in calculation easily.
• Radial velocity in the analysis is verified according to comparison with simulation therefore other extractive parameters accuracy like radial and hoop stresses depended to radial velocity is sufficient.
• Generally it can be stated that analysis technique introduced in this paper is accurate and the application of material models in the process of expansion recall is complete according to the simulation results and a small amount of errors. So this model can be used to predict the results before testing experimental.
References