CN113553739A - Method for calculating explosion output characteristics of mixed explosive - Google Patents

Method for calculating explosion output characteristics of mixed explosive Download PDF

Info

Publication number
CN113553739A
CN113553739A CN202110761644.2A CN202110761644A CN113553739A CN 113553739 A CN113553739 A CN 113553739A CN 202110761644 A CN202110761644 A CN 202110761644A CN 113553739 A CN113553739 A CN 113553739A
Authority
CN
China
Prior art keywords
calculation
equation
chemical reaction
formula
constructing
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
CN202110761644.2A
Other languages
Chinese (zh)
Other versions
CN113553739B (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.)
Xian Modern Chemistry Research Institute
Original Assignee
Xian Modern Chemistry Research Institute
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 Xian Modern Chemistry Research Institute filed Critical Xian Modern Chemistry Research Institute
Priority to CN202110761644.2A priority Critical patent/CN113553739B/en
Publication of CN113553739A publication Critical patent/CN113553739A/en
Application granted granted Critical
Publication of CN113553739B publication Critical patent/CN113553739B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/20—Design optimisation, verification or simulation
    • G06F30/28—Design optimisation, verification or simulation using fluid dynamics, e.g. using Navier-Stokes equations or computational fluid dynamics [CFD]
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/20—Design optimisation, verification or simulation
    • G06F30/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
    • 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
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2113/00—Details relating to the application field
    • G06F2113/08—Fluids
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/14—Force analysis or force optimisation, e.g. static or dynamic forces

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Computing Systems (AREA)
  • General Engineering & Computer Science (AREA)
  • Geometry (AREA)
  • Evolutionary Computation (AREA)
  • Chemical & Material Sciences (AREA)
  • Computer Hardware Design (AREA)
  • Bioinformatics & Cheminformatics (AREA)
  • Crystallography & Structural Chemistry (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Chemical Kinetics & Catalysis (AREA)
  • Bioinformatics & Computational Biology (AREA)
  • Analytical Chemistry (AREA)
  • Algebra (AREA)
  • Fluid Mechanics (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Mathematical Physics (AREA)
  • Pure & Applied Mathematics (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

The invention provides a method for calculating explosion output characteristics of a mixed explosive, which is used for calculating the explosion output characteristics of the mixed explosive based on a fluid dynamics equation under a two-dimensional Euler coordinate system, and comprises the following calculation steps: step 100, constructing a calculation model; step 200, constructing a state equation; step 300, constructing an explosive chemical reaction model; step 400, constructing an elastic-plastic model; step 500, calculating mass migration; step 600, reassign flow field variables. The calculation method comprises a plurality of chemical reaction models, can effectively describe the detonation reaction process of the mixed explosive, is suitable for the energy output characteristic calculation problem of a plurality of different processes, and has wide applicability. The calculation method has reliable calculation accuracy. The invention can conveniently and repeatedly calculate the number of the multiple types and times, is beneficial to improving the working efficiency and reducing the working intensity of personnel.

Description

Method for calculating explosion output characteristics of mixed explosive
Technical Field
The invention belongs to the field of mixed explosives, relates to numerical calculation of detonation energy release rules of mixed explosives, and particularly relates to a method for calculating explosion output characteristics of mixed explosives.
Background
The evaluation of the detonation energy release characteristics of the mixed explosive is an important aspect for measuring the power of the mixed explosive. On one hand, researchers in related fields need to carry out corresponding design work according to actual working conditions; on the other hand, the user also needs to make reasonable choices according to the energy characteristics of the mixed explosive.
The method is the most direct evaluation method for testing the explosion output characteristic parameters of the mixed explosives by adopting various experimental techniques. However, as the mixed explosive belongs to dangerous chemicals, no small potential risk is involved in sample preparation and experimental measurement. And the economic cost is higher when carrying out the experiment measurement often, still has the big problem of test data acquisition degree of difficulty.
In order to reduce the dependence on experiments and improve the design level, the explosion energy release characteristics of the mixed explosive need to be designed and calculated. The detonation energy action process of the mixed explosive can be physically regarded as a flow process changing along with time, and the detonation reaction of the mixed explosive is described by combining various chemical reaction models. In the basic principle framework of fluid dynamics, the change rule of various physical quantities in an explosion flow field can be obtained through a numerical calculation method, so that the general process and characteristic rule of the energy output of the mixed explosive are reflected.
The Chinese patent with the patent number ZL201610034329.9 discloses a high-energy high-frequency explosive source design method, which comprises the step of calculating by using finite element models for explosion of different explosives. However, the invention only simplifies the explosive explosion process to process a seismic source process, and does not perform detailed calculation on the energy output characteristics of the mixed explosive.
Therefore, it is necessary to invent a method for calculating the detonation output characteristics of the mixed explosive, which is based on the general principle of fluid dynamics and calculates the detonation energy characteristic release and work output rules of the mixed explosive.
Disclosure of Invention
Aiming at the defects in the prior art, the invention aims to provide a method for calculating the explosion output characteristics of a mixed explosive, and solve the technical problem that the mixed explosive explosion flow field containing a chemical reaction process under Euler coordinates cannot be calculated in the prior art.
In order to solve the technical problems, the invention adopts the following technical scheme:
a method for calculating explosion output characteristics of a mixed explosive is based on a fluid dynamic equation under a two-dimensional Euler coordinate system and is used for calculating the explosion output characteristics of the mixed explosive, and comprises the following calculation steps:
step 100, constructing a calculation model:
constructing a mass conservation equation as shown in a formula 1.1 under a two-dimensional Euler coordinate system, a momentum conservation equation as shown in a formula 1.2 and a formula 1.3, and an energy conservation equation as shown in a formula 1.4;
conservation of mass equation:
Figure BDA0003150095330000021
conservation of momentum in the x direction equation:
Figure BDA0003150095330000022
conservation of momentum in the y-direction equation:
Figure BDA0003150095330000023
energy conservation equation:
Figure BDA0003150095330000031
formula 1.1-formula 1.4:
rho is density;
u is the speed;
e is energy;
s is stress offset;
p is the pressure intensity;
t is time;
alpha is a symmetric coefficient, alpha is 2 under a cylindrical coordinate system, and alpha is 1 under a rectangular coordinate system;
subscripts x and y are the components of the respective directions;
step 200, constructing a state equation:
the physical property change rule of the detonation product of the mixed explosive is described by adopting an HOM state equation, and the expression form is shown as formula 1.5:
Figure BDA0003150095330000032
in formula 1.5:
pCJthe pressure of detonation products in the CJ state;
eCJinternal energy of detonation product under CJ state;
TCJis the temperature of detonation products under CJ conditions;
subscript CJ indicates CJ status as a reference;
Vgis the volume of detonation products;
coefficients A, B, C, D and E, K, L, M, N and O, and Q, R, S, U and X are fitting coefficients of the respective physical quantities, respectively;
step 300, constructing an explosive chemical reaction model:
the method comprises an Arrhenius chemical reaction model, a Taylor wave chemical reaction model, a square wave impact chemical reaction model, a CJ volume combustion chemical reaction model and a Forest-Fire heterogeneous combustion chemical reaction model;
step 400, constructing an elastic-plastic model:
when the yield strength of the material is exceeded, an elastic-plastic model is used to describe the stress state of the object, expressed in the form of equation 1.6:
Figure BDA0003150095330000041
in formula 1.6:
f represents the principal stress;
Y0indicates the yield strength of the material;
step 500, mass migration is calculated:
transferring the mass according to the motion speed of the flow field, wherein the transferred mass flows into a new computational grid;
step 600, redistributing flow field variables:
flow field variables including density, mass fraction, x-direction velocity, y-direction velocity, and internal energy are redistributed according to the mass transfer calculation in step 500.
The invention also has the following technical characteristics:
in step 300, the Forest-Fire heterogeneous combustion chemical reaction model further comprises a linear volume combustion chemical reaction model.
Step 300 also includes adding a viscosity coefficient.
Compared with the prior art, the invention has the following technical effects:
the calculation method comprises a plurality of chemical reaction models, can effectively describe the detonation reaction process of the mixed explosive, is suitable for the energy output characteristic calculation problem of a plurality of different processes, and has wide applicability.
(II) the calculation method has reliable calculation precision.
(III) by applying the invention, the calculation of various types and multiple times can be conveniently and repeatedly carried out, the working efficiency is improved, and the working intensity of personnel is reduced.
Drawings
FIG. 1 is a flow chart of a method for calculating the explosion output characteristics of the mixed explosive.
FIG. 2 is a cloud of calculated pressure results of mixed explosive EX1 in example 1 when it was applied to a steel cylinder in an air medium.
FIG. 3 is a cloud of calculated pressure results of mixed explosive EX2 detonation in an aqueous medium in example 2.
FIG. 4 is a cloud chart of the calculated pressure during the process of converting combustion to detonation of the mixed explosive EX3 in example 3.
The present invention will be explained in further detail with reference to examples.
Detailed Description
In the present invention, it is to be noted that:
all physical quantities are in units commonly used in the art.
CJ refers to the detonation wave CJ theory, namely the planar one-dimensional fluid dynamics theory about detonation waves proposed by Chapman and Jouget respectively in the beginning of the 20 th century.
The HOM equation of state is a thermodynamic equation of state that describes the physical state of the unreacted mixed explosive and the detonation products of the mixed explosive.
Aiming at the defects in the prior art, the invention discloses a method for calculating the explosion output characteristics of a mixed explosive, which is based on the general principle of fluid dynamics and adopts a conservation equation under a two-dimensional Euler coordinate to calculate the explosion reaction flow field of the mixed explosive containing a chemical reaction process, thereby reflecting the main rule of the energy output characteristics of the explosion process of the mixed explosive.
From the analysis of the calculation principle, the method for calculating the explosion output characteristics of the mixed explosive disclosed by the invention solves the problem of calculating the explosion flow field of the mixed explosive containing the chemical reaction process under the Euler coordinate, and can be applied to numerical calculation under a two-dimensional rectangular coordinate system or a two-dimensional cylindrical coordinate system.
For each computation time step, it can be generalized to three computation segments.
The first calculation step is that according to the state equation of the explosion product of the mixed explosive, the chemical reaction process is calculated to obtain the pressure change and the stress state change in the explosion flow field;
the second calculation step, according to the pressure and stress change condition obtained in the first calculation step, driving the flow field to generate migration, temporarily neglecting the influence of a migration item, and calculating and processing the flow field speed and energy by combining the conservation relation of fluid dynamics;
and a third calculation step, wherein the mass transfer process is omitted in the second calculation step, the mass transfer is calculated by using a donor-acceptor method, the mass is transferred to a new part, and then each physical quantity in the flow field is redistributed.
And completing the calculation of the current time step by the sequential connection of the three calculation links. And performing reciprocating circulation on the three calculation links so as to complete the calculation of the whole flow field change under the condition of predicted time. And finally obtaining a calculation result of the output characteristic of the mixed explosive.
The three calculation links can be further refined into six calculation stages:
the first stage is as follows: the calculation of pressure and stress changes includes:
(a) calculating a state equation of a detonation product of the mixed explosive;
(b) calculating the chemical reaction process of the mixed explosive;
(c) and calculating the offset stress.
And a second stage: calculation of elastoplasticity comprising:
(d) calculating the flow field trial speed;
(e) and (4) calculating the elasticity and the plasticity.
And a third stage: calculation of energy, comprising:
(f) calculating the internal energy of detonation products;
(g) and (4) calculating the total energy of the flow field.
A fourth stage: a calculation of mass transfer comprising:
(h) and calculating mass transfer.
The fifth stage: redistribution of flow field variables, comprising:
(i) redistribution of flow field variables.
The sixth stage: elastoplastic flow comprising:
(j) calculating the elastoplastic flow of the flow field;
(k) various forms of construction of partial derivatives.
From the aspect of execution operation analysis, when the mixed explosive explosion output characteristic calculation is specifically executed, three specific steps of calculation parameter file construction, solving calculation and calculation result visualization processing are included.
The calculation parameter file construction is combined with a specific practical problem, and is used for constructing a model parameter file which is used for solving calculation and can be identified, wherein the model parameter file comprises material property parameters, state equation parameters, calculation grid parameters, initial condition parameters, reaction model parameters, calculation control parameters and result output parameters.
The solving calculation is a numerical calculation process of the specified problem and is used for obtaining all calculation results of the mixed explosive detonation reaction flow field.
The calculation result visualization processing is to perform graphical processing on the specified calculation result data by combining with the user requirements to form a distribution cloud picture and a time course curve.
Specifically, the calculation comprises the following calculation steps:
step 100, constructing a calculation model, and constructing a fluid mechanics mass, momentum and energy conservation equation under a two-dimensional Euler coordinate system, wherein the equations are respectively shown as formula 1, formula 2, formula 3 and formula 4:
conservation of mass equation:
Figure BDA0003150095330000081
conservation of momentum in the x direction equation:
Figure BDA0003150095330000082
conservation of momentum in the y-direction equation:
Figure BDA0003150095330000083
energy conservation equation:
Figure BDA0003150095330000084
in equations 1 to 4, ρ is density, U is velocity, e is energy, S is stress offset, p is pressure, t is time, α is a symmetry coefficient, α is 2 in a cylindrical coordinate system, α is 1 in a rectangular coordinate system, and subscripts x and y are components in corresponding directions.
Under the Euler coordinate system, the conservation of mass equation is satisfied automatically.
Step 200, constructing a state equation, describing the physical property change of the detonation product of the mixed explosive by adopting a gaseous HOM state equation, wherein the form is shown as a formula 5:
Figure BDA0003150095330000085
in formula 5, pCJIs the pressure of detonation products in CJ conditions, eCJIs the internal energy of detonation product in CJ state, TCJThe temperature of the detonation product in CJ condition, the subscript CJ denoting the CJ condition as a reference, VgThe coefficients A, B, C, D and E, K, L, M, N and O and Q, R, S, U and X are fitting coefficients for the respective physical quantities, which are the volumes of detonation products.
For the third expression in the expression 5, the partial derivative is calculated from the polynomial function relationship of the temperature and the volume, and the Gruneisen coefficient expression of the gas phase detonation product can be obtained, as shown in the expression 6:
Figure BDA0003150095330000091
the state change of other solids in the flow field is described by adopting a solid state form HOM state equation, and the form is shown as formula 7:
Figure BDA0003150095330000092
in formula 6, THTemperature under the Hugoniot impact curve, pHIs the pressure under the Hugoniot impact curve, eHIs the internal energy of the Hugoniot impact curve, VsIs the solid volume, ΓsGruneisen coefficient for solids, CVFor the isochoric heat capacity of the solid, the subscript H indicates the coefficients F, G, H, I and J are polynomial function fit coefficients of temperature versus solid volume, referenced to the Hugoniot impact regime.
If expression 8 is taken as a reference for the relationship of the Hugoniot impact curve of a solid:
Us＝c+supformula 8
Then the expression for the impact pressure and the impact energy of the solid can be derived, as shown in equation 9:
Figure BDA0003150095330000093
in formulas 8 and 9, UsIs the velocity of the shock wave upIs the particle velocity, c and s are first order fitting coefficients, V0Initial volume of solid, pH、eHAnd VsHas the same meaning as formula 7.
For a mixed explosive, the HOM equation of state parameters in the unreacted state, i.e. in the solid, and the detonation products, i.e. in the gas, need to be given separately.
For other solids or liquids in the flow field, only the HOM equation of state parameters of the solid form need be given.
And 300, constructing an explosive chemical reaction model, wherein the explosive chemical reaction model comprises an Arrhenius chemical reaction model, a Taylor wave chemical reaction model, a square wave impact chemical reaction model, a CJ volume combustion chemical reaction model and a Forest-Fire heterogeneous combustion chemical reaction model.
The expression of the Arrhenius chemical reaction model is shown in formula 10:
Figure BDA0003150095330000101
in formula 10, WiRepresenting the mass fraction of the mixed explosive in the ith calculation grid, delta t is the calculation time step length, Z is a pre-index factor, EaFor activation energy, R is the gas universal constant, TiFor the temperature in the ith computational grid, superscripts n and n +1 represent the computation time steps.
The expression of the Taylor wave chemical reaction model is shown in formula 11:
Figure BDA0003150095330000102
in formula 11, UCJ、DCJ、pCJAnd VCJRespectively representing the particle velocity, shock wave velocity, pressure and volume under CJ state, gamma is the adiabatic index of detonation products, rho0Is an initial density, V0Is the initial volume, p0V0＝1。
The expression of the square wave impact chemical reaction model is shown as formula 12:
Figure BDA0003150095330000111
the symbols in formula 12 have the same meanings as in formula 11.
The expression of the CJ volume combustion chemical reaction model is shown in formula 13:
Figure BDA0003150095330000112
the meanings of the symbols in formula 13 are shown in formula 10 and formula 11.
The expression of the Forest-Fire heterogeneous combustion chemical reaction model is shown as formula 14:
Figure BDA0003150095330000113
in formula 14, WiΔ t and the superscripts n and n +1 have the same meaning as in equation 10, the factor a0、a1、a2.., and amFitting coefficients to polynomials of powers of various orders with respect to the pressure p.
Further, the Forest-Fire heterogeneous fuel chemical reaction model also comprises a linear volume combustion chemical reaction model.
The expression of the linear volume combustion chemical reaction model is shown in formula 15:
Figure BDA0003150095330000114
in formula 15:
w is the mass fraction of the mixed explosive, (S)0/V0) Represents a combustion usable specific surface area, q represents a shape factor, q-0 represents a planar constant combustion, q-1/2 represents a cylindrical pellet combustion, q-2/3 represents a spherical pellet combustion; c represents a combustion linearity coefficient, and b represents a combustion pressure index. The minus sign on the right of the equal sign indicates that the mass fraction decreases as the chemical reaction proceeds.
Further, in step 300, a step of adding a viscosity coefficient is included, wherein the viscosity coefficient includes artificial viscosity or real viscosity, so that the calculation result is stable.
Step 400, constructing an elastic-plastic model, and when the yield strength of the material is exceeded, describing the stress state of the object by using the elastic-plastic model, wherein the expression form is shown as formula 16:
Figure BDA0003150095330000121
in formula 16, f represents the principal stress; y is0Indicating the yield strength of the material, S is the stress offset, and subscripts x and y are the components in the respective directions.
And 500, calculating mass migration, wherein the mass is migrated according to the motion speed of the flow field, and the migrated mass flows into a new calculation grid.
Mass migration was calculated using the donor-acceptor method. Mass migration in the donor and acceptor grids was calculated from the mass fraction of the material in the acceptor grid. The density and energy of the components corresponding to the migrated mass in the donor mesh are calculated from the equation of state in the first stage. Only one material reduction is allowed in one calculation time step.
For j migration from donor grid (i-1) to acceptor grid (i, j) along x-direction, in the case of two-dimensional cylindrical symmetry coordinates, if migration terms of energy, mass fraction, x-direction momentum, y-direction momentum, mass and density in the grid are expressed in DE, DW, DPU, DPV, DM and DMASS, respectively, their expressions are shown in equations 17-22:
Figure BDA0003150095330000122
Figure BDA0003150095330000123
Figure BDA0003150095330000124
Figure BDA0003150095330000125
DMi,j(DMASS) (GEO) formula 21
Figure BDA0003150095330000131
In equations 17 to 22, subscripts i, j denote calculation grid numbers, i denotes the x direction, j denotes the y direction, superscripts n denotes the calculation time step, UxSpeed in x direction, UyIs the y-direction velocity and ρ is the density.
In equations 17 to 22, GEO represents a shape factor, and the expression is shown in equation 23:
Figure BDA0003150095330000132
in equation 23, i represents the x-direction calculation grid number, and the number is counted in order from i to 1.
In equation 22, the discriminant is expressed by Δ, and the expression is shown in equation 24:
Figure BDA0003150095330000133
in equation 24, subscripts i, j denote calculation grid numbers, i denotes the x direction, j denotes the y direction, superscripts n denotes the calculation time step, UxThe speed in the x direction, Δ t is the calculation time step, and Δ x is the calculation space step in the x direction.
Step 600, redistributing the flow field variables, and redistributing the flow field variables including density, mass fraction, x-direction velocity, y-direction velocity and internal energy according to the mass migration calculation result in step 500.
According to the mass transfer calculation result, corresponding mass transfer terms are added to corresponding physical quantities, and expressions of density, mass fraction, x-direction speed, y-direction speed and internal energy are respectively shown as formulas 25 to 29:
Figure BDA0003150095330000134
Figure BDA0003150095330000135
Figure BDA0003150095330000136
Figure BDA0003150095330000137
Figure BDA0003150095330000138
in formulae 25 to 29, each symbol has the same meaning as in formulae 17 to 22.
And after redistribution, the flow field is continuously enabled to flow, the calculation of the current calculation step is completed, the corresponding calculation result is output and stored, circulation is continuously carried out until the predicted calculation time is reached, and then result visualization processing is carried out. Thereby completing the calculation of the explosion output characteristics of the mixed explosive.
The invention is based on the fluid dynamics principle, and carries out numerical calculation on the mixed explosive explosion reaction flow field containing chemical reaction under the framework of Euler coordinate system, so as to reflect the general process and rule of the change of the explosion output characteristic of the mixed explosive. When the specific technical scheme is involved, on the basis of a fluid mechanics mass, momentum and energy conservation equation, a plurality of chemical reaction models are adopted to describe the reaction process of the mixed explosive, an HOM state equation is used to describe the change rule of physical properties of detonation gas products and other solid materials in a flow field, the elasto-plastic mechanical response of various materials is considered, and a control equation is solved through a numerical method, so that the calculation result of the change of each physical variable in the flow field along with time is obtained. When numerical calculation is carried out, a trial-correction method is used for carrying out loop iteration to ensure that the precision of a calculation result meets the requirement; the calculation process is also stabilized by adjusting the viscosity coefficient. When mass transfer is involved, processing in the donor-acceptor approach reflects the intrinsic nature of mass transfer and redistribution in the flow field in the euler coordinate frame.
The method for calculating the explosion output characteristics of the mixed explosive is based on the general principle of fluid dynamics, carries out numerical calculation on a mixed explosive explosion reaction flow field containing a chemical reaction process under Euler coordinates, reflects the law of the explosion output characteristics of the mixed explosive, and comprises the steps of constructing a calculation model, constructing a state equation, constructing a chemical reaction model, constructing an elasto-plastic model, calculating mass transfer and redistributing flow field variables. The invention adopts various chemical reaction models to describe the explosion reaction process of the mixed explosive and is suitable for the energy output characteristic calculation of the mixed explosive in various different processes.
The following embodiments of the present invention are provided, and it should be noted that the present invention is not limited to the following embodiments, and all equivalent changes based on the technical solutions of the present invention are within the protection scope of the present invention.
Example 1:
according to the technical scheme, the mixed explosive EX1 is used for calculating the steel cylinder explosion effect in the air medium by using a mixed explosive explosion output characteristic calculation method.
The mixed explosive EX1 comprises the following components in parts by weight: hexogen (RDX): polyisobutylene (PIB): dioctyl sebacate (DOP): oil (Oil) ═ 91: 5.3: 2.1: 1.6.
the model contains three materials: mix explosive, air and steel.
Step 100, constructing a calculation model which comprises mass, momentum and energy conservation equations, wherein 50 calculation grids are arranged in the x direction, 100 calculation grids are arranged in the y direction, the total number of the calculation grids is 5000, the initial temperature is 300K, the initial pressure is 0.1MPa, the calculation time step length is 0.05 mu s, and the calculation total time length is 20 mu s;
step 200, constructing a state equation, wherein the mixed explosive is described by adopting a solid phase HOM state equation when not reacting, and the parameters of the state equation are as follows: the coefficients F ═ 0.0000, G ═ -9.0419, H ═ -71.3185, I ═ -125.2050, and J ═ -92.0424;
the detonation reaction product is described by adopting a gas-phase HOM equation of state, and the parameters of the equation of state are as follows: the coefficients a-3.5391, B-2.5774, C-0.2601, D-0.0139 and E-0.0114, the coefficients K-1.6191, L-0.5215, M-0.0678, N-0.0043 and O-1.0468 × 10-4The coefficients Q-7.3642, R-0.4937, S-0.0292, U-0.0330 and X-0.0115;
air is described using the gas phase HOM equation of state, with a-4.5060, B-1.2755, C-0.0037, D-0.0124, and E-0.0021, K-1.6266, L-0.0905, M-0.0027, N-5.4358 × 10-5And O ═ 1.5852 × 10-6Coefficient Q is 8.2264, R ═ -0.2515, S ═ -0.0134, U ═ 0.0141, and X ═ 0.0022;
describing steel by adopting a solid phase HOM state equation; the coefficient F is 0.0000, G is-3.8238X 103，H＝-7.0321×103、I＝-4.8267×103And J ═ 1.4668 × 103；
Step 300, constructing an explosive chemical reaction model, describing the mixed explosive by adopting a CJ volume fuel chemical reaction model, wherein the volume V is in a CJ stateCJ0.407, viscosity coefficient 0.01;
step 400, constructing an elastic-plastic model, wherein the shear modulus of steel is 78GPa, and the yield strength is 690 MPa;
step 500, calculating mass migration;
step 600, reassign flow field variables.
The embodiment can calculate the explosion effect output characteristic of the mixed explosive in the air medium, the calculation precision meets the actual requirement, and the design level is improved.
Example 2:
according to the technical scheme, the explosion output characteristic of the mixed explosive EX2 in the water medium is calculated by using a mixed explosive explosion output characteristic calculation method.
The mixed explosive EX2 comprises the following components in parts by weight: hexogen (RDX): trinitrotoluene (TNT): aluminum powder (Al): wax (Wax) ═ 40: 38: 17: 5.
the model contains two materials: the explosive and water are mixed.
Step 100, constructing a calculation model which comprises mass, momentum and energy conservation equations, wherein the number of calculation grids is 100 in the x direction, the number of calculation grids is 200 in the y direction, the total number is 20000, the initial temperature is 300K, the initial pressure is 0.1MPa, the calculation time step length is 0.05 mu s, and the total calculation time length is 30 mu s;
step 200, constructing a state equation, wherein the mixed explosive is described by adopting a solid phase HOM state equation when not reacting, and the parameters of the state equation are as follows: the coefficients F ═ 0.0000, G ═ -8.6662, H ═ -58.3138, I ═ -69.7160, and J ═ -8.2010;
detonation reaction products are employedDescribing a gas-phase HOM state equation, wherein the parameters of the state equation are as follows: the coefficients a-3.6653, B-2.4671, C-0.2285, D-0.0605 and E-0.0193, the coefficients K-1.1592, L-0.5304, M-0.0956, N-0.0091 and O-0.3444 × 10-3The coefficient Q is 7.5771, R is-0.4388, S is 0.0923, and U is 2.3668 × 10-3And X ═ 3.2434 × 10-2；
Describing water by adopting a solid phase HOM state equation; the coefficients F-0.0000, G-5.7206, H-0.6926, I-8.8139 and J-36.0120;
step 300, constructing an explosive chemical reaction model, describing the mixed explosive by adopting a CJ volume fuel chemical reaction model, wherein the volume V is in a CJ stateCJ0.357, viscosity coefficient 0.2;
step 400, constructing an elastic-plastic model;
step 500, calculating mass migration;
step 600, reassign flow field variables.
The embodiment can calculate the output characteristic of the mixed explosive when exploding in water, and the calculation precision meets the actual requirement and is beneficial to improving the design level.
Example 3:
according to the technical scheme, the method for calculating the combustion-to-detonation characteristic of the mixed explosive EX3 is used for calculating the combustion-to-detonation characteristic of the mixed explosive.
The mixed explosive EX3 comprises the following components in parts by weight: HMX: fluororubber (F-Rubber) 95: 5.
the model contains two materials: mixing explosive and aluminum.
Step 100, constructing a calculation model which comprises mass, momentum and energy conservation equations, wherein the number of calculation grids is 100 in the x direction, the number of calculation grids is 200 in the y direction, the total number is 20000, the initial temperature is 300K, the initial pressure is 0.1MPa, the calculation time step length is 0.05 mu s, and the total calculation time length is 25 mu s;
step 200, constructing a state equation, wherein the mixed explosive is described by adopting a solid phase HOM state equation when not reacting, and the parameters of the state equation are as follows: the coefficients F ═ 0.0000, G ═ -16.4687, H ═ -105.4688, I ═ -181.0623, and J ═ -129.3990;
the detonation reaction product is described by adopting a gas-phase HOM equation of state, and the parameters of the equation of state are as follows: the coefficients a-3.8489, B-2.6707, C-0.2268, D-0.0780 and E-0.0337, the coefficients K-1.5864, L-0.5244, M-0.0899, N-0.0077 and O-0.2605 × 10-3The coefficient Q is 7.1167, R is-0.5500, S is 0.0596, U is 0.0112 and X is-0.0101;
aluminum is described by a solid phase HOM state equation; the coefficients F ═ 0.0000, G ═ -20.7547, H ═ -115.6178, I ═ -177.7626, and J ═ -114.2528;
300, constructing an explosive chemical reaction model, describing the mixed explosive by adopting a Forest-Fire heterogeneous combustion and linear volume combustion chemical reaction model, wherein the Forest-Fire heterogeneous combustion chemical reaction model coefficient a0＝-7.7583×105、a1＝7.7280×105、a2＝-3.0193×105、a3＝5.9730×104、a4＝-6.4367×103、a5411.1578 and a6-15.2642, coefficient of linear volume combustion chemical reaction model, (S)0/V0)＝75.0，q＝2/3；c＝0.0092，b＝2/3；
Step 400, constructing an elastic-plastic model; the shear modulus of the aluminum is 29GPa, and the yield strength is 260 MPa;
step 500, calculating mass migration;
step 600, reassign flow field variables.
According to the embodiment, the characteristic of the mixed explosive of the combustion-detonation transition can be calculated, the calculation precision meets the actual requirement, and the improvement of the design level is facilitated.
Although the invention has been described herein with reference to a number of illustrative embodiments thereof, it should be understood that numerous other modifications and embodiments can be devised by those skilled in the art that will fall within the spirit and scope of the principles of this disclosure.

Claims (3)

1. A method for calculating explosion output characteristics of a mixed explosive is based on a fluid dynamic equation under a two-dimensional Euler coordinate system and is characterized by comprising the following calculation steps of:
step 100, constructing a calculation model:
constructing a mass conservation equation as shown in a formula 1.1 under a two-dimensional Euler coordinate system, a momentum conservation equation as shown in a formula 1.2 and a formula 1.3, and an energy conservation equation as shown in a formula 1.4;
conservation of mass equation:
Figure FDA0003150095320000011
conservation of momentum in the x direction equation:
Figure FDA0003150095320000012
conservation of momentum in the y-direction equation:
Figure FDA0003150095320000013
energy conservation equation:
Figure FDA0003150095320000014
formula 1.1-formula 1.4:
rho is density;
u is the speed;
e is energy;
s is stress offset;
p is the pressure intensity;
t is time;
alpha is a symmetric coefficient, alpha is 2 under a cylindrical coordinate system, and alpha is 1 under a rectangular coordinate system;
subscripts x and y are the components of the respective directions;
step 200, constructing a state equation:
the physical property change rule of the detonation product of the mixed explosive is described by adopting an HOM state equation, and the expression form is shown as formula 1.5:
Figure FDA0003150095320000021
in formula 1.5:
pCJthe pressure of detonation products in the CJ state;
eCJinternal energy of detonation product under CJ state;
TCJis the temperature of detonation products under CJ conditions;
subscript CJ indicates CJ status as a reference;
Vgis the volume of detonation products;
coefficients A, B, C, D and E, K, L, M, N and O, and Q, R, S, U and X are fitting coefficients of the respective physical quantities, respectively;
step 300, constructing an explosive chemical reaction model:
the method comprises an Arrhenius chemical reaction model, a Taylor wave chemical reaction model, a square wave impact chemical reaction model, a CJ volume combustion chemical reaction model and a Forest-Fire heterogeneous combustion chemical reaction model;
step 400, constructing an elastic-plastic model:
when the yield strength of the material is exceeded, an elastic-plastic model is used to describe the stress state of the object, expressed in the form of equation 1.6:
Figure FDA0003150095320000031
in formula 1.6:
f represents the principal stress;
Y0indicates the yield strength of the material;
step 500, mass migration is calculated:
transferring the mass according to the motion speed of the flow field, wherein the transferred mass flows into a new computational grid;
step 600, redistributing flow field variables:
flow field variables including density, mass fraction, x-direction velocity, y-direction velocity, and internal energy are redistributed according to the mass transfer calculation in step 500.
2. The method for calculating the explosion output characteristics of the mixed explosive according to claim 1, wherein in the step 300, the Forest-Fire heterogeneous combustion chemical reaction model further comprises a linear volume combustion chemical reaction model.
3. The method for calculating the explosion output characteristics of the mixed explosive according to claim 1, wherein the step 300 further comprises a step of adding a viscosity coefficient.
CN202110761644.2A 2021-07-06 2021-07-06 Mixed explosive explosion output characteristic calculation method Active CN113553739B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110761644.2A CN113553739B (en) 2021-07-06 2021-07-06 Mixed explosive explosion output characteristic calculation method

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110761644.2A CN113553739B (en) 2021-07-06 2021-07-06 Mixed explosive explosion output characteristic calculation method

Publications (2)

Publication Number Publication Date
CN113553739A true CN113553739A (en) 2021-10-26
CN113553739B CN113553739B (en) 2023-05-12

Family

ID=78131348

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110761644.2A Active CN113553739B (en) 2021-07-06 2021-07-06 Mixed explosive explosion output characteristic calculation method

Country Status (1)

Country Link
CN (1) CN113553739B (en)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114678078A (en) * 2022-03-11 2022-06-28 南京航空航天大学 CO (carbon monoxide)2-CH4Reforming foam reactor and optimization design method thereof
CN114861396A (en) * 2022-03-30 2022-08-05 西北核技术研究所 Mathematical model considering fine sand impact compression characteristic and heat absorption phase change and modeling method

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP2237175A1 (en) * 2009-04-01 2010-10-06 Livermore Software Technology Corporation Explosion simulation in finite element analysis
CN105676273A (en) * 2016-01-19 2016-06-15 北京理工大学 Design method for high-energy high-frequency explosive seismic source
CN106503379A (en) * 2016-10-28 2017-03-15 京工博创（北京）科技有限公司 A kind of gas burst emulation mode that is reacted based on adaptive simplifying with grid subdivision
CN109214082A (en) * 2018-09-03 2019-01-15 北京理工大学 A kind of high resolution numerical simulation method of near field underwater blast wave load
CN112699618A (en) * 2020-12-18 2021-04-23 赣江新区澳博颗粒科技研究院有限公司 Numerical simulation method for in-situ leaching process of ionic rare earth ore
CN112800601A (en) * 2021-01-19 2021-05-14 中国人民解放军陆军工程大学 Dissimilar metal explosion composite optimal combination parameter calculation method

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP2237175A1 (en) * 2009-04-01 2010-10-06 Livermore Software Technology Corporation Explosion simulation in finite element analysis
CN105676273A (en) * 2016-01-19 2016-06-15 北京理工大学 Design method for high-energy high-frequency explosive seismic source
CN106503379A (en) * 2016-10-28 2017-03-15 京工博创（北京）科技有限公司 A kind of gas burst emulation mode that is reacted based on adaptive simplifying with grid subdivision
CN109214082A (en) * 2018-09-03 2019-01-15 北京理工大学 A kind of high resolution numerical simulation method of near field underwater blast wave load
CN112699618A (en) * 2020-12-18 2021-04-23 赣江新区澳博颗粒科技研究院有限公司 Numerical simulation method for in-situ leaching process of ionic rare earth ore
CN112800601A (en) * 2021-01-19 2021-05-14 中国人民解放军陆军工程大学 Dissimilar metal explosion composite optimal combination parameter calculation method

Non-Patent Citations (5)

* Cited by examiner, † Cited by third party
Title
LIANG GE等: "Investigation of underwater explosion near composite structures using a combined RKDG-FEM approach" *
刘志跃等: "炸药落锤试验点火机理的数值分析", 《计算力学学报》 *
吴宗铎等: "基于Mie-Grüneisen状态方程的水下爆炸数值模拟" *
管怀安: "《爆破与安全》", 31 March 1984 *
蒋国平: "非均质凝聚炸药二维冲击起爆实验与Lagrange分析研究", 《中国优秀硕士学位论文全文数据库 基础科学辑》 *

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114678078A (en) * 2022-03-11 2022-06-28 南京航空航天大学 CO (carbon monoxide)2-CH4Reforming foam reactor and optimization design method thereof
CN114861396A (en) * 2022-03-30 2022-08-05 西北核技术研究所 Mathematical model considering fine sand impact compression characteristic and heat absorption phase change and modeling method

Also Published As

Publication number Publication date
CN113553739B (en) 2023-05-12

Similar Documents

Publication Publication Date Title
CN113553739A (en) Method for calculating explosion output characteristics of mixed explosive
Bammann et al. Failure in ductile materials using finite element methods
Abouaf et al. Finite element simulation of hot isostatic pressing of metal powders
CN107590336A (en) Method for numerical simulation of the fuel gas pipeline leakage to internal flow
CN110750933B (en) Accurate interface tracking processing method for coupling Lagrange particles and Euler method
Ren et al. Experimental study on dynamic buckling of submerged grid-stiffened cylindrical shells under intermediate-velocity impact
CN106777457A (en) A kind of solid rocket motor grain Structural Reliability Evaluation software systems
Yu et al. Aerodynamic analysis of projectile in gun system firing process
Goyal et al. Smooth particle hydrodynamics for bird-strike analysis using LS-DYNA
Tianbao et al. Large scale high precision computation for explosion and impact problems
Xu et al. Numerical method and simplified analytical model for predicting the blast load in a partially confined chamber
Mansouri et al. Buckling assessment of imperfect cylindrical shells under axial loads using a gep technique
Menshov et al. Numerical modeling of wave processes accompanying combustion of inhomogeneously distributed composite propellant
Sielicki et al. Implementation of sapper-blast-module, a rapid prediction software for blast wave properties
Zhou et al. Design of scaled model for dynamic characteristics of stiffened cylindrical shells based on equivalent similar method
CN110083946A (en) A kind of modified method of multistate model based on unconstrained optimization model
Brovko et al. AA Il’yushin’s scientific heritage and development of his ideas in mechanics
Čarnogurská et al. Modelling the flow conditions in the tunnel and its reduced model
Hu et al. A coupled computational framework for the transient heat transfer in a circular pipe heated internally with expanding heat sources
Valko et al. Simulation of an Air Shock Wave Using the Equations of State for the Jones–Wilkins–Lee Detonation Products
Xi et al. Prediction of Spatial Peak Overpressure Profile of Air Blast Shocks Using Multiple Linear Regression and Artificial Neural Network
Duan et al. On the transient numerical simulation of solid rocket motor by coupling quasi one-dimension internal flow with three-dimension propellant grain combustion
Baratta et al. Nonlinear versus linear stress-strain relations in engine turbulence modeling under swirl and squish flow conditions
Herreinstein et al. Piston motion mathematical model in the tube influenced by burning gas
Simoens et al. Analysis of the dynamic response of a controlled detonation chamber

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