CN114639450A - Method for calculating ballistic characteristic parameters in mixed charge under plasma ignition - Google Patents

Method for calculating ballistic characteristic parameters in mixed charge under plasma ignition Download PDF

Info

Publication number
CN114639450A
CN114639450A CN202210216931.XA CN202210216931A CN114639450A CN 114639450 A CN114639450 A CN 114639450A CN 202210216931 A CN202210216931 A CN 202210216931A CN 114639450 A CN114639450 A CN 114639450A
Authority
CN
China
Prior art keywords
gunpowder
projectile
powder
equation
ballistic
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Pending
Application number
CN202210216931.XA
Other languages
Chinese (zh)
Inventor
李海元
张皓婷
栗保明
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Nanjing University of Science and Technology
Original Assignee
Nanjing University of Science and Technology
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Nanjing University of Science and Technology filed Critical Nanjing University of Science and Technology
Priority to CN202210216931.XA priority Critical patent/CN114639450A/en
Publication of CN114639450A publication Critical patent/CN114639450A/en
Pending legal-status Critical Current

Links

Images

Classifications

    • G—PHYSICS
    • G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16C—COMPUTATIONAL CHEMISTRY; CHEMOINFORMATICS; COMPUTATIONAL MATERIALS SCIENCE
    • G16C20/00—Chemoinformatics, i.e. ICT specially adapted for the handling of physicochemical or structural data of chemical particles, elements, compounds or mixtures
    • G16C20/10—Analysis or design of chemical reactions, syntheses or processes

Landscapes

  • Chemical & Material Sciences (AREA)
  • Analytical Chemistry (AREA)
  • Chemical Kinetics & Catalysis (AREA)
  • Crystallography & Structural Chemistry (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Bioinformatics & Cheminformatics (AREA)
  • Bioinformatics & Computational Biology (AREA)
  • Computing Systems (AREA)
  • Theoretical Computer Science (AREA)
  • Plasma Technology (AREA)

Abstract

The invention provides a method for calculating ballistic characteristic parameters in mixed charge under plasma ignition. Determining a pulse power supply discharge parameter; establishing a gunpowder burning rate formula under the plasma enhancement effect; determining shape functions of the two gunpowders; establishing a projectile motion equation and an electrothermal chemical emission internal trajectory equation; and substituting the input electric energy into an electrothermal chemical launching internal trajectory equation, combining a gunpowder burning rate formula and a shape function to carry out time propulsion solution, and finally obtaining the in-chamber pressure and the projectile velocity curve in the projectile discharging process by using a projectile motion equation and the electrothermal chemical launching internal trajectory equation. The method solves the problem that only a single charge structure is considered in the existing electrothermal chemical launch inner trajectory model, and can quickly and effectively calculate the in-bore pressure and the projectile velocity curve in the process of discharging the electrothermal chemical launch projectile out of the bore under the mixed charge structure.

Description

Method for calculating ballistic characteristic parameters in mixed charge under plasma ignition
Technical Field
The invention relates to an electrothermal chemical emission technology, in particular to a calculation method of internal ballistic characteristic parameters under plasma ignition considering a mixed charge structure.
Background
The electrothermal chemical launching technology is an energy mixed launching technology, and can provide partial energy by energetic materials and also can complete the propelling and launching of the projectile by the action of electric energy. A plasma generator is used in the launching device to replace a conventional device for ignition, a pulse forming network is utilized to output strong pulse current to electrodes at two ends of a capillary, working medium in the capillary is heated by electrode discharge, the working medium is ionized to generate high-temperature and high-pressure plasma and is injected into a combustion chamber at high speed, the high-temperature and high-pressure plasma and gas generated by gunpowder combustion are pushed into a bullet together to advance in an accelerating manner, the bullet obtains high initial speed, the muzzle kinetic energy is improved by 15% -35% compared with that of a traditional artillery, and low-cost engineering transformation can be realized on the traditional artillery.
An existing electrothermal chemical emission inner ballistic calculation model is mainly established based on a single charge mode and mainly simulates the influence of the change of the discharge condition of a pulse network on the inner ballistic performance. The document 'plasma enhancement numerical simulation of an electrothermal chemical gun' combines the characteristics of pulse formation network discharge, adds plasmas into an inner trajectory basic equation in the form of source items, establishes an engineering practical inner trajectory calculation model of the electrothermal chemical gun, simulates the influence of discharge voltage and capacitance change on the inner trajectory performance, but does not research the influence caused by charge structure change.
The charging structure is a main means for improving weapon power, increasing the initial speed of the projectile and inhibiting harmful phenomena of shooting when the overall structure of the artillery is fixed. By adjusting the charging structure, the contradiction among various performances of the weapon can be adjusted, and the requirement of comprehensive performance of the trajectory can be better met. In the literature, cross-cut rod-shaped and coated granular propellant mixed charging constant-volume combustion performance, 9/19 cross-cut rod-shaped propellant and 7/19 granular propellant coated with TiO2 high-molecular flame retardant are mixed in different proportions to perform a closed exploder experiment, the influence of the mixing proportion on the combustion performance of gunpowder is researched, and the experiment shows that the mixed charging structure can remarkably adjust the combustion performance of the gunpowder, adjust the pressure curve in a chamber and strengthen the pressure platform effect.
The effective means for improving the power of the electrothermal chemical cannon is to change the plasma jet parameters and the charging structure to improve the combustion and inner trajectory performance of gunpowder. From the technical data searched at present, an internal trajectory calculation model which can simultaneously consider the mixed charge structure and the plasma enhancement effect in the electrothermal chemical emission technology is not found.
Disclosure of Invention
The invention aims to solve the problem that only a single charge structure is considered in the existing electrothermal chemical emission inner trajectory model.
The technical solution for realizing the purpose of the invention is as follows: a method for calculating ballistic characteristic parameters in a mixed charge under plasma ignition comprises the following steps:
step 1, determining a pulse power supply input electric energy calculation model;
step 2, establishing a gunpowder burning rate formula under the plasma enhancement effect;
step 3, establishing a shape function of the gunpowder under the condition of mixing the two gunpowders;
step 4, establishing a projectile motion equation and an electrothermal chemical emission internal trajectory equation under mixed charging;
and 5, substituting the input electric energy into an electrothermal chemical launch internal trajectory equation, combining a powder burning rate formula and a shape function to carry out time propulsion solution, and finally obtaining an in-chamber pressure and a pellet speed curve in the process of discharging the pellet by using a pellet motion equation and the electrothermal chemical launch internal trajectory equation.
Further, step 1, determining the input electric energy of the pulse power supply according to the discharge energy and the electric energy utilization efficiency of the pulse power supply, wherein the specific calculation method comprises the following steps:
E0＝Etcpl
in the formula, E0For pulse power supply, EtDischarging energy for a pulsed power supply, cplThe electric energy utilization efficiency is improved; etDetermined by the discharge current and voltage conditions; c. CplThe experimental result of the closed exploder is obtained by the following steps:
Figure BDA0003533568920000021
Figure BDA0003533568920000022
wherein u (t) and i (t) are the discharge voltage and current of the pulse power supply, V0To seal the volume of the detonator, pmIs the maximum pressure in the explosive device, omega is the charge, alpha is the residual volume, c1Is the heat loss coefficient, f is the explosive power, and k is the specific heat ratio of the fuel gas.
Further, step 2, establishing a gunpowder burning rate formula under the plasma enhancement effect, wherein the specific method comprises the following steps:
Figure BDA0003533568920000023
wherein u is the burning rate of gunpowder, e1Is the arc thickness of gunpowder, Z is the relative burned thickness of gunpowder, u1Is a burning rate factor, n is an exponential factor, p is the pressure in the combustion chamber, betaeAs an electric power enhancing factor, PeIs the electrical power input to the generator.
Further, step 3, calculating a shape function of the gunpowder under the mixture of the two gunpowders, wherein the specific method comprises the following steps:
gunpowder shape function:
Figure BDA0003533568920000024
in the formula, N is the number of gunpowder types, the lower corner mark N represents the related parameters of the nth gunpowder,
Figure BDA0003533568920000025
the proportion of the burnt powder is the proportion of the burnt powder,
Figure BDA0003533568920000031
and wnThe burned proportion and the powder mass of the nth powder are obtained; modeling using porous gunpowder
Figure BDA0003533568920000032
The split point correlation should be considered in the correlation calculation:
Figure BDA0003533568920000033
in the formula, ZnThe relative burnt thickness, Z, of the n-th powderknIs the n-th gunpowderRelative burned thickness, chi, at full burn-out after splittingn、λn、μnThe parameter of the gunpowder type before the splitting of the nth gunpowder is shown, the lower corner mark s represents the related parameter after the splitting of the porous gunpowder, and chisn、λsnI.e. the parameter of the explosive type after the nth gunpowder is split, Zi、ZkiCalculated from the following formulas, respectively:
Figure BDA0003533568920000034
Figure BDA0003533568920000035
in the formula, enHalf of the burnt thickness of the n-th powder at a certain moment, e1nIs half of the original thickness of the gunpowder, rnThe equivalent inscribed circle radius in the section when the nth porous gunpowder is classified is related to the shape of the gunpowder;
the powder type parameters are calculated by the following formulas:
Figure BDA0003533568920000036
in the formula (I), the compound is shown in the specification,
Figure BDA0003533568920000037
2e12c denotes the charge width and 2a denotes the charge length, which are the original thicknesses of the charges.
Further, step 4, establishing a projectile motion equation and an electrothermal chemical launching internal trajectory equation under mixed charging, wherein the specific method comprises the following steps:
the motion equation of the projectile:
Figure BDA0003533568920000038
Figure BDA0003533568920000039
wherein v represents the projectile velocity, l is the barrel length, p is the projectile bottom pressure,
Figure BDA00035335689200000310
is the secondary work coefficient, S is the bore cross-sectional area, and m is the shot mass;
an electrothermal chemical emission internal trajectory equation under mixed charge:
neglecting the mass of the plasma jet, adding the plasma jet as an energy source term into the energy equation, E0The input electric energy of the pulse power supply is represented, and the basic equation of the trajectory in the electrothermal chemical emission is represented as follows:
Figure BDA0003533568920000041
in the formula, N is the number of gunpowder types,
Figure BDA0003533568920000044
indicating the reduction of the free volume of the chamber, fnIs the explosive power of the nth powder, wn is the mass of the nth powder,
Figure BDA0003533568920000042
the burned ratio of the nth powder and theta is the thermodynamic coefficient.
Further, step 5, substituting the input electric energy into an electrothermal chemical emission internal trajectory equation, combining a powder burning rate formula and a shape function to carry out time propulsion solution, and finally obtaining an in-bore pressure and a pellet speed curve in the process of discharging the pellet from the bore by using a pellet motion equation and the electrothermal chemical emission internal trajectory equation, wherein the specific method comprises the following steps:
step 51: constant determination
Determining the discharge voltage and current of a pulse power supply by combining an actual emission model, determining the quality of a shot, the sectional area of a bore, the volume of a charge chamber, the thermal coefficient, the secondary work coefficient, the starting pressure of the shot and the length of a gun barrel, and determining the arc thickness, the length and the width of selected gunpowder particles;
step 52: setting a calculation step size
Setting a calculation step length, and dividing the full trajectory into 100-200 points;
step 53: calculation of initial value
v0＝0；
Figure BDA0003533568920000043
l0＝0；
In the formula, the subscript 0 is the initial variable value, v0And l0I.e. the initial projectile velocity and the length of the projectile travel along the bore, psi0The ratio of the burnt powder at the beginning, delta is the powder charge density, rhoPIs the density of the powder, p0The initial pressure in the chamber is usually atmospheric pressure;
step 54: ballistic cycle calculation
Substituting the calculation result of the electric energy input by the pulse power supply into an electrothermal chemical launching inner trajectory equation, calculating the burnt thickness of the gunpowder by combining a gunpowder burning rate formula along with the time step advance, calculating the burnt proportion of the gunpowder by combining the gunpowder burning rate formula, substituting the calculation result into a projectile motion equation and an electrothermal chemical launching inner trajectory equation to solve the pressure in a chamber and the projectile speed, and finally outputting the pressure in the chamber and the projectile speed curve in the projectile discharging process.
A system for calculating the inner ballistic characteristic parameters of mixed charges under plasma ignition realizes the rapid calculation of an inner ballistic curve based on the inner ballistic calculation method considering the structure of the mixed charges and the enhancement effect of plasma.
A computer arrangement comprising a memory, a processor and a computer program stored on the memory and executable on the processor, the processor when executing the computer program, enabling a fast calculation of an inner ballistic curve based on the inner ballistic calculation method taking into account the mixed charge structure and the plasma enhancement effect.
A computer-readable storage medium, on which a computer program is stored which, when being executed by a processor, enables a fast calculation of an inner ballistic curve based on the inner ballistic calculation method taking into account the mixed charge structure and the plasma enhancement effect.
Compared with the prior art, the invention has the remarkable advantages that: 1) compared with other electrothermal chemical emission inner trajectory calculation models, the calculation method provided by the invention considers the influence of the mixed charge structure on the inner trajectory performance, and fills the blank of the existing model; 2) compared with the limitation brought by experimental observation, the invention can obtain a complete trajectory curve, characteristic point pressure and projectile velocity, and can more comprehensively analyze the change rule of the trajectory process in electrothermal chemical emission; 3) the projectile discharging speed under each proportion can be obtained by setting different mixing proportions of gunpowder, and compared with a conventional method for determining the optimal mixing proportion of a charging structure, the technical method provided by the invention can quickly obtain the optimal mixing proportion under the condition of ensuring effectiveness and precision.
Drawings
Figure 1 is a graph of current versus voltage discharge for a selected plasma generator of the present invention.
FIG. 2 is a schematic view of a 30mm electrothermal chemical device.
Figure 3 is a pressure image of single base powder charge and 15/7 granular powder at different mix ratios.
Figure 4 is a velocity image of single base gunpowder and 15/7 granular powder at different mix ratios.
Detailed Description
In order to make the objects, technical solutions and advantages of the present application more apparent, the present application is described in further detail below with reference to the accompanying drawings and embodiments. It should be understood that the specific embodiments described herein are merely illustrative of the present application and are not intended to limit the present application.
Based on a classical inner ballistic model, the mixed charge structure and the modification of the plasma enhancement effect are introduced, a calculation model of ballistic characteristic parameters in the mixed charge under plasma ignition is established, inner ballistic characteristic point parameters and ballistic curve change rules can be quickly and effectively calculated, and the problem that only a single charge structure is considered in the conventional electrothermal chemical emission inner ballistic model is solved. The method specifically comprises the following steps:
step 1, determining a pulse power supply discharge parameter.
Referring to fig. 1, the input power of the pulse power supply is affected by the discharge energy and the power utilization efficiency, which can be expressed as:
E0＝Etcpl
in the formula, E0For pulse power supply, EtDischarge of energy for a pulsed power supply, cplThe electric energy utilization efficiency is improved; etDetermined by the discharge current and voltage conditions; c. CplThe experimental result of the closed exploder is obtained by the following steps:
Figure BDA0003533568920000061
Figure BDA0003533568920000062
wherein u (t) and i (t) are the discharge voltage and current of the pulse power supply, V0For sealing the volume of the explosion chamber, pmIs the maximum pressure in the explosive device, omega is the charge, alpha is the residual volume, c1Is the heat loss coefficient, f is the explosive power, and k is the specific heat ratio of the fuel gas.
And 2, establishing a gunpowder burning speed formula under the plasma enhancement effect.
Introducing electric power burning rate enhancement factor beta on the basis of conventional gunpowder burning rate formulaeThe firing rate of the gunpowder under the action of the plasma can be expressed as:
Figure BDA0003533568920000063
wherein u is the burning rate of gunpowder, e1Is the arc thickness of gunpowder, Z is the relative burned thickness of gunpowder, u1Is a burning rate factor, n is an exponential factor, p is the pressure in the combustion chamber, betaeAs electrical power enhancing factor, PeIs the electrical power input to the generator.
Step 3, determining shape functions of the two gunpowders;
gunpowder shape function:
Figure BDA0003533568920000064
in the formula, N is the number of gunpowder types, the lower corner mark N represents the related parameters of the nth gunpowder,
Figure BDA0003533568920000065
the proportion of the burnt powder is the proportion of the burnt powder,
Figure BDA0003533568920000066
and wnThe burned proportion and the powder mass of the nth powder are obtained; modeling using porous gunpowder
Figure BDA0003533568920000067
The split point correlation should be considered in the correlation calculation:
Figure BDA0003533568920000068
in the formula, ZnThe relative burnt thickness, Z, of the n-th powderknIs the relative burned thickness, chi, of the nth powder after splitting and completely burningn、λn、μnThe parameter of the gunpowder type before the splitting of the nth gunpowder is shown, the lower corner mark s represents the related parameter after the splitting of the porous gunpowder, and chisn、λsnI.e. the parameter of the explosive type after the nth gunpowder is split, Zi、ZkiCalculated by the following formula, respectively:
Figure BDA0003533568920000071
Figure BDA0003533568920000072
in the formula, enHalf of the burnt thickness of the n-th powder at a certain moment, e1nHalf of the original thickness of the powder, rnThe equivalent inscribed circle radius in the cross section when the n-th porous powder is sorted is related to the shape of the powder.
The powder type parameters are calculated by the following formulas:
Figure BDA0003533568920000073
wherein
Figure BDA0003533568920000074
In the formula, 2e12c denotes the charge width and 2a denotes the charge length, which are the original thicknesses of the charges.
Selecting single-base tubular gunpowder and 15/7 granular gunpowder, wherein the parameters of the types of the single-base tubular gunpowder and the 15/7 granular gunpowder are as follows
TABLE 1 medicine type parameter table
Parameter(s) Single base tubular medicine Granular medicine
Type parameter χ i 1 0.7164
Drug type parameter lambda i 0 0.2246
Drug type parameter mui 0 -0.0269
Type parameter χ si 1 1.7616
Drug type parameter lambdasi 0 -0.4323
Step 4, establishing a projectile motion equation and an electrothermal chemical emission internal trajectory equation;
the motion equation of the projectile:
Figure BDA0003533568920000075
Figure BDA0003533568920000076
wherein v represents the projectile velocity, l is the barrel length, p is the projectile bottom pressure,
Figure BDA0003533568920000077
is the secondary work coefficient, S is the bore cross-sectional area, and m is the projectile mass.
An electrothermal chemical emission internal trajectory equation under mixed charge:
ignoring the mass of the plasma jet, adding the plasma jet as an energy source term into the energy equation, E0The input electric energy of the pulse power supply is represented, and the basic equation of the trajectory in the electrothermal chemical emission is represented as follows:
Figure BDA0003533568920000081
in the formula, N is the number of gunpowder types,
Figure BDA0003533568920000082
indicating the reduction of the free volume of the chamber, fnThe explosive power of the nth gunpowder, wnThe mass of the nth gunpowder is,
Figure BDA0003533568920000083
the burned ratio of the nth powder and theta is the thermodynamic coefficient.
Step 5, substituting the input electric energy into an electrothermal chemical launch internal trajectory equation, combining a gunpowder burning rate formula and a shape function to carry out time propulsion solution, and finally obtaining an in-chamber pressure and a projectile velocity curve in the projectile discharging process by utilizing a projectile motion equation and the electrothermal chemical launch internal trajectory equation;
step 51: constant determination
And determining the discharge voltage and current of a pulse power supply by combining an actual emission model, determining the quality of the projectile, the sectional area of a bore, the volume of a charge chamber, the thermal coefficient, the secondary power coefficient, the starting pressure of the projectile and the length of a gun barrel, and determining the parameters such as the arc thickness of the selected gunpowder particles, the length of the gunpowder, the width of the gunpowder and the like.
Referring to fig. 2, the electrothermal chemical launcher with 30mm caliber is shown in the figure, the length of the cannon body tube is 2.75m, the volume of the medicine chamber is 356mL, the starting pressure of the projectile is 30MPa, and the mass of the projectile is 68.6 g.
Step 52: setting a calculation step size
In general, the calculation requirement can be satisfied by dividing the full trajectory into 100-200 points.
Step 53: calculation of initial value
v0＝0；
Figure BDA0003533568920000084
l0＝0；
In the formula, the subscript 0 is the initial variable value, v0And l0I.e. initial ammunitionThe speed of the projectile and the length of the projectile in the bore, psi0The ratio of the burnt powder at the beginning, delta is the powder charge density, rhoPIs the density of the powder, p0The initial pressure in the bore is usually atmospheric pressure.
Step 54: ballistic cycle calculation
Substituting the calculation result of the electric energy input by the pulse power supply into an electrothermal chemical launching inner trajectory equation, calculating the burnt thickness of the gunpowder by combining a gunpowder burning rate formula along with the time step advance, calculating the burnt ratio of the gunpowder by combining the gunpowder burning rate formula, substituting the calculation result into a projectile motion equation and an electrothermal chemical launching inner trajectory equation to solve the pressure in a chamber and the speed of the projectile, and finally outputting the pressure in the chamber and the speed curve of the projectile in the process of the projectile leaving the chamber.
As shown in fig. 3 and 4, in the 30mm electrothermal chemical launching device, when the single-base tubular medicine and the 15/7 granular medicine are mixed in different proportions, the initial speed of the projectile tends to increase and then decrease along with the increase of the proportion of the mixed 15/7 granular medicine, and the mixing proportion of 15% can enable the projectile to have the maximum outlet speed.
The invention further provides a system for calculating the inner ballistic characteristic parameters of the mixed charge under plasma ignition, and the inner ballistic characteristic point parameters and the ballistic curve are quickly calculated based on the inner ballistic calculation method considering the mixed charge structure and the plasma enhancement effect.
A computer arrangement comprising a memory, a processor and a computer program stored on the memory and executable on the processor, the processor when executing the computer program enabling a fast calculation of internal ballistic characteristic point parameters and ballistic curve based on the internal ballistic calculation method taking into account mixed charge structure and plasma enhancement.
A computer-readable storage medium, on which a computer program is stored which, when being executed by a processor, enables a fast calculation of internal ballistic characteristic point parameters and ballistic curves based on the internal ballistic calculation method taking into account mixed charge structure and plasma enhancement.
The technical features of the above embodiments can be arbitrarily combined, and for the sake of brevity, all possible combinations of the technical features in the above embodiments are not described, but should be considered as the scope of the present specification as long as there is no contradiction between the combinations of the technical features.
The above-mentioned embodiments only express several embodiments of the present application, and the description thereof is more specific and detailed, but not construed as limiting the scope of the invention. It should be noted that, for a person skilled in the art, several variations and modifications can be made without departing from the concept of the present application, which falls within the scope of protection of the present application. Therefore, the protection scope of the present patent shall be subject to the appended claims.

Claims (9)

1. A method for calculating ballistic characteristic parameters in a mixed charge under plasma ignition comprises the following steps:
step 1, determining a pulse power supply input electric energy calculation model;
step 2, establishing a gunpowder burning rate formula under the plasma enhancement effect;
step 3, establishing a shape function of the gunpowder under the condition of mixing the two gunpowders;
step 4, establishing a projectile motion equation and an electrothermal chemical emission internal trajectory equation under mixed charging;
and 5, substituting the input electric energy into an electrothermal chemical launch internal trajectory equation, combining a powder burning rate formula and a shape function to carry out time propulsion solution, and finally obtaining an in-chamber pressure and a pellet speed curve in the process of discharging the pellet by using a pellet motion equation and the electrothermal chemical launch internal trajectory equation.
2. The method for calculating the ballistic characteristic parameters of the mixed charge under plasma ignition according to claim 1, wherein in the step 1, the input electric energy of the pulse power supply is determined according to the discharge energy and the electric energy utilization efficiency of the pulse power supply, and the specific calculation method comprises the following steps:
E0＝Etcpl
in the formula, E0For pulse power supply, EtDischarging energy for a pulsed power supply, cplThe electric energy utilization efficiency is improved; etDetermined by the discharge current and voltage conditions; c. CplThe experimental result of the closed exploder is obtained by the following steps:
Figure FDA0003533568910000011
Figure FDA0003533568910000012
wherein u (t) and i (t) are the discharge voltage and current of the pulse power supply, V0For sealing the volume of the explosion chamber, pmIs the maximum pressure in the explosive device, omega is the charge, alpha is the residual volume, c1Is the heat loss coefficient, f is the explosive power, and k is the specific heat ratio of the fuel gas.
3. The method for calculating the ballistic characteristic parameters of the mixed charge under plasma ignition according to claim 1, wherein the step 2 of establishing a powder burning rate formula under the plasma enhancement effect comprises the following specific steps:
Figure FDA0003533568910000013
wherein u is the burning rate of gunpowder, e1Arc thickness of gunpowder, Z relative burned thickness of gunpowder, u1Is a burning rate factor, n is an exponential factor, p is the pressure in the combustion chamber, betaeAs electrical power enhancing factor, PeIs the electrical power input to the generator.
4. The method for calculating the ballistic characteristic parameters of the mixed charge under plasma ignition according to claim 1, wherein the step 3 of calculating the shape function of the two types of gunpowder under mixed gunpowder comprises the following specific steps:
gunpowder shape function:
Figure FDA0003533568910000021
in the formula, N is the number of gunpowder types, the lower corner mark N represents the related parameters of the nth gunpowder,
Figure FDA0003533568910000022
the proportion of the burnt powder is the proportion of the burnt powder,
Figure FDA0003533568910000023
and wnThe burned proportion and the powder mass of the nth powder are obtained; modeling using porous gunpowder
Figure FDA0003533568910000024
The split point correlation should be considered in the correlation calculation:
Figure FDA0003533568910000025
in the formula, ZnThe relative burnt thickness, Z, of the n-th powderknIs the relative burned thickness, chi, of the nth powder after splitting and completely burnedn、λn、μnThe parameter of the gunpowder type before the splitting of the nth gunpowder is shown, the lower corner mark s represents the related parameter after the splitting of the porous gunpowder, and chisn、λsnI.e. the parameter of the explosive type after the nth gunpowder is split, Zi、ZkiCalculated by the following formula, respectively:
Figure FDA0003533568910000026
Figure FDA0003533568910000027
in the formula，enHalf of the burnt thickness of the n-th powder at a certain moment, e1nIs half of the original thickness of the gunpowder, rnThe equivalent inscribed circle radius in the section when the nth porous gunpowder is classified is related to the shape of the gunpowder;
the powder type parameters are calculated by the following formulas:
Figure FDA0003533568910000028
in the formula (I), the compound is shown in the specification,
Figure FDA0003533568910000029
2e12c denotes the charge width and 2a denotes the charge length, which are the original thicknesses of the charges.
5. The method for calculating the internal ballistic characteristic parameters of the mixed charge under plasma ignition according to claim 1, wherein in the step 4, a projectile motion equation and an internal ballistic equation of electrothermal chemical emission under the mixed charge are established, and the specific method is as follows:
the motion equation of the projectile:
Figure FDA0003533568910000031
Figure FDA0003533568910000032
wherein v represents the projectile velocity, l is the barrel length, p is the projectile bottom pressure,
Figure FDA0003533568910000033
is the secondary work coefficient, S is the bore cross-sectional area, and m is the shot mass;
an electrothermal chemical emission internal trajectory equation under mixed charge:
ignoring the mass of the plasma jet, treating the plasma jet as aEnergy source terms added to the energy equation, E0The input electric energy of the pulse power supply is expressed, and the basic equation of the inner trajectory of the electrothermal chemical emission is expressed as follows:
Figure FDA0003533568910000034
in the formula, N is the number of gunpowder types,
Figure FDA0003533568910000035
showing the reduction of the free volume of the chamber, fnThe explosive power of the nth powder, wnThe mass of the nth gunpowder is,
Figure FDA0003533568910000036
the nth powder is burnt, and theta is a thermal coefficient.
6. The method for calculating the internal ballistic characteristic parameters of mixed charge under plasma ignition according to claim 1, wherein in step 5, the input electric energy is substituted into an electrothermal chemical projectile internal ballistic equation, time propulsion solution is performed by combining a powder burning rate formula and a shape function, and an in-chamber pressure and a projectile velocity curve in the projectile discharging process are finally obtained by using a projectile motion equation and the electrothermal chemical projectile internal ballistic equation, and the specific method is as follows:
step 51: constant determination
Determining the discharge voltage and current of a pulse power supply by combining an actual emission model, determining the quality of a shot, the sectional area of a bore, the volume of a charge chamber, the thermal coefficient, the secondary work coefficient, the starting pressure of the shot and the length of a gun barrel, and determining the arc thickness, the length and the width of selected gunpowder particles;
step 52: setting a calculation step size
Setting a calculation step length, and dividing the whole trajectory into 100-200 points;
step 53: calculation of initial value
v0＝0；
Figure FDA0003533568910000037
l0＝0；
In the formula, the lower subscript 0 is the initial variable value, v0And l0I.e. the initial projectile velocity and the length of the projectile travel along the bore, psi0The ratio of the burnt powder at the beginning, delta is the powder charge density, rhoPIs the density of the powder, p0The initial pressure in the chamber is usually atmospheric pressure;
step 54: ballistic cycle calculation
Substituting the calculation result of the electric energy input by the pulse power supply into an electrothermal chemical launching inner trajectory equation, calculating the burnt thickness of the gunpowder by combining a gunpowder burning rate formula along with the time step advance, calculating the burnt ratio of the gunpowder by combining the gunpowder burning rate formula, substituting the calculation result into a projectile motion equation and an electrothermal chemical launching inner trajectory equation to solve the pressure in a chamber and the speed of the projectile, and finally outputting the pressure in the chamber and the speed curve of the projectile in the process of the projectile leaving the chamber.
7. A system for calculating the inner ballistic characteristic parameters of a mixed charge under plasma ignition, characterized in that the rapid calculation of the inner ballistic curve is achieved on the basis of the inner ballistic calculation method taking into account the structure of the mixed charge and the plasma enhancement according to any one of claims 1 to 6.
8. A computer arrangement comprising a memory, a processor and a computer program stored on the memory and executable on the processor, when executing the computer program, enabling a fast calculation of an inner ballistic curve based on the inner ballistic calculation method taking into account the mixed charge structure and the plasma enhancement as claimed in any one of claims 1-6.
9. A computer-readable storage medium, on which a computer program is stored which, when being executed by a processor, enables a fast calculation of an inner ballistic curve based on the inner ballistic calculation method taking into account the mixed charge structure and the plasma-enhanced action of any one of claims 1-6.
CN202210216931.XA 2022-03-06 2022-03-06 Method for calculating ballistic characteristic parameters in mixed charge under plasma ignition Pending CN114639450A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202210216931.XA CN114639450A (en) 2022-03-06 2022-03-06 Method for calculating ballistic characteristic parameters in mixed charge under plasma ignition

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202210216931.XA CN114639450A (en) 2022-03-06 2022-03-06 Method for calculating ballistic characteristic parameters in mixed charge under plasma ignition

Publications (1)

Publication Number Publication Date
CN114639450A true CN114639450A (en) 2022-06-17

Family

ID=81947052

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202210216931.XA Pending CN114639450A (en) 2022-03-06 2022-03-06 Method for calculating ballistic characteristic parameters in mixed charge under plasma ignition

Country Status (1)

Country Link
CN (1) CN114639450A (en)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115563730A (en) * 2022-09-23 2023-01-03 南京理工大学 Calculation method of trajectory data elements in light high-low pressure gun considering backseat movement

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115563730A (en) * 2022-09-23 2023-01-03 南京理工大学 Calculation method of trajectory data elements in light high-low pressure gun considering backseat movement
CN115563730B (en) * 2022-09-23 2023-10-31 南京理工大学 Calculation method of trajectory data in light high-low pressure artillery considering backseat movement

Similar Documents

Publication Publication Date Title
US3253511A (en) Launching process and apparatus
CN107025325A (en) A kind of ionization seed optimization method for making powder gas generate plasma
CN114639450A (en) Method for calculating ballistic characteristic parameters in mixed charge under plasma ignition
CN110287599A (en) A kind of propellant powder generates the construction method of plasma rule model
Cao et al. Design optimization for a launching system with novel structure
JP2000055597A (en) Launching powder
Guo et al. Numerical Simulation of Hypervelocity Launch Process of Cased Telescoped Ammunition Armor-piercing Projectile
Wang et al. Design optimization and simulation of interior ballistic measurement technology based on parameters of large-caliber launch fuel
Rogaev et al. Rifle shot upgrading by using model paste-like propellant
Sidorov et al. Analysis of the grain powder combustion under electrothermal chemical technology of shot at different initial temperatures
Chen et al. Experimental research on the launching system of auxiliary charge with filter cartridge structure
Surma Recoilless Gun System as a Particular Form of General Interior Ballistics Model of Gun Propellant Systems
Li et al. Numerical simulation of the two-phase flow of a constant-pressure rear-spray low recoil weapon with two chambers
Iorga et al. Design and Testing of an Unguided Rocket with Thermobaric Warhead for Multiple Launcher System
CAO et al. Interior Ballistic Modeling and Charge Structure Optimization of Cased Telescoped Ammunition
RU2807726C1 (en) Modernized malevanny ammunition and method of its manufacturing
Ishchenko et al. Determination of the propellant combustion law under ballistic experiment conditions
Rashad et al. Numerical simulation of interior ballistics for large caliber guided projectile naval gun
Zecevic et al. Specific Design Features of Solid Propellant Rocket Motors for Shoulder-Launched Weapon Systems
Cao et al. Effect of ignition delay runaway of traveling charge on gun firing safety
KR102269204B1 (en) Projectile containing ramjet engine
RU2150080C1 (en) Rocket
RU2812632C1 (en) Projectile
Ma et al. Study on micro recoil mechanism of the weapon with a nozzle and two chambers separated by a partition
Xie et al. A fast calculation model for interior ballistics simulation of a recoilless launch system with high chamber pressure

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination