CN108627629A - Method for determining overpressure peak value of axial and radial shock waves of cylindrical explosive - Google Patents

Method for determining overpressure peak value of axial and radial shock waves of cylindrical explosive Download PDF

Info

Publication number
CN108627629A
CN108627629A CN201810436867.XA CN201810436867A CN108627629A CN 108627629 A CN108627629 A CN 108627629A CN 201810436867 A CN201810436867 A CN 201810436867A CN 108627629 A CN108627629 A CN 108627629A
Authority
CN
China
Prior art keywords
explosive
overpressure
detonation
axial
cylinder
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
CN201810436867.XA
Other languages
Chinese (zh)
Other versions
CN108627629B (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.)
Army Engineering University of PLA
Original Assignee
Army Engineering University of PLA
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 Army Engineering University of PLA filed Critical Army Engineering University of PLA
Priority to CN201810436867.XA priority Critical patent/CN108627629B/en
Publication of CN108627629A publication Critical patent/CN108627629A/en
Application granted granted Critical
Publication of CN108627629B publication Critical patent/CN108627629B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N33/00Investigating or analysing materials by specific methods not covered by groups G01N1/00 - G01N31/00
    • G01N33/22Fuels; Explosives
    • G01N33/227Explosives, e.g. combustive properties thereof
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01LMEASURING FORCE, STRESS, TORQUE, WORK, MECHANICAL POWER, MECHANICAL EFFICIENCY, OR FLUID PRESSURE
    • G01L5/00Apparatus for, or methods of, measuring force, work, mechanical power, or torque, specially adapted for specific purposes
    • G01L5/14Apparatus for, or methods of, measuring force, work, mechanical power, or torque, specially adapted for specific purposes for measuring the force of explosions; for measuring the energy of projectiles

Landscapes

  • Chemical & Material Sciences (AREA)
  • Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Analytical Chemistry (AREA)
  • Medicinal Chemistry (AREA)
  • Food Science & Technology (AREA)
  • Biochemistry (AREA)
  • General Health & Medical Sciences (AREA)
  • Combustion & Propulsion (AREA)
  • Immunology (AREA)
  • Pathology (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)
  • Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)

Abstract

The invention discloses a method for determining an overpressure peak value of axial and radial shock waves of a cylindrical explosive, which comprises the following steps of: (10) obtaining original parameters: acquiring original parameters of the cylindrical explosive, including the overall dimension, the mass, the explosive type and the overpressure calculation distance; (20) calculating intermediate parameters: calculating intermediate parameters of the cylindrical explosive according to the original parameters, wherein the intermediate parameters comprise equivalent TNT mass, proportional burst length and length-diameter ratio; (30) determining an overpressure peak value: and calculating to obtain the overpressure peak values of the axial shock waves and the radial shock waves of the cylindrical explosive according to the intermediate parameters and an overpressure peak value fitting formula. The method for determining the overpressure peak value of the axial and radial shock waves of the cylindrical explosive is not only suitable for cylindrical explosives with one end being detonated, but also suitable for cylindrical explosives with the center being detonated, and meanwhile, when the length-diameter ratio is changed, the result is accurate and reliable.

Description

Cylinder explosive axial-radial shock wave overpressure peak determines method
Technical field
The invention belongs to explosive performance test ing technical fields, are especially applicable not only to the cylinder explosive of one end detonation, Be also applied for center detonation cylinder explosive, meanwhile, when draw ratio becomes, as a result accurately, reliable cylinder explosive axial direction Radial impact wave peak overpressure determines method.
Background technology
When being designed with the structure and design explosive test scheme of antiknock requirement, or estimate the influence size of the attack of terrorism When, the power for estimating explosive charge is highly important.The explosive of cylindrical shape is common, cylinder on civilian and military There is very big difference in the peak overpressure that body explosive charge generates, on different spatial positions especially when apart from closer It waits, peak overpressure is exactly that the air pressure of explosion generation subtracts the maximum value after atmospheric pressure.It summarizes and has the study found that ratio It is quick-fried away from approximately greater than 4m/kg1/3After, explosive shape can just ignore the influence that peak overpressure is distributed substantially, you can with The shock wave overpressure peak generated with the spherical explosive of phase homogenous quantities is for estimating cylinder explosive charge power.
The shock wave overpressure peak predictor formula of existing cylinder explosive is mainly by being fitted a large amount of experimental data It gets, is mainly used for estimating the shock wave overpressure peak at different distance in axially and radially both direction.
Existing cylinder peak overpressure predictor formula is substantially polynomial formRadially very Rare formula considers the influence of draw ratio (L/D), the influence of no formula consideration draw ratio in axial direction.The cylinder of different draw ratios The shock wave overpressure peak distribution that body explosive generates is clearly to have differences, draw ratio for radial peak overpressure influence compared with It is small, it is very big for axial peak overpressure influence, but due to the discreteness of field test data, length is considered in fitting formula Diameter than influence it is relatively difficult, the reason for causing test data discreteness big has that experimental enviroment is complicated and changeable, explosive property is unstable It is fixed etc..In addition the existing experiment for being fitted cylinder explosive peak overpressure predictor formula is all one end detonation, without center The measurement of detonation is tested, also without the peak overpressure predictor formula of center detonation.
Therefore, problem of the existing technology is:It is fried that existing shock wave overpressure peak determines that method is only used for cylinder The peak overpressure prediction of medicine one end detonation, it is impossible to be used in the peak value prediction of center detonation;Existing axial direction peak overpressure prediction technique In do not account for the influence of draw ratio, when draw ratio changes, as a result accuracy is poor.
Invention content
The purpose of the present invention is to provide cylinder explosive axial-radial shock wave overpressure peaks to determine method, is not only applicable in In the cylinder explosive of one end detonation, and it is also applied for the cylinder explosive of center detonation, meanwhile, when draw ratio becomes, knot Fruit is accurate, reliable.
Realize that the technical solution of the object of the invention is:
A kind of cylinder explosive axial-radial shock wave overpressure peak determines method, includes the following steps:
(10) initial parameter obtains:Obtain cylinder explosive initial parameter, including appearance and size, quality, explosive type, Superpressure calculates distance;
(20) intermediate parameters calculate:According to initial parameter, the intermediate parameters of cylinder explosive are calculated, including equivalent TNT mass, ratio are quick-fried away from, draw ratio;
(30) peak overpressure determines:According to intermediate parameters and peak overpressure fitting formula, cylinder explosive axis is calculated To with radial impact wave peak overpressure.
Compared with prior art, the present invention its remarkable advantage is:
1, applied widely:The peak overpressure predictor formula of the cylinder explosive of center detonation axially and radially is provided, It is applicable not only to the cylinder explosive of one end detonation, and is also applied for the cylinder explosive of center detonation;
2, result is more acurrate:The influence that draw ratio is considered in axial peak overpressure predictor formula, when draw ratio becomes, knot Fruit is accurate, reliable.
Description of the drawings
Fig. 1 is that axially and radially shock wave overpressure peak determines the main flow chart of method to cylinder explosive of the present invention.
Fig. 2 is to use the radial peak overpressure of cylinder explosive one end detonation and draw ratio determined by the method for the present invention for 1/ The comparison diagram of 1 test data
Fig. 3 is to use the radial peak overpressure of cylinder explosive one end detonation and draw ratio determined by the method for the present invention for 4/ The comparison diagram of 1 test data
Fig. 4 is using cylinder explosive one end detonation axial proximal peak overpressure and draw ratio determined by the method for the present invention For the comparison diagram of 1/1 test data
Fig. 5 is using cylinder explosive one end detonation axial proximal peak overpressure and draw ratio determined by the method for the present invention For the comparison diagram of 4/1 test data
Fig. 6 is using cylinder explosive one end detonation axial distal end peak overpressure and draw ratio determined by the method for the present invention For the comparison diagram of 1/1 test data
Fig. 7 is using cylinder explosive one end detonation axial distal end peak overpressure and draw ratio determined by the method for the present invention For the comparison diagram of 4/1 test data
Fig. 8 is to use the radial peak overpressure of cylinder explosive center detonation and draw ratio determined by the method for the present invention for 1/ The comparison diagram of 1 to 4/1 analogue data
Fig. 9 is to use the axial peak overpressure of cylinder explosive center detonation and draw ratio determined by the method for the present invention for 1/ The comparison diagram of 1 to 4/1 analogue data
Figure 10 is the schematic diagram of cylinder explosive difference measurement direction.
Figure 11 is the detonation of cylinder explosive center and one end detonation schematic diagram.
Specific implementation mode
As shown in Figure 1, cylinder explosive axial-radial shock wave overpressure peak of the present invention determines method, including walk as follows Suddenly:
(10) initial parameter obtains:Obtain cylinder explosive initial parameter, including appearance and size, quality, explosive type, Superpressure calculates distance;
(20) intermediate parameters calculate:According to initial parameter, the intermediate parameters of cylinder explosive are calculated, including equivalent TNT mass, ratio are quick-fried away from, draw ratio;
It is described that ((20) intermediate parameters calculate step:
(21) equivalent TNT Mass Calculations:According to the corresponding TNT equivalent coefficients k of explosive type, equivalent TNT is calculated as follows Quality W,
W=k × m,
In formula, m is the quality of cylinder explosive, unit kg;
(22) distance R is calculated according to superpressure, it is quick-fried away from Z to be calculated as follows ratio,
In formula, W is equivalent TNT mass, and R is distance of the pressure measurement point apart from explosive geometric center, unit m.
(30) peak overpressure determines:According to intermediate parameters and peak overpressure fitting formula, cylinder explosive axis is calculated To with radial impact wave peak overpressure.
(30) peak overpressure determines that step includes:
(31) radial peak overpressure calculates:The radial impact of the detonation of cylinder explosive center and one end detonation is calculated as follows Wave peak overpressure, unit kPa,
(32) detonation axial peak overpressure in center calculates:The axial impact wave of center detonation cylinder explosive is calculated as follows Peak overpressure, unit kPa,
(33) detonation axial distal end peak overpressure in one end calculates:Be calculated as follows one end detonation cylinder explosive axially away from The peak overpressure of fire point one end:
PAxis, one end=1.4PAxis, center,
(34) detonation axial proximal peak overpressure in one end calculates:Be calculated as follows one end detonation cylinder explosive axially adjacent to The peak overpressure of fire point one end:
PAxis, one end=0.5PAxis, center,
In above-mentioned formula, Z is that ratio is quick-fried away from draw ratio L/D, L is the height of cylinder, and D is the straight of cylinder bottom surface circle Diameter, a, b are mediants, and e is scientific notation symbol, P0It is atmospheric pressure.
It is described that ((20) intermediate parameters calculate step:
(21) equivalent TNT Mass Calculations:According to the corresponding TNT equivalent coefficients k of explosive type, equivalent TNT is calculated as follows Quality W,
W=k × m,
In formula, m is the quality of cylinder explosive, unit kg;
(22) distance R is calculated according to superpressure, it is quick-fried away from Z to be calculated as follows ratio,
In formula, W is equivalent TNT mass, and R is distance of the pressure measurement point apart from explosive geometric center, unit m.
To verify effectiveness of the invention and accuracy, below with using the different modes detonation of the method for the present invention determination The experimental data that peak overpressure is respectively 1/1 and 4/1 with draw ratio compares.
Figure 10 is the schematic diagram of cylinder explosive difference measurement direction.
Figure 11 is the detonation of cylinder explosive center and one end detonation schematic diagram.
Fig. 2 is to use the radial peak overpressure of cylinder explosive one end detonation and draw ratio determined by the method for the present invention for 1/ The comparison diagram of 1 test data
Fig. 3 is to use the radial peak overpressure of cylinder explosive one end detonation and draw ratio determined by the method for the present invention for 4/ The comparison diagram of 1 test data
Fig. 4 is using cylinder explosive one end detonation axial proximal peak overpressure and draw ratio determined by the method for the present invention For the comparison diagram of 1/1 test data
Fig. 5 is using cylinder explosive one end detonation axial proximal peak overpressure and draw ratio determined by the method for the present invention For the comparison diagram of 4/1 test data
Fig. 6 is using cylinder explosive one end detonation axial distal end peak overpressure and draw ratio determined by the method for the present invention For the comparison diagram of 1/1 test data
Fig. 7 is using cylinder explosive one end detonation axial distal end peak overpressure and draw ratio determined by the method for the present invention For the comparison diagram of 4/1 test data
Fig. 8 is to use the radial peak overpressure of cylinder explosive center detonation and draw ratio determined by the method for the present invention for 1/ The comparison diagram of 1 to 4/1 analogue data
Fig. 9 is to use the axial peak overpressure of cylinder explosive center detonation and draw ratio determined by the method for the present invention for 1/ The comparison diagram of 1 to 4/1 analogue data
Test data bibliography includes:
1、Plooster M.N.:Blast effects from cylindrical explosive charges: Experimental measurements.Report NWC TP 6382.Naval Report Centre,China Lake, California 93555.1982.
2、Knock C,Davies N,Reeves T.Predicting Blast Waves from the Axial Direction of a Cylindrical Charge[J].Propellants Explosives Pyrotechnics, 2015,40(2):169-179.
3、Knock C,Davies N.Predicting the Impulse from the Curved Surface of Detonating Cylindrical Charges[J].Propellants Explosives Pyrotechnics,2011,36 (2):105-109.
4、Knock C,Davies N.Blast waves from cylindrical charges[J].Shock Waves,2013, 23(4):337-343.
Existing peak overpressure determines that the radial prediction that method detonates for one end is relatively ripe, but for axis To peak overpressure prediction do not account for the influence of draw ratio, a part of reason may be since test data discreteness is big, no It is mixed in together with the corresponding data of draw ratio, it is difficult to observe the rule of draw ratio influence and propose accurate predictor formula. Can be seen that from Fig. 2 to Fig. 7 using determined by the method for the present invention cylinder explosive one end detonate axial proximal, distal end, with And radial peak overpressure is sufficiently close to existing test data, in addition to the cylinder explosive axial distal end that draw ratio is 1/1 in Fig. 6 It tests peak overpressure and is more than formula calculated value, possible reason has that experimental enviroment is complicated and changeable, explosive property is unstable, sensor Precision difference, data processing human error etc., axis when illustrating the method for the present invention for determining the variation of cylinder explosive draw ratio To peak overpressure be accurately and reliably.
The present invention is applicable not only to the cylinder explosive of one end detonation, and the cylinder for being also applied for center detonation is fried Medicine.Due to the perfect test data for not having cylinder explosive center to detonate in open source literature, using occupying international military industry The numerical simulation software AUTODYN in 80% or more market is simulated, and is verified using one end detonation test data, is obtained Reliable numerical model, the test data for calculating center detonation propose new prediction accordingly for analyzing superpressure changing rule Formula.From Fig. 8 and Fig. 9 as can be seen that using detonation different draw ratios in cylinder explosive center determined by the method for the present invention Axial-radial peak overpressure be sufficiently close to analog result, illustrate the present invention determination method it is accurate and reliable.

Claims (3)

1. a kind of cylinder explosive axial-radial shock wave overpressure peak determines method, which is characterized in that include the following steps:
(10) initial parameter obtains:Obtain the initial parameter of cylinder explosive, including appearance and size, quality, explosive type, superpressure Calculate distance;
(20) intermediate parameters calculate:According to initial parameter, the intermediate parameters of cylinder explosive, including equivalent TNT matter is calculated Amount, ratio are quick-fried away from, draw ratio;
(30) peak overpressure determines:According to intermediate parameters and peak overpressure fitting formula, be calculated cylinder explosive it is axial and Radial impact wave peak overpressure.
2. peak overpressure according to claim 1 determines method, which is characterized in that (30) peak overpressure determines step Including:
(31) radial peak overpressure calculates:The radial impact wave that the detonation of cylinder explosive center and one end detonation is calculated as follows is super Voltage crest value, unit kpa,
(32) detonation axial peak overpressure in center calculates:The axial impact wave superpressure of center detonation cylinder explosive is calculated as follows Peak value, unit kpa,
(33) detonation axial distal end peak overpressure in one end calculates:One end is calculated as follows and detonates cylinder explosive axially away from detonation The peak overpressure of point one end:
PAxis, one end=1.4PAxis, center,
(34) detonation axial proximal peak overpressure in one end calculates:One end is calculated as follows and detonates cylinder explosive axially adjacent to detonation The peak overpressure of point one end:
PAxis, one end=0.5PAxis, center,
In above-mentioned formula, Z be ratio it is quick-fried away from, in draw ratio L/D, L is the height of cylinder, and D is cylinder bottom surface diameter of a circle, a, B is mediant, and e is scientific notation symbol, P0It is atmospheric pressure.
3. peak overpressure according to claim 2 determines method, which is characterized in that described ((20) intermediate parameters calculate step Suddenly include:
(21) equivalent TNT Mass Calculations:According to the corresponding TNT equivalent coefficients k of explosive type, equivalent TNT mass is calculated as follows W,
W=k × m,
In formula, m is the quality of cylinder explosive, unit kg;
(22) ratio is quick-fried away from calculating:It is quick-fried away from Z that ratio is calculated as follows,
In formula, R is distance of the pressure measurement point apart from explosive geometric center, and unit is rice.
CN201810436867.XA 2018-05-09 2018-05-09 Method for determining overpressure peak value of axial and radial shock waves of cylindrical explosive Active CN108627629B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201810436867.XA CN108627629B (en) 2018-05-09 2018-05-09 Method for determining overpressure peak value of axial and radial shock waves of cylindrical explosive

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201810436867.XA CN108627629B (en) 2018-05-09 2018-05-09 Method for determining overpressure peak value of axial and radial shock waves of cylindrical explosive

Publications (2)

Publication Number Publication Date
CN108627629A true CN108627629A (en) 2018-10-09
CN108627629B CN108627629B (en) 2020-06-19

Family

ID=63692313

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201810436867.XA Active CN108627629B (en) 2018-05-09 2018-05-09 Method for determining overpressure peak value of axial and radial shock waves of cylindrical explosive

Country Status (1)

Country Link
CN (1) CN108627629B (en)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109975142A (en) * 2019-04-30 2019-07-05 公安部第一研究所 A kind of non-contact explosion wave superpressure test macro of plate product and method

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2010068328A1 (en) * 2008-12-10 2010-06-17 Continental Disc Corporation Controlling the rated burst pressure of a rupture disc through the use of control scores on the disc dome
CN103558118A (en) * 2013-10-11 2014-02-05 西安近代化学研究所 Internal explosion experiment based test method for quantitative evaluation on thermobaric effect of thermobaric explosive
CN104568613A (en) * 2014-12-19 2015-04-29 西安近代化学研究所 Ship cabin indoor explosion structure damage evaluation method based on equivalent scale model
CN104743086A (en) * 2015-02-14 2015-07-01 中国人民解放军海军工程大学 Ship whipping motion response prediction method and system under underwater explosion situation
CN108009385A (en) * 2017-12-29 2018-05-08 中国人民解放军陆军工程大学 Method for estimating equivalent explosive amount of explosive for simulating centrifugal machine underwater explosion test
CN105389415B (en) * 2015-10-16 2018-08-10 华侨大学 A kind of prediction technique of column charge blasting vibration particle peak velocity

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2010068328A1 (en) * 2008-12-10 2010-06-17 Continental Disc Corporation Controlling the rated burst pressure of a rupture disc through the use of control scores on the disc dome
CN103558118A (en) * 2013-10-11 2014-02-05 西安近代化学研究所 Internal explosion experiment based test method for quantitative evaluation on thermobaric effect of thermobaric explosive
CN104568613A (en) * 2014-12-19 2015-04-29 西安近代化学研究所 Ship cabin indoor explosion structure damage evaluation method based on equivalent scale model
CN104743086A (en) * 2015-02-14 2015-07-01 中国人民解放军海军工程大学 Ship whipping motion response prediction method and system under underwater explosion situation
CN105389415B (en) * 2015-10-16 2018-08-10 华侨大学 A kind of prediction technique of column charge blasting vibration particle peak velocity
CN108009385A (en) * 2017-12-29 2018-05-08 中国人民解放军陆军工程大学 Method for estimating equivalent explosive amount of explosive for simulating centrifugal machine underwater explosion test

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
刘磊等: "柱形装药水下爆炸远场冲击波压力峰值分布", 《振动与冲击》 *
马珊珊等: "装药形状对冲击波超压峰值的影响", 《首届兵器工程大会论文集》 *

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109975142A (en) * 2019-04-30 2019-07-05 公安部第一研究所 A kind of non-contact explosion wave superpressure test macro of plate product and method

Also Published As

Publication number Publication date
CN108627629B (en) 2020-06-19

Similar Documents

Publication Publication Date Title
CN108984996B (en) Aluminum powder reactivity-based aluminum-containing explosive JWL state equation parameter calculation method
Cheng et al. Detonation characteristics of emulsion explosives sensitized by MgH 2
CN110823435A (en) Explosion impulse testing device, system and method
CN108627629A (en) Method for determining overpressure peak value of axial and radial shock waves of cylindrical explosive
CN115221814A (en) Method for predicting chemical explosion shock wave intensity of hydrogen pipeline
CN106248282A (en) A kind of wall pressure measurement apparatus of resistance to HI high impact
CN104215949B (en) A kind of target RCS scaling method based on AGC voltage correction
Chiquete et al. Detonation performance experiments and modeling for the DAAF-based high explosive PBX 9701
CN116183164B (en) Test system and model building method for simulating warehouse explosion
CN105626291B (en) Internal combustion engine and its control method
Hornung et al. On the question of universality of imploding shock waves
CN103743638A (en) Simulated composite loading device for explosive wave
CN107430097A (en) For calibrating the instrument of supersonic inspection device
Wang et al. Study on pressure reconstruction method of explosion shock wave
RU2739303C1 (en) Method of determining reacted explosive charge mass
Baker et al. Theory and calibration of JWL and JWLB thermodynamic equations of state
RU2566417C1 (en) Method of graduation of pressure sensors for sensing of air shock waves
Apparao et al. TNT equivalency of unconfined aerosols of propylene oxide
Fiebig Influence of the inter teeth volumes on the noise generation in external gear pumps
Elkarous et al. Investigation on gas pressure measurement inside small caliber weapons with piezoelectric transducers
Xin et al. Numerical simulation of underwater explosion loads
TRANĂ et al. Numerical simulation study on the ring fragmentation
Bai et al. Study on the shockwave pressure testing method for detonating cords
Thangadurai et al. Numerical study of blast wave generation through shock tube and its comparison with free-field blast wave parameters
RU2490589C1 (en) Method to determine initiating capacity of remote ammunition and related device

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