CN108984996A - A kind of aluminum composition JWL state equation calculation method of parameters based on aluminium powder degree of reaction - Google Patents
A kind of aluminum composition JWL state equation calculation method of parameters based on aluminium powder degree of reaction Download PDFInfo
- Publication number
- CN108984996A CN108984996A CN201811136556.8A CN201811136556A CN108984996A CN 108984996 A CN108984996 A CN 108984996A CN 201811136556 A CN201811136556 A CN 201811136556A CN 108984996 A CN108984996 A CN 108984996A
- Authority
- CN
- China
- Prior art keywords
- state equation
- aluminum composition
- aluminium powder
- jwl
- degree
- 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
Links
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/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
Abstract
The present invention provides a kind of aluminum composition JWL state equation calculation method of parameters based on aluminium powder degree of reaction, this method utilizes unreacted aluminium powder in inert substance substitution aluminum composition detonation process, integrated use detonation theory calculates, water shooting numerical simulation and water shooting are tested, that is, can determine aluminum composition detonation product JWL state equation parameter.This method has many advantages, such as that accuracy is high, with strong points, good reliability, and error big defect when overcoming theoretical calculation prediction aluminum composition JWL state equation parameter has certain use value for the design of aluminum composition formula and explosion property assessment.
Description
Technical field
The present invention relates to a kind of JWL state equation calculation method of parameters more particularly to a kind of containing based on aluminium powder degree of reaction
Al explosive JWL state equation calculation method.
Background technique
JWL state equation more can accurately describe the swollen of detonation products as a kind of dynamics state equation
Swollen driving acting process, is widely used at present in various detonation numbered analog simulation software for calculation.Determine explosive JWL state
There are two ways to equation, first is that can accurately obtain JWL state equation parameter by cylinder test, another is to utilize heat
Mechanics program determines to calculate the constant entropy expansion data of detonation product.The method of cylinder test is although accurate, testing expenses
And the time is at high cost, and under different densities, different compositions, the JWL state equation parameter of detonation products is different.Therefore it manages
The assessment in formulation designs stage is carried out with regard to extremely important by acquisition JWL state equation parameter is calculated.
For non-ideal constraints such as aluminum compositions, the detonation reaction time is longer (having reached ms grades), and there are a fixed widths
The reaction zone of degree.Therefore, the constant entropy expansion data of aluminum composition are directly calculated not using detonation parameter calculation of thermodynamics program
Meet truth.Only consider aluminium powder degree of reaction when, the aluminum composition JWL state equation parameter of calculating can just obtain compared with
For accurate result.
" the RDX base aluminum composition JWL state equation parameter prediction based on KHT program studies [J] " Beijing Institute of Technology is learned
Report, 2013,33 (3): 239-243 discloses a kind of based on KHT programmed algorithm calculating RDX base aluminum composition JWL state equation ginseng
Several methods.Program is based on KHT equation, from explosive composition, solves explosive using detonation product free energy minimization principle
Constant entropy expansion data, then it is fitted to the form of JWL state equation.Calculated result shows that the JWL state equation that the algorithm acquires misses
Difference about 15%, error is larger, is unsatisfactory for the needs assessed aluminum composition.
Summary of the invention
For deficiency existing for existing aluminum composition JWL state equation calculation method of parameters, the present invention provides one kind and is based on
The aluminum composition JWL state equation calculation method of aluminium powder degree of reaction, for estimating the JWL state equation parameter of aluminum composition.
The present invention provides a kind of aluminum composition JWL state equation calculation method of parameters based on aluminium powder degree of reaction, calculates stream
Journey is as shown in Figure 1, comprising the following steps:
Step S1 chooses inert substance, replaces the part aluminium powder in aluminum composition using inert substance, design aluminium powder is different
The aluminum composition formula of degree of reaction；
Step S2, by the aluminum composition of differential responses degree formula molecular formula, density, the enthalpy of formation, detonation product thermodynamics
Function and solid state equation parameter are input in BKW detonation parameter calculation of thermodynamics program, and application program calculates aluminium powder not
With the constant entropy expansion curve of degree of reaction aluminum composition；
Step S3, by different constant entropy expansion curve data fits at the form of JWL state equation parameter；
Step S4 constructs explosive underwater explosion computation model, choosing using nonlinear Finite Element Numerical simulation software AUTODYN
It selects material model and modifies the JWL state equation parameter of explosive material, fire point is set and the solution time is each to 100g~10kg
The aluminum composition formula of differential responses degree carries out Underwater Explosion test simulation, as shown in Figure 2.Solve different moments, different location
The shock wave pressure curve at place, as shown in Figure 3.Shock wave pressure peak value is obtained by shock wave pressure time-history curves；
Step S5, the one group of shock wave pressure peak value and experiment value that will acquire compare and analyze, and choose calculated value and reality
Test the consistent JWL state equation parameter for being determined as the aluminum composition of value；
Step S1 inert substance of the invention preferentially selects lithium fluoride (LiF) material.LiF is as inert material, as fried
The group timesharing of medicine is not involved in the chemical reaction of explosive, while its density is close with aluminium powder, can replace in mass ratio for 1:1.
The present invention calculate when substance thermodynamic function and solid state equation be preferably respectively adopted entropy-temperature equation with
Cowan state equation describes.
The preferred BKW equation of detonation parameter calculation of thermodynamics program in step S2 of the invention, preferred α=0.5 of equation coefficient,
β=0.16, θ=400, κ=10.909, the constant entropy expansion for solving detonation product using the free energy minimization principle of detonation product are bent
Line number evidence.
The preferred AUTODYN software of the Nonlinear Finite meta software that the present invention uses, simulation Underwater Explosion test use fried
The ball shape pack of the preferred 1kg~5kg of medicine shape, the preferred unlimited water area in the waters of use utilize Eulers algorithm, Shock shape
The description of state equation, establishes the preferably one-dimensional wedge shaped element of model of use, unit size preferably 1~2mm, 6000~8000 lists
Member.It solves at the quick-fried heart 1m~5m of position detection point preferred distance, solves moment preferred 0ms~6ms.
Advantages of the present invention: compared with prior art, the present invention is the JWL shape based on aluminium powder degree of reaction fitting aluminum composition
State equation parameter, the reliable error of calculated result is small, and data usability is good；The explosive JWL state equation parameter that the present invention fits
It can be used for the detonation numerical simulation calculation of explosive, reduce dependence and requirement to experimental technique.
Detailed description of the invention
Fig. 1 is a kind of aluminum composition JWL state equation calculation flow chart based on aluminium powder degree of reaction.
Fig. 2 is water shooting computation model.
Fig. 3 is to calculate the shock wave pressure time-history curves away from quick-fried heart 1m observation point.
Specific embodiment
Below with reference to embodiment, the present invention will be further described, but the present invention is not limited by the following example.
Embodiment 1
Step S1, from the entropy searched in online JANAF thermodynamic data table under LiF (cr) temperature range 0K~2500K
Value, and entropy-temperature equation is fitted to according to these data.LiF (cr) is searched from online LASL impact Hugoniot database
Shock Compression parameter, and be fitted to the form of Cowan solid state equation；
Step S2, [heisuojin (RDX)/TNT/LiF/ wax mass ratio is 45/ to the aluminum composition each component for being 0 by degree of reaction
30/20/5] molecular formula, density, the enthalpy of formation, detonation product and the inert substance of fitting (H2O、CO2、N2、H2、O2、CO、NH3、
H、NO、OH、CH4, C, LiF) thermodynamic function parameter and solid state equation parameter be input to the detonation based on BKW equation ginseng
In number calculation of thermodynamics program, the constant entropy expansion curve of the aluminum composition under the degree of reaction is calculated；
Step S3, by the aluminum composition constant entropy expansion curve data fit under the degree of reaction at JWL state equation form；
Step S4 starts AUTODYN software, creates a new computation model file, creates 1kg spheroidal charging unlimited
The bang model in waters inputs the JWL parameter of explosive.Model use one-dimensional wedge shaped element, unit size 1mm, totally 6000
Unit.Water uses Eulers algorithm, is described using Shock state equation, is being that 5 are arranged at 1,2,3,4,5m apart from quick-fried heart d
Observation point.
Obtained shock wave pressure curve is analyzed and processed, obtains shock wave pressure peak value by step S5.By shock wave
Pressure peak is compared with the underwater explosive test data of 1kg spherical shape explosive pack, compares the meter acquired under the degree of reaction
The error of the value that counts and experiment value；
Step S6 repeats step S2~S5, solves the aluminum composition [heisuojin (RDX)/TNT/ that aluminium powder degree of reaction is 10%
Al/LiF/ wax mass ratio is 45/30/2/18/5].Based on being carried out to the explosive under differential responses degree under the formula by above step
Calculate, aluminium powder degree of reaction increases by 10% every time, solve always to aluminium powder degree of reaction be 100% aluminum composition [heisuojin (RDX)/
TNT/Al/ wax mass ratio is 45/30/20/5].
The explosive under differential responses degree under the formula is calculated by above step, until the calculated result that acquires with
Experiment value is consistent, that is, acquires the degree of reaction and JWL state equation parameter of the aluminium powder under the formula.
Embodiment 2
Step S1, from the entropy searched in online JANAF thermodynamic data table under LiF (cr) temperature range 0K~2500K
Value, and entropy-temperature equation is fitted to according to these data.LiF (cr) is searched from online LASL impact Hugoniot database
Shock Compression parameter, and be fitted to the form of Cowan solid state equation；
Step S2, [heisuojin (RDX)/TNT/LiF/ wax mass ratio is 40/ to the aluminum composition each component for being 0 by degree of reaction
38/17/5] molecular formula, density, the enthalpy of formation, detonation product and the inert substance of fitting (H2O、CO2、N2、H2、O2、CO、NH3、
H、NO、OH、CH4, C, LiF) thermodynamic function parameter and solid state equation parameter be input to the detonation based on BKW equation ginseng
In number calculation of thermodynamics program, the constant entropy expansion curve of the aluminum composition under the degree of reaction is calculated；
Step S3, by the aluminum composition constant entropy expansion curve data fit under the degree of reaction at JWL state equation form；
Step S4 starts AUTODYN software, creates a new computation model file, creates 1kg spheroidal charging unlimited
The bang model in waters inputs the JWL parameter of explosive.Model use one-dimensional wedge shaped element, unit size 1mm, totally 6000
Unit.Water uses Eulers algorithm, is described using Shock state equation, is being that 5 are arranged at 1,2,3,4,5m apart from quick-fried heart d
Observation point.
Obtained shock wave pressure curve is analyzed and processed, obtains shock wave pressure peak value by step S5.By shock wave
Pressure peak is compared with the underwater explosive test data of 1kg spherical shape explosive pack, compares the meter acquired under the degree of reaction
The error of the value that counts and experiment value；
Step S6 repeats step S2~S5, solves the aluminum composition [heisuojin (RDX)/TNT/ that aluminium powder degree of reaction is 10%
Al/LiF/ wax mass ratio is 40/38/1.7/15.3/5].The explosive under differential responses degree under the formula is carried out by above step
It calculates, aluminium powder degree of reaction increases by 10% every time, solves the aluminum composition [heisuojin for being 100% to aluminium powder degree of reaction always
(RDX)/TNT/Al/ wax mass ratio is 40/38/17/5].
The explosive under differential responses degree under the formula is calculated by above step, until the calculated result that acquires with
Experiment value is consistent, that is, acquires the degree of reaction and JWL state equation parameter of the aluminium powder under the formula.
Embodiment 3
Step S1, from the entropy searched in online JANAF thermodynamic data table under LiF (cr) temperature range 0K~2500K
Value, and entropy-temperature equation is fitted to according to these data.LiF (cr) is searched from online LASL impact Hugoniot database
Shock Compression parameter, and be fitted to the form of Cowan solid state equation；
Step S2, [heisuojin (RDX)/TNT/LiF/ wax mass ratio is the aluminum composition each component for being 0 by degree of reaction
60.8/19/17/3.2 molecular formula, density, the enthalpy of formation, detonation product and the inert substance of fitting (H]2O、CO2、N2、H2、O2、
CO、NH3、H、NO、OH、CH4, C, LiF) thermodynamic function parameter and solid state equation parameter be input to based on BKW equation
In detonation parameter calculation of thermodynamics program, the constant entropy expansion curve of the aluminum composition under the degree of reaction is calculated；
Step S3, by the aluminum composition constant entropy expansion curve data fit under the degree of reaction at JWL state equation form；
Step S4 starts AUTODYN software, creates a new computation model file, creates 1kg spheroidal charging unlimited
The bang model in waters inputs the JWL parameter of explosive.Model use one-dimensional wedge shaped element, unit size 1mm, totally 6000
Unit.Water uses Eulers algorithm, is described using Shock state equation, is being that 5 are arranged at 1,2,3,4,5m apart from quick-fried heart d
Observation point.
Obtained shock wave pressure curve is analyzed and processed, obtains shock wave pressure peak value by step S5.By shock wave
Pressure peak is compared with the underwater explosive test data of 1kg spherical shape explosive pack, compares the meter acquired under the degree of reaction
The error of the value that counts and experiment value；
Step S6 repeats step S2~S5, solves the aluminum composition [heisuojin (RDX)/TNT/ that aluminium powder degree of reaction is 10%
Al/LiF/ wax mass ratio is 60.8/19/1.7/15.3/3.2].By above step to the explosive under differential responses degree under the formula
It is calculated, aluminium powder degree of reaction increases by 10% every time, solves the aluminum composition [heisuojin for being 100% to aluminium powder degree of reaction always
(RDX)/TNT/Al/ wax mass ratio is 60.8/19/17/3.2].
The explosive under differential responses degree under the formula is calculated by above step, until the calculated result that acquires with
Experiment value error is no more than 10%, that is, acquires the degree of reaction and JWL state equation parameter of the aluminium powder under the formula.
Utilize the shock wave for the JWL state equation simulation aluminum composition underwater explosion that embodiment 1- embodiment 3 is calculated
Can, list file names with same recipe using the calculated result of algorithm before improving, as shown in table 1.
1 shock wave energy calculated result of table
As a result: the error calculated under 3 kinds of embodiments is respectively 2.9%, 2.0%, 8.1%.Than the error before improving
Substantially reduce.
Conclusion: JWL state equation parameter simulation explosive underwater blast wave energy calculated value is calculated based on aluminium powder degree of reaction
It is substantially reduced with test value error.The precision of aluminum composition underwater explosion numerical simulation is substantially improved.
Claims (4)
1. a kind of aluminum composition JWL state equation calculation method of parameters based on aluminium powder degree of reaction, which is characterized in that utilize inertia
Substance substitute part aluminium powder design aluminium powder differential responses degree aluminum composition, the calculation method the following steps are included:
Step S1 chooses inert substance, replaces the part aluminium powder in aluminum composition using inert substance, designs aluminium powder differential responses
The aluminum composition formula of degree；
Step S2, by the aluminum composition of differential responses degree formula molecular formula, density, the enthalpy of formation, detonation product thermodynamic function
It is input in BKW detonation parameter calculation of thermodynamics program with solid state equation parameter, it is different anti-that application program calculates aluminium powder
The constant entropy expansion curve of response aluminum composition；
Step S3, by different constant entropy expansion curve data fits at the form of JWL state equation parameter；
Step S4 constructs explosive underwater explosion computation model using nonlinear Finite Element Numerical simulation software AUTODYN, selects material
Material model and the JWL state equation parameter for modifying explosive material, are arranged fire point and the solution time is variant to 100g~10kg
The aluminum composition formula of degree of reaction carries out Underwater Explosion test simulation, the shock wave pressure at solution different moments, different location
Peak value；
Step S5, the one group of shock wave pressure peak-data and experiment value that will acquire compare and analyze, and choose calculated value and reality
Test the consistent JWL state equation parameter for being determined as the aluminum composition of value.
2. JWL state equation calculation method of parameters according to claim 1, it is characterised in that select inert material in step S1
Matter substitutes unreacted aluminium powder to design the aluminum composition formula of differential responses degree.
3. JWL state equation calculation method of parameters according to claim 1, it is characterised in that detonation product in step S2
Thermodynamic function is described using entropy-temperature equation description, solid state equation using Cowan solid state equation.
4. JWL state equation calculation method of parameters according to claim 2, it is characterised in that inert substance using LiF come
Description.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811136556.8A CN108984996B (en) | 2018-09-28 | 2018-09-28 | Aluminum powder reactivity-based aluminum-containing explosive JWL state equation parameter calculation method |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811136556.8A CN108984996B (en) | 2018-09-28 | 2018-09-28 | Aluminum powder reactivity-based aluminum-containing explosive JWL state equation parameter calculation method |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN108984996A true CN108984996A (en) | 2018-12-11 |
| CN108984996B CN108984996B (en) | 2022-12-13 |
Family
ID=64544118
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201811136556.8A Active CN108984996B (en) | 2018-09-28 | 2018-09-28 | Aluminum powder reactivity-based aluminum-containing explosive JWL state equation parameter calculation method |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN108984996B (en) |
Cited By (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN112946229A (en) * | 2021-01-29 | 2021-06-11 | 北京理工大学 | Method for acquiring performance of aluminum-containing explosive based on cylinder-sheet device |
| CN113177348A (en) * | 2021-05-25 | 2021-07-27 | 北京理工大学 | Quasi-isentropic model construction method and working performance evaluation method for aluminum-containing explosive |
| CN113593650A (en) * | 2021-07-06 | 2021-11-02 | 西安近代化学研究所 | Method for calculating detonation parameters of mixed explosive |
| CN114861481A (en) * | 2022-03-30 | 2022-08-05 | 西北核技术研究所 | Calculation method for any-configuration explosion source ground motion steady-state effect |
| CN116046840A (en) * | 2022-12-13 | 2023-05-02 | 西南交通大学 | Experimental device for aluminum-containing explosive product state equation under millisecond time scale |
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20030037692A1 (en) * | 2001-08-08 | 2003-02-27 | Liqing Liu | Use of aluminum in perforating and stimulating a subterranean formation and other engineering applications |
| CN102592053A (en) * | 2012-01-09 | 2012-07-18 | 北京理工大学 | Method for determining detonation parameters of liquid explosive |
-
2018
- 2018-09-28 CN CN201811136556.8A patent/CN108984996B/en active Active
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20030037692A1 (en) * | 2001-08-08 | 2003-02-27 | Liqing Liu | Use of aluminum in perforating and stimulating a subterranean formation and other engineering applications |
| CN102592053A (en) * | 2012-01-09 | 2012-07-18 | 北京理工大学 | Method for determining detonation parameters of liquid explosive |
Non-Patent Citations (1)
| Title |
|---|
| 项大林等: "基于KHT程序的RDX基含铝***JWL状态方程参数预测研究", 《北京理工大学学报》 * |
Cited By (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN112946229A (en) * | 2021-01-29 | 2021-06-11 | 北京理工大学 | Method for acquiring performance of aluminum-containing explosive based on cylinder-sheet device |
| CN112946229B (en) * | 2021-01-29 | 2023-03-14 | 北京理工大学 | Method for acquiring performance of aluminum-containing explosive based on cylinder-sheet device |
| CN113177348A (en) * | 2021-05-25 | 2021-07-27 | 北京理工大学 | Quasi-isentropic model construction method and working performance evaluation method for aluminum-containing explosive |
| CN113177348B (en) * | 2021-05-25 | 2022-07-12 | 北京理工大学 | Quasi-isentropic model construction method and working performance evaluation method for aluminum-containing explosive |
| CN113593650A (en) * | 2021-07-06 | 2021-11-02 | 西安近代化学研究所 | Method for calculating detonation parameters of mixed explosive |
| CN113593650B (en) * | 2021-07-06 | 2023-09-12 | 西安近代化学研究所 | Mixed explosive detonation parameter calculation method |
| CN114861481A (en) * | 2022-03-30 | 2022-08-05 | 西北核技术研究所 | Calculation method for any-configuration explosion source ground motion steady-state effect |
| CN116046840A (en) * | 2022-12-13 | 2023-05-02 | 西南交通大学 | Experimental device for aluminum-containing explosive product state equation under millisecond time scale |
| CN116046840B (en) * | 2022-12-13 | 2023-08-29 | 西南交通大学 | Experimental device for aluminum-containing explosive product state equation under millisecond time scale |
Also Published As
| Publication number | Publication date |
|---|---|
| CN108984996B (en) | 2022-12-13 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN108984996A (en) | A kind of aluminum composition JWL state equation calculation method of parameters based on aluminium powder degree of reaction | |
| Trzciński et al. | Detonation properties of 1, 1-diamino-2, 2-dinitroethene (DADNE) | |
| Wharton et al. | Airblast TNT equivalence for a range of commercial blasting explosives | |
| CN102750425B (en) | Simulation method of texture evolvement of heat affected zone in welding process | |
| Fried et al. | BKWC: an empirical BKW parametrization based on cylinder test data | |
| Zavaliangos et al. | Thermo-elasto-viscoplasticity of isotropic porous metals | |
| Vrabec et al. | Vapor–liquid equilibria of the ternary mixture CH4+ C2H6+ CO2 from molecular simulation | |
| Maynadier et al. | Thermo-mechanical description of phase transformation in Ni-Ti Shape Memory Alloy | |
| Hara et al. | Experimental and Theoretical Investigation of the Heat of Combustion of RDX-based Propellants | |
| Brommesson et al. | Experimental observations and modelling of cyclic and relaxation behaviour of the Ni-based superalloy Haynes 282 | |
| Sanchidrián et al. | Calculation of the energy of explosives with a partial reaction model. Comparison with cylinder test data | |
| Berry et al. | Anelastic behavior of hydrogenated amorphous metals | |
| Opěla et al. | New model predicting flow curves in wide range of thermomechanical conditions of 38MnVS6 steel | |
| Zhang et al. | Prediction of the explosion effect of aluminized explosives | |
| Belaada et al. | A melt-cast composition based on NTO and FOX-7 | |
| CN112033998A (en) | Thermal insulation acceleration calorimetry-based explosive material thermal stability grading method | |
| Alister et al. | Viscoplastic and temperature behavior of Zn–Cu–Ti alloy sheets: experiments, characterization, and modeling | |
| CN109632613B (en) | Estimation method for safe service life of JHLDU-1 explosive charge of press-loaded penetration explosive | |
| 위성욱 et al. | Prediction of Low-Cycle Fatigue Life of In738LC Using Plastic Strain Energy Density | |
| Biss et al. | On the use of composite charges to determine insensitive explosive material properties at the laboratory scale | |
| Santacreu et al. | Thermomechanical fatigue of stainless steels for automotive exhaust systems | |
| Hershkowitz et al. | Evaluation by a modified cylinder test of metal acceleration by non-ideal explosives containing ammonium nitrate | |
| Pei et al. | High temperature constitutive model of q345B steel | |
| Baker et al. | Stable detonation velocities for aluminized explosives | |
| JP2020066049A (en) | Forging structure prediction method, forging structure prediction device, and plastically processed product manufacturing method |
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 |