[go: up one dir, main page]

CN114813835A - Evaluation method for rapid burning of explosive - Google Patents

Evaluation method for rapid burning of explosive Download PDF

Info

Publication number
CN114813835A
CN114813835A CN202210379372.4A CN202210379372A CN114813835A CN 114813835 A CN114813835 A CN 114813835A CN 202210379372 A CN202210379372 A CN 202210379372A CN 114813835 A CN114813835 A CN 114813835A
Authority
CN
China
Prior art keywords
model
test
rapid
fuel
charge
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.)
Granted
Application number
CN202210379372.4A
Other languages
Chinese (zh)
Other versions
CN114813835B (en
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.)
Beijing Institute of Technology BIT
Original Assignee
Beijing Institute of Technology BIT
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 Beijing Institute of Technology BIT filed Critical Beijing Institute of Technology BIT
Priority to CN202210379372.4A priority Critical patent/CN114813835B/en
Publication of CN114813835A publication Critical patent/CN114813835A/en
Application granted granted Critical
Publication of CN114813835B publication Critical patent/CN114813835B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N25/00—Investigating or analyzing materials by the use of thermal means
    • G01N25/50—Investigating or analyzing materials by the use of thermal means by investigating flash-point; by investigating explosibility
    • G01N25/54—Investigating or analyzing materials by the use of thermal means by investigating flash-point; by investigating explosibility by determining explosibility
    • F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
    • F42—AMMUNITION; BLASTING
    • F42B—EXPLOSIVE CHARGES, e.g. FOR BLASTING, FIREWORKS, AMMUNITION
    • F42B35/00—Testing or checking of ammunition
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01J—MEASUREMENT OF INTENSITY, VELOCITY, SPECTRAL CONTENT, POLARISATION, PHASE OR PULSE CHARACTERISTICS OF INFRARED, VISIBLE OR ULTRAVIOLET LIGHT; COLORIMETRY; RADIATION PYROMETRY
    • G01J5/00—Radiation pyrometry, e.g. infrared or optical thermometry
    • G01J5/0014—Radiation pyrometry, e.g. infrared or optical thermometry for sensing the radiation from gases, flames
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N25/00—Investigating or analyzing materials by the use of thermal means
    • G01N25/20—Investigating or analyzing materials by the use of thermal means by investigating the development of heat, i.e. calorimetry, e.g. by measuring specific heat, by measuring thermal conductivity
    • G01N25/22—Investigating or analyzing materials by the use of thermal means by investigating the development of heat, i.e. calorimetry, e.g. by measuring specific heat, by measuring thermal conductivity on combustion or catalytic oxidation, e.g. of components of gas mixtures
    • G—PHYSICS
    • G06—COMPUTING OR CALCULATING; COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/20—Design optimisation, verification or simulation

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Chemical & Material Sciences (AREA)
  • Immunology (AREA)
  • Pathology (AREA)
  • General Engineering & Computer Science (AREA)
  • Theoretical Computer Science (AREA)
  • General Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Biochemistry (AREA)
  • Analytical Chemistry (AREA)
  • Health & Medical Sciences (AREA)
  • Geometry (AREA)
  • Combustion & Propulsion (AREA)
  • Chemical Kinetics & Catalysis (AREA)
  • Spectroscopy & Molecular Physics (AREA)
  • Evolutionary Computation (AREA)
  • Computer Hardware Design (AREA)
  • Investigating Or Analyzing Materials Using Thermal Means (AREA)

Abstract

本发明公开了一种炸药快速烤燃的评估方法，具体步骤为：根据快速烤燃试验系统建立含空气域的快速烤燃计算模型进行仿真；所述含空气域的快速烤燃计算模型包括：试验弹模型、空气域、燃料入口及燃料出口；建立不含空气域的快速烤燃计算模型，所述不含空气域的快速烤燃计算模型仅包括试验弹模型；最终得到热通量边界与对应的装药模型的直径之间的散点对应关系；根据装药模型的任一直径及计算得到的热通量边界，将该热通量边界加载到试验弹模型的外壁面进行计算，获取装药模型内部温度变化、点火时间及点火温度，对炸药快速烤燃特性进行评估；本发明能够对现有的含空气域快速烤燃计算进行优化，实现无空气域炸药快速烤燃的评估。

Figure 202210379372

The invention discloses an evaluation method for rapid roasting of explosives. The specific steps are: establishing a rapid roasting calculation model containing air domain according to a rapid roasting test system for simulation; the rapid roasting calculation model containing air domain comprises: Test bomb model, air domain, fuel inlet and fuel outlet; establish a rapid roasting calculation model without air domain, and the rapid roasting calculation model without air domain only includes the test bomb model; finally obtain the heat flux boundary and Scatter correspondence between the diameters of the corresponding charge model; according to any diameter of the charge model and the calculated heat flux boundary, load the heat flux boundary onto the outer wall of the test projectile model for calculation, and obtain The internal temperature change, ignition time and ignition temperature of the charge model are used to evaluate the rapid burning characteristics of explosives; the present invention can optimize the existing rapid burning calculation in the air-containing region to realize the evaluation of the rapid burning of the explosive without air.

Figure 202210379372

Description

Evaluation method for rapid burning of explosive
Technical Field
The invention belongs to the technical field of explosive safety evaluation, and particularly relates to an evaluation method for rapid roasting and burning of an explosive.
Background
The modern war has higher and higher requirements on the thermal safety of weapons and ammunitions, and the firearms and ammunitions must be absolutely safe in the processes of manufacturing, storing, transporting and using, cannot cause serious accidents such as ignition and explosion and cannot reduce the safety of the firearms and ammunitions in the service life cycle. Therefore, the research on the safety of the ammunition under thermal stimulation has very important significance.
The fire test is an important means for testing and evaluating the thermal safety of ammunition and is divided into fast fire and slow fire according to different heating environments. The rapid fire-burning mainly simulates the thermal response characteristic of ammunition in a fire environment, and tests and numerical simulation are important methods for researching the rapid fire-burning of ammunition. Because the burning test needs a large test field, the test result is greatly influenced by the environmental conditions, and the test device is mainly used for the test research of prototype bullets.
With the development of a numerical simulation technology, numerical simulation gradually becomes an important research method for rapid burning, more researches are carried out in China, wherein the dynamite burning test and the numerical simulation are carried out based on a mass _ flux method [ J ] weapon equipment engineering newspaper, 2020,41(8):1-6 ] based on RBOE explosives, and a calculation model comprising a projectile body, an air domain, a fuel inlet and a fuel outlet is established on the basis of a physical diagram of a burning test device, as shown in figure 1. And loading the programmed linear temperature-time historical curve and mass flow rate to a fuel inlet by adopting a mass flow inlet boundary condition, reacting in a charging area according to an autothermal reaction rule, and applying a reaction kinetic equation to the charging area through a subprogram. The projectile was then heated by high temperature gas radiation and thermal convection, the temperature history of the test and simulated flame environments is shown in fig. 2, and the mass flow rate was adjusted so that the calculated ignition time was substantially consistent with the test ignition time.
The existing method for calculating the fast burning value has the defects that:
(1) the mass flow rate and the linear temperature-time history are loaded to a fuel inlet, and the calculation is less consistent with the temperature history curve of the test environment in a temperature rise section, so that the flame temperature growth process cannot be well described.
(2) The air domain needs to be established in the fast-baking numerical calculation, and in order to simulate the flame combustion environment, the size of the air domain needs to be far larger than that of the projectile body, so that the calculation amount is multiplied, and the calculation efficiency is reduced.
(3) The existing numerical calculation of the quick air-containing domain combustion-baking value is not directly loaded with the temperature on the surface of an elastomer, but is heated through the radiation and convection action of high-temperature fuel gas, the calculation result is greatly influenced by factors such as grid size and air domain size, and the calculation precision is poor.
Disclosure of Invention
In view of this, the invention provides an evaluation method for explosive fast-burning, which can optimize the existing air-containing domain fast-burning calculation and realize the evaluation of explosive fast-burning without an air domain.
The invention is realized by the following technical scheme:
an evaluation method for fast burning of explosives comprises the following specific steps:
step one, building a rapid burning test system and carrying out a rapid burning test;
establishing a rapid combustion-supporting calculation model containing an air domain according to the rapid combustion-supporting test system for simulation; the fast-burning calculation model of the air-containing domain comprises the following steps: the test bomb comprises a test bomb model, an air domain, a fuel inlet and a fuel outlet;
calculating the charge models in the test bomb models in the fast combustion calculation model containing the air domain in different sizes, namely, calculating to obtain the ignition time A of the charge models under different charge model diameters by changing the diameters of the charge models;
establishing a rapid combustion-baking calculation model without an air domain, wherein the rapid combustion-baking calculation model without the air domain only comprises a test bomb model; the outer wall surface of the cartridge case model of the test cartridge model adopts a heat flux boundary as a heating boundary; according to the calculation model of the rapid combustion without the air domain, under the condition that the diameter of the calculation model of the charge is the same as that of the charge model in the third step, calculating the ignition time B of the charge model, and adjusting the size of a heat flux boundary to enable the ignition time A in the third step to be equal to the ignition time B; finally, obtaining the corresponding relation of the heat flux boundary and the corresponding diameter of the charging model, and fitting the scattering point into a curve;
and fifthly, loading the heat flux boundary to the outer wall surface of the test bomb model for calculation according to any diameter of the explosive charging model and the heat flux boundary obtained by calculation according to the curve obtained in the fourth step, obtaining the internal temperature change, the ignition time and the ignition temperature of the explosive charging model, and evaluating the fast burning characteristics of the explosive.
Further, in step one, the testing system comprises: the device comprises a fuel oil pool, a test bomb, a thermocouple, a bracket, a suspension rod, an ignition device and a temperature recorder; the test bomb comprises: the cartridge case and the explosive filled in the cartridge case; the test bomb is fixed on the support through the suspension rod, the fuel pool is positioned under the test bomb, the distance between the lowest part of the test bomb and the upper surface of fuel oil in the fuel pool is 300mm, and the fuel oil in the fuel pool is connected with the ignition device; the thermocouple is arranged outside the test bomb, the thermocouple is positioned on an axis extension line of the test bomb, the closest distance between the thermocouple and the outer end face of the test bomb is 50mm, and the temperature recorder is electrically connected with the thermocouple;
when a rapid burning test is carried out, the ignition device ignites the fuel oil, the fuel oil in the fuel oil pool is combusted to generate flame and high-temperature fuel gas, a flame field is formed around the test bomb, the charge of the test bomb is heated in the flame field, and when the temperature of the flame field is higher than the burning point of the charge of the test bomb, the charge is ignited when the highest temperature of the charge of the test bomb is higher than the burning point; the thermocouple measures the ambient temperature of the flame field where the test bomb is located when fuel oil in the fuel oil pool is combusted, and the temperature recorder records the change of the ambient temperature measured by the thermocouple along with time to obtain the change curve of the ambient temperature along with time.
Further, in step two, the test bomb model is identical to the test bomb, and includes: the shell model and a charging model filled in the shell model; the fuel inlet is high-temperature fuel gas generated by fuel oil combustion of the fuel oil pool, the boundary condition of the fuel inlet is a mass flow inlet boundary, and the mass flow inlet boundary adopts a temperature piecewise function and mass flow rate to fit the test environment temperature; the distance between the lowest part of the test bomb model and the fuel inlet is the same as the distance between the lowest part of the test bomb and the upper surface of the fuel oil pool; the boundary condition of the fuel outlet is a pressure outlet boundary; and simulating the condition of a flame field generated by fuel oil combustion of the fuel oil pool by adjusting the fuel gas mass flow rate and the temperature piecewise function, calculating to obtain a change curve of the environmental temperature of the flame field where the test bomb model is located along with time, and if the difference value of the change curve of the environmental temperature of the flame field where the test bomb model is located along with time and the change curve of the environmental temperature obtained by the test along with time at the same time point is in a set range, indicating that the established calculation model can describe the flame field of the rapid roasting test system.
Further, in the fourth step, a single-phase exponential decay function in Origin software is adopted to fit the corresponding dispersion relation between the heat flux boundary and the diameter of the corresponding charge model.
Further, in the fourth step, Matlab software is adopted to fit the dispersion point correspondence between the heat flux boundary and the diameter of the corresponding charge model.
Further, the scatter point correspondence between the heat flux boundary and the diameter of the corresponding charge model is fit to the following expression after the curve:
Figure BDA0003591616450000031
where Φ is the heat flux boundary and D is the diameter of the charge model.
Further, in the fourth step, the scatter point correspondence between the heat flux boundary and the diameter of the corresponding charge model is fit to the following expression after the curve:
Figure BDA0003591616450000032
where Φ is the heat flux boundary and D is the diameter of the charge model.
Further, before the third step, grid division needs to be performed on the rapid combustion-baking calculation model containing the air domain established in the second step, the grid size of the divided air domain is 100mm, and the grid size of the divided test bomb model is 2 mm.
Has the advantages that: the method for evaluating the rapid burning of the explosive provided by the invention optimizes the existing rapid burning calculation method on one hand and provides a new rapid burning calculation method on the other hand.
(1) For the existing fast burning calculation method of the air-containing domain, the temperature piecewise function and the mass flow rate are adopted for fitting the test environment temperature at the mass flow inlet boundary, namely, the temperature piecewise function is adopted for fitting the test environment temperature and the temperature piecewise function is combined with the gas mass flow rate to obtain the inlet boundary condition of the fast burning calculation model of the air-containing domain, so that the problem that the environment temperature calculated by the existing calculation method is poor in coincidence with the test environment temperature is solved.
(2) According to the method, a rapid burning calculation model only considering the test bomb is established, a heat flux boundary is used as a boundary condition for rapid burning numerical calculation, namely the heat flux boundary is used as a heating boundary, so that the calculation model is simplified, and the efficiency of rapid burning numerical simulation calculation is improved; the problems that an existing rapid combustion-baking calculation model containing an air domain is large in calculation amount, complex in influence factors and the like are solved.
(3) The method adopts the single-phase exponential decay function in Origin software to fit the relation between the heat flux boundary and the diameter of the charge model to obtain the relation between the heat flux boundary and the diameter of the charge model, so that the test and calculation quantity can be reduced, the proper heating boundary can be selected according to the size of the charge model to carry out fast burning calculation, and the numerical calculation efficiency of fast burning is greatly improved.
(4) The invention carries out grid convergence analysis, the air domain adopts a large-size grid, the test bomb model adopts a small-size grid, the grids with different sizes are calculated, and the calculation convergence is ensured, thereby improving the calculation efficiency on the premise of ensuring the calculation precision.
Drawings
FIG. 1 is a diagram of a computational model of an air-containing domain in the background art;
FIG. 2 is a graph of environmental temperature for testing and simulation in the background art;
FIG. 3 is a block diagram of a rapid ignition test system according to the present invention;
FIG. 4 is a graph of the environmental temperature of the test and simulation of the present invention over time;
FIG. 5 is a computational model of a fast-burn air-containing domain of the present invention;
fig. 6 is a graph of mesh convergence analysis, (a) is a graph of temperature versus time for different under-mesh charge models for the air domain (b) is a graph of temperature versus time for different under-mesh charge models for the test bomb model;
FIG. 7 is a computational model of a fast-fire without air domains according to the present invention;
FIG. 8 is a graph of the relationship between the heat flux boundary of the present invention and the diameter of the corresponding charge pattern;
the method comprises the following steps of 1-fuel inlet, 2-test bomb model, 3-fuel outlet, 4-air domain, 5-cartridge case model, 6-charging model, 7-outer wall surface, 8-fuel pool, 9-thermocouple, 10-bracket, 11-suspension rod, 12-ignition device, 13-temperature recorder, 14-test bomb, 15-cartridge case and 16-charging.
Detailed Description
The invention is described in detail below by way of example with reference to the accompanying drawings.
The embodiment provides an evaluation method for rapid burning of an explosive, which comprises the following specific steps:
step 1, building a rapid burning test system and developing a rapid burning test, as shown in fig. 3, the test system comprises: the device comprises a fuel oil pool 8, a test bomb 14, a thermocouple 9, a bracket 10, a suspension rod 11, an ignition device 12 and a temperature recorder 13; the test bomb 14 includes: the bullet comprises a bullet shell 15 and a charge 16 filled in the bullet shell 15, wherein the length-diameter ratio of the charge 16 is 2, and the wall thickness of the bullet shell 15 is 6 mm; the test bomb 14 is fixed on the support 10 through the suspension rod 11, the fuel pool 8 is located right below the test bomb 14, the distance between the lowest position of the test bomb 14 and the upper surface of fuel of the fuel pool 8 is 300mm, and the fuel in the fuel pool 8 is connected with the ignition device 12; the thermocouple 9 is positioned and installed outside the test bomb 14 through a mounting rack and rock wool, the thermocouple 9 is located on an axis extension line of the test bomb 14, the closest distance between the thermocouple 9 and the outer end face of the test bomb 14 is 50mm, and the temperature recorder 13 is electrically connected with the thermocouple 9;
when a rapid burning test is carried out, the ignition device 12 ignites the fuel oil, the fuel oil in the fuel oil pool 8 is burnt to generate flame and high-temperature fuel gas, a flame field is formed around the test bomb 14, the temperature of the charge 16 of the test bomb 14 is increased in the flame field, and when the highest temperature of the charge 16 of the test bomb 14 is higher than the burning point of the charge 16, the charge 16 is ignited; the thermocouple 9 measures the ambient temperature of the flame field where the test bomb 14 is located when the fuel oil in the fuel oil pool 8 is combusted, and the temperature recorder 13 records the change of the ambient temperature measured by the thermocouple 9 along with time to obtain the change curve of the tested ambient temperature along with time, as shown by the solid line in fig. 4;
step 2, establishing a rapid combustion-baking calculation model containing an air domain according to the rapid combustion-baking test system for simulation, as shown in fig. 5, wherein the rapid combustion-baking calculation model containing the air domain comprises: the test bomb model 2, the air domain 4, the fuel inlet 1 and the fuel outlet 3; the test bomb model 2 is identical to the test bomb 14 and comprises the following components: a shell model 5 and a charging model 6 filled in the shell model 5; the fuel inlet 1 is high-temperature fuel gas generated by fuel oil combustion of the fuel oil pool 8, the boundary condition of the fuel inlet 1 is a mass flow inlet boundary, the fuel gas mass flow rate of the fuel inlet 1 is specified, the boundary temperature is calculated by adopting a piecewise function, namely the temperature piecewise function to fit the test environment temperature, and finally the temperature piecewise function and the mass flow rate are adopted to fit the test environment temperature on the mass flow inlet boundary; the distance between the lowest part of the test bomb model 2 and the fuel inlet 1 is the same as the distance between the lowest part of the test bomb 14 and the upper surface of the fuel oil pool 8, and the distances are 300 mm; the boundary condition of the fuel outlet 3 is a pressure outlet boundary; simulating the condition of a flame field generated by fuel oil combustion of the fuel oil pool 8 by adjusting the fuel gas mass flow rate and the temperature piecewise function, and calculating to obtain a change curve of the environmental temperature of the flame field where the test bomb model 2 is located along with time, as shown by a dotted line in fig. 4, if the dotted line is basically consistent with the change curve of the environmental temperature along with time (namely a solid line in fig. 4) (namely the difference value of the environmental temperatures of the dotted line and the solid line at the same time point is in a set range), it indicates that the established calculation model can describe the flame field of the rapid combustion test system;
step 3, carrying out grid division on the rapid combustion-baking calculation model containing the air domain established in the step 2, and carrying out grid convergence analysis; because the mesh size of the fast combustion-baking calculation model containing the air domain has a large influence on the calculation result, the mesh convergence calculation needs to be performed on different areas, the calculation of the mesh sizes of 150mm, 100mm and 70mm is performed on the air domain 4, and the ignition time of the charge model 6 is gradually delayed and has a convergence trend along with the reduction of the mesh size, as shown in fig. 6 (a); the calculation of the grid size of the test bomb model 2 to be 2.5mm, 2mm and 1.5mm is respectively carried out, the ignition time of the explosive charging model 6 is gradually delayed along with the reduction of the grid size, and the explosive charging model has a convergence trend, as shown in fig. 6 (b); in order to balance the calculation efficiency and the calculation precision, the size of the grid for dividing the air domain 4 is 100mm, and the size of the grid for dividing the test bomb model 2 is 2 mm;
step 4, calculating different sizes of the charge models 6 in the test bomb model 2 according to the quick combustion-baking calculation model of the air-containing domain after the grid division is finished, namely, obtaining ignition time A of the charge models 6 under different diameters D of the charge models 6 by changing the diameter D of the charge models 6;
step 5, establishing a rapid combustion-baking calculation model without an air domain, as shown in fig. 7, wherein the rapid combustion-baking calculation model without the air domain only comprises a test bomb model 2; the outer wall surface 7 of the cartridge case model 5 of the test cartridge model 2 adopts a heat flux boundary as a heating boundary; according to the fast burning calculation model without the air domain, under the condition that the diameter of the charge model 6 in the step 4 is the same, calculating the ignition time B of the charge model 6, and adjusting the size of a heat flux boundary to enable the ignition time A in the step 4 to be equal to the ignition time B; finally, the relationship between the heat flux boundary and the diameter D of the corresponding charge model 6 is obtained, as shown in the scatter diagram in fig. 8;
step 6, fitting the scatter diagram in the graph 8 by adopting a single-phase exponential decay function or Matlab software in Origin software to obtain a relational expression between the heat flux boundary phi and the diameter D of the charging model 6, wherein the relational expression is as follows:
Figure BDA0003591616450000061
and 7, loading the heat flux boundary phi to the outer wall surface 7 of the test bomb model 2 for calculation according to any diameter D of the explosive charging model 6 and the corresponding heat flux boundary phi obtained by calculation in the formula (1), acquiring information such as internal temperature change, ignition time and ignition temperature of the explosive charging model 6, and evaluating the rapid ignition characteristic of the explosive.
In summary, the above description is only a preferred embodiment of the present invention, and is not intended to limit the scope of the present invention. Any modification, equivalent replacement, or improvement made within the spirit and principle of the present invention should be included in the protection scope of the present invention.

Claims (8)

1.一种炸药快速烤燃的评估方法，其特征在于，该评估方法的具体步骤如下：1. an evaluation method for rapid roasting and burning of explosive, is characterized in that, the concrete steps of this evaluation method are as follows: 步骤一，搭建快速烤燃试验系统，开展快速烤燃试验；Step 1, build a rapid burning test system and carry out rapid burning test; 步骤二，根据快速烤燃试验系统建立含空气域的快速烤燃计算模型进行仿真；所述含空气域的快速烤燃计算模型包括：试验弹模型(2)、空气域(4)、燃料入口(1)及燃料出口(3)；Step 2: According to the rapid roasting test system, establish a rapid roasting calculation model containing an air domain for simulation; the rapid roasting calculation model containing an air domain includes: a test bomb model (2), an air domain (4), a fuel inlet (1) and fuel outlet (3); 步骤三，对含空气域的快速烤燃计算模型中的试验弹模型(2)内的装药模型(6)进行不同尺寸的计算，即通过改变装药模型(6)的直径，通过计算获得不同装药模型(6)直径下的装药模型(6)的点火时间A；In step 3, different sizes are calculated for the charge model (6) in the test bomb model (2) in the rapid roasting calculation model with air domain, that is, by changing the diameter of the charge model (6), the calculation is obtained. Ignition time A of the charge model (6) under different charge model (6) diameters; 步骤四，建立不含空气域的快速烤燃计算模型，所述不含空气域的快速烤燃计算模型仅包括试验弹模型(2)；试验弹模型(2)的弹壳模型(5)的外壁面(7)采用热通量边界作为加热边界；根据该不含空气域的快速烤燃计算模型，在与步骤三中的装药模型(6)的直径相同的条件下，计算装药模型(6)的点火时间B，并通过调节热通量边界大小，使得步骤三中的点火时间A与所述点火时间B相等；最终得到热通量边界与对应的装药模型(6)的直径之间的散点对应关系，并将该散点拟合为曲线；Step 4: Establish a rapid roasting calculation model without air domain, and the rapid roasting calculation model without air domain only includes the test bomb model (2); the outer shell of the shell model (5) of the test bomb model (2) The wall surface (7) adopts the heat flux boundary as the heating boundary; according to the fast roasting calculation model without the air domain, under the same conditions as the diameter of the charge model (6) in step 3, the charge model ( 6), and by adjusting the size of the heat flux boundary, the ignition time A in step 3 is equal to the ignition time B; finally, the difference between the heat flux boundary and the diameter of the corresponding charge model (6) is obtained. The corresponding relationship between the scatter points, and the scatter points are fitted as a curve; 步骤五，根据装药模型(6)的任一直径及根据步骤四的曲线计算得到的热通量边界，将该热通量边界加载到试验弹模型(2)的外壁面(7)进行计算，获取装药模型(6)内部温度变化、点火时间及点火温度，对炸药快速烤燃特性进行评估。Step 5, according to any diameter of the charge model (6) and the heat flux boundary calculated according to the curve of step 4, load the heat flux boundary onto the outer wall surface (7) of the test bomb model (2) for calculation. , obtain the internal temperature change, ignition time and ignition temperature of the charge model (6), and evaluate the rapid burn-off characteristics of the explosive. 2.如权利要求1所述的一种炸药快速烤燃的评估方法，其特征在于，在步骤一中，所述试验系统包括：燃油池(8)、试验弹(14)、热电偶(9)、支架(10)、悬挂杆(11)、点火装置(12)及温度记录仪(13)；所述试验弹(14)包括：弹壳(15)和装填在弹壳(15)内的装药(16)；试验弹(14)通过悬挂杆(11)固定在支架(10)上，燃油池(8)位于试验弹(14)的正下方，且试验弹(14)的最低处与燃油池(8)的燃油上表面之间的距离为300mm，燃油池(8)内的燃油与点火装置(12)连接；热电偶(9)安装在试验弹(14)的外部，且热电偶(9)位于所述试验弹(14)的轴线延长线上，所述热电偶(9)与试验弹(14)的外端面的最近距离为50mm，所述温度记录仪(13)与热电偶(9)电性连接；2. the evaluation method of a kind of explosive quick roasting burning as claimed in claim 1, is characterized in that, in step 1, described test system comprises: fuel pool (8), test bomb (14), thermocouple (9) ), a bracket (10), a suspension rod (11), an ignition device (12) and a temperature recorder (13); the test bomb (14) includes: a shell (15) and a charge loaded in the shell (15) (16); the test bomb (14) is fixed on the bracket (10) through the suspension rod (11), the fuel pool (8) is located directly under the test bomb (14), and the lowest part of the test bomb (14) is connected to the fuel pool The distance between the upper surfaces of the fuel (8) is 300mm, and the fuel in the fuel pool (8) is connected to the ignition device (12); the thermocouple (9) is installed outside the test bomb (14), and the thermocouple (9) ) is located on the extension line of the axis of the test bomb (14), the closest distance between the thermocouple (9) and the outer end face of the test bomb (14) is 50mm, and the temperature recorder (13) is connected to the thermocouple (9). ) electrical connection; 进行快速烤燃试验时，点火装置(12)点燃所述燃油，燃油池(8)内的燃油燃烧产生火焰及高温燃气，并在试验弹(14)的周围形成火焰场，试验弹(14)的装药(16)在火焰场内升温，当火焰场的温度高过试验弹(14)的装药(16)的燃点时当试验弹(14)的装药(16)的最高温度高于其燃点时，装药(16)点火；热电偶(9)测量在燃油池(8)的燃油燃烧时，试验弹(14)所处火焰场的环境温度，温度记录仪(13)记录热电偶(9)所测量的环境温度随时间的变化，得到试验的环境温度随时间的变化曲线。During the rapid roasting test, the ignition device (12) ignites the fuel oil, the fuel oil in the fuel pool (8) burns to generate flame and high-temperature gas, and a flame field is formed around the test bomb (14), and the test bomb (14) When the temperature of the flame field is higher than the ignition point of the charge (16) of the test bomb (14), the maximum temperature of the charge (16) of the test bomb (14) is higher than When it ignites, the charge (16) ignites; the thermocouple (9) measures the ambient temperature of the flame field where the test bomb (14) is located when the fuel in the fuel pool (8) burns, and the temperature recorder (13) records the thermocouple (9) Variation of the measured ambient temperature with time, and obtain the variation curve of the ambient temperature of the test with time. 3.如权利要求2所述的一种炸药快速烤燃的评估方法，其特征在于，在步骤二中，所述试验弹模型(2)与所述试验弹(14)完全相同，包括：弹壳模型(5)和装填在弹壳模型(5)内的装药模型(6)；所述燃料入口(1)为燃油池(8)的燃油燃烧产生的高温燃气，燃料入口(1)的边界条件为质量流量入口边界，质量流量入口边界采用温度分段函数和质量流率拟合试验环境温度；试验弹模型(2)的最低处距燃料入口(1)距离与试验弹(14)的最低处与燃油池(8)的燃油上表面之间的距离相同；燃料出口(3)的边界条件为压力出口边界；通过调节所述燃气质量流率和温度分段函数来模拟燃油池(8)的燃油燃烧产生的火焰场的情况，计算得到试验弹模型(2)所处火焰场的环境温度随时间的变化曲线，若该计算得到试验弹模型(2)所处火焰场的环境温度随时间的变化曲线与试验得到的环境温度随时间的变化曲线在同一时间点的环境温度差值在设定范围内，说明所建立的计算模型能够描述所述快速烤燃试验系统的火焰场。3. the evaluation method of a kind of explosive quick roasting burning as claimed in claim 2, is characterized in that, in step 2, described test bomb model (2) is identical with described test bomb (14), comprises: shell case The model (5) and the charge model (6) loaded in the cartridge case model (5); the fuel inlet (1) is the high-temperature gas generated by the combustion of fuel in the fuel pool (8), and the boundary condition of the fuel inlet (1) is the mass flow inlet boundary, and the mass flow inlet boundary adopts the temperature piecewise function and mass flow rate to fit the test ambient temperature; the distance between the lowest point of the test bomb model (2) and the fuel inlet (1) and the lowest point of the test bomb (14) The distance between the upper surface of the fuel oil of the fuel pool (8) is the same; the boundary condition of the fuel outlet (3) is the pressure outlet boundary; the fuel pool (8) is simulated by adjusting the gas mass flow rate and the temperature piecewise function. In the case of the flame field generated by the combustion of fuel oil, the variation curve of the ambient temperature of the flame field where the test bomb model (2) is located with time can be calculated. The ambient temperature difference between the variation curve and the ambient temperature variation curve obtained by the experiment at the same time point is within the set range, indicating that the established calculation model can describe the flame field of the rapid roasting test system. 4.如权利要求1-3任一项所述的一种炸药快速烤燃的评估方法，其特征在于，在步骤四中，采用Origin软件中的单相指数衰减函数，对热通量边界与对应的装药模型(6)的直径之间的散点对应关系进行拟合。4. the evaluation method of a kind of explosive quick roasting burning as described in any one of claim 1-3, it is characterized in that, in step 4, adopt the single-phase exponential decay function in Origin software, to heat flux boundary and The scatter correspondence between the diameters of the corresponding charge model (6) is fitted. 5.如权利要求1-3任一项所述的一种炸药快速烤燃的评估方法，其特征在于，在步骤四中，采用Matlab软件对热通量边界与对应的装药模型(6)的直径之间的散点对应关系进行拟合。5. the evaluation method of a kind of explosive quick roasting burning as described in any one of claim 1-3, it is characterised in that in step 4, adopt Matlab software to heat flux boundary and corresponding charge model (6) The scatter correspondence between the diameters of , is fitted. 6.如权利要求4所述的一种炸药快速烤燃的评估方法，其特征在于，在步骤四中，热通量边界与对应的装药模型(6)的直径之间的散点对应关系拟合为曲线后的表达式如下：6. the evaluation method of a kind of explosive fast roasting burning as claimed in claim 4, is characterized in that, in step 4, the scatter point correspondence between the diameter of heat flux boundary and corresponding charge model (6) The expression after fitting to a curve is as follows:
Figure FDA0003591616440000021
Figure FDA0003591616440000021
 其中，Φ为热通量边界，D为装药模型(6)的直径。Among them, Φ is the heat flux boundary, and D is the diameter of the charge model (6). 7.如权利要求5所述的一种炸药快速烤燃的评估方法，其特征在于，在步骤四中，热通量边界与对应的装药模型(6)的直径之间的散点对应关系拟合为曲线后的表达式如下：7. the evaluation method of a kind of explosive fast roasting burning as claimed in claim 5, is characterized in that, in step 4, the scatter point correspondence between the diameter of heat flux boundary and corresponding charge model (6) The expression after fitting to a curve is as follows:
Figure FDA0003591616440000022
Figure FDA0003591616440000022
 其中，Φ为热通量边界，D为装药模型(6)的直径。Among them, Φ is the heat flux boundary, and D is the diameter of the charge model (6). 8.如权利要求1-3任一项所述的一种炸药快速烤燃的评估方法，其特征在于，在步骤三之前，需要对步骤二建立的含空气域的快速烤燃计算模型进行网格划分，划分空气域(4)网格尺寸为100mm，划分试验弹模型(2)的网格尺寸为2mm。8. the evaluation method of a kind of explosive quick roasting burning as described in any one of claim 1-3, it is characterized in that, before step 3, need to carry out network to the rapid roasting burning calculation model of the air-containing domain established in step 2. The grid size of the divided air domain (4) is 100 mm, and the grid size of the divided test bomb model (2) is 2 mm.
CN202210379372.4A 2022-04-12 2022-04-12 A method for evaluating rapid cooking of explosives Active CN114813835B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202210379372.4A CN114813835B (en) 2022-04-12 2022-04-12 A method for evaluating rapid cooking of explosives

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202210379372.4A CN114813835B (en) 2022-04-12 2022-04-12 A method for evaluating rapid cooking of explosives

Publications (2)

Publication Number Publication Date
CN114813835A true CN114813835A (en) 2022-07-29
CN114813835B CN114813835B (en) 2025-06-24

Family

ID=82534372

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202210379372.4A Active CN114813835B (en) 2022-04-12 2022-04-12 A method for evaluating rapid cooking of explosives

Country Status (1)

Country Link
CN (1) CN114813835B (en)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN121093637A (en) * 2025-11-06 2025-12-09 江汉大学 A method for predicting the detonation velocity decay of emulsion explosives during thermal aging
CN121093637B (en) * 2025-11-06 2026-02-10 江汉大学 A method for predicting the detonation velocity decay of emulsion explosives during thermal aging

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4622922A (en) * 1984-06-11 1986-11-18 Hitachi, Ltd. Combustion control method
CN106895750A (en) * 2017-02-20 2017-06-27 北京理工大学 A kind of ammunition fast cook off test experimental system
CN107167490A (en) * 2017-05-12 2017-09-15 北京理工大学 A kind of fast cook off test experimental rig for medium-and-large-sized ammunition
CN111912303A (en) * 2020-07-22 2020-11-10 北京理工大学 Thermal stimulation ammunition charging ignition and reaction growth overall process testing system
CN112818575A (en) * 2021-01-26 2021-05-18 北京理工大学 Macro-micro numerical analysis method for explosive burning

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4622922A (en) * 1984-06-11 1986-11-18 Hitachi, Ltd. Combustion control method
CN106895750A (en) * 2017-02-20 2017-06-27 北京理工大学 A kind of ammunition fast cook off test experimental system
CN107167490A (en) * 2017-05-12 2017-09-15 北京理工大学 A kind of fast cook off test experimental rig for medium-and-large-sized ammunition
CN111912303A (en) * 2020-07-22 2020-11-10 北京理工大学 Thermal stimulation ammunition charging ignition and reaction growth overall process testing system
CN112818575A (en) * 2021-01-26 2021-05-18 北京理工大学 Macro-micro numerical analysis method for explosive burning

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
李文凤;余永刚;叶锐;: "底排药快速烤燃特性的数值模拟", 含能材料, no. 10, 31 December 2016 (2016-12-31) *

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN121093637A (en) * 2025-11-06 2025-12-09 江汉大学 A method for predicting the detonation velocity decay of emulsion explosives during thermal aging
CN121093637B (en) * 2025-11-06 2026-02-10 江汉大学 A method for predicting the detonation velocity decay of emulsion explosives during thermal aging

Also Published As

Publication number Publication date
CN114813835B (en) 2025-06-24

Similar Documents

Publication Publication Date Title
CN114813834B (en) A method for evaluating the slow cooking of explosives
Togashi et al. Numerical simulation of long-duration blast wave evolution in confined facilities
Li et al. Effects of Charge Size on Slow Cook‐Off Characteristics of AP/HTPB Composite Propellant in Base Bleed Unit
CN109100147A (en) A kind of experimental provision for assessing solid engines slow cook-off responsiveness
CN114813835A (en) Evaluation method for rapid burning of explosive
CN117592389A (en) Gas temperature-pressure change rule calculation method for semi-closed exploder
CN108801084B (en) Projectile body exhaust slow-release environment simulation test device
Zhang et al. Study of in-cylinder heat transfer boundary conditions for diesel engines under variable altitudes based on the CHT model
CN119378434B (en) A modeling method for energy conversion and work done by a projectile in a launch path
CN108562507B (en) Experimental method and device for anti-ablation wear of rapid-fire weapon barrel
Bougamra et al. Multiphase CFD simulation of solid propellant combustion in a small gun chamber
CN112666216A (en) Heat flow dust cloud ignition temperature characteristic testing device and testing method thereof
CN109632613B (en) A method for estimating safe service life of press-packed penetrating explosive JHLDu-1
Wei et al. Simulation study on the impact response of barrels with different rifling profiles during bullet engraving
Wei et al. Mechanical–thermal–chemical coupled modeling techniques and numerical simulation methods for the firearm igniting process
CN211856417U (en) Multi-temperature-zone hot baking test device for energetic material product
Chen et al. Numerical simulation and optimized design of a launching system with a novel charge structure
CN118070586A (en) Method for estimating slow baking combustion reaction intensity of full-size explosive charge
He et al. Compression fracture, heat transfer, and uncertain optimization of launch propellant charge
CN111046531B (en) Baking combustion numerical simulation calculation method for HTPE type solid engine
Moussa et al. An Approach to Predict the Slow Cook-off Response of Confined & Vented Munitions Based on Small Scale Tests.
CN117350048B (en) Energy safety threshold assessment method for fuel system based on transient flame expansion process
CN113969758B (en) Oil pipe tubular column cutting device with pressure balancer
Cooper Transient temperature distribution and thermo-elastic stress in gun tubes
CN110987456A (en) Engine simulation device for ignition device test

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
GR01 Patent grant
GR01 Patent grant