CN110866296A - Method for rapidly evaluating damage capability of explosive shock waves to fixed target - Google Patents
Method for rapidly evaluating damage capability of explosive shock waves to fixed target Download PDFInfo
- Publication number
- CN110866296A CN110866296A CN201910929225.8A CN201910929225A CN110866296A CN 110866296 A CN110866296 A CN 110866296A CN 201910929225 A CN201910929225 A CN 201910929225A CN 110866296 A CN110866296 A CN 110866296A
- Authority
- CN
- China
- Prior art keywords
- volume
- fixed target
- shock wave
- intersecting
- damage
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
- 230000035939 shock Effects 0.000 title claims abstract description 48
- 239000002360 explosive Substances 0.000 title claims abstract description 27
- 238000000034 method Methods 0.000 title claims abstract description 18
- 238000004364 calculation method Methods 0.000 claims abstract description 10
- 238000004880 explosion Methods 0.000 claims abstract description 6
- 238000011156 evaluation Methods 0.000 claims description 10
- 238000000342 Monte Carlo simulation Methods 0.000 claims description 7
- 239000000126 substance Substances 0.000 claims description 4
- 239000004576 sand Substances 0.000 claims description 2
- 230000000694 effects Effects 0.000 description 3
- 230000007547 defect Effects 0.000 description 2
- 239000000463 material Substances 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000005422 blasting Methods 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
Images
Landscapes
- Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)
Abstract
The invention discloses a method for rapidly evaluating the damage capability of an explosive shock wave to a fixed target, which comprises the following steps: s1, initializing the parameters; s2, dispersing the effective volume formed by the blast shock wave, and determining the intersecting volume V of the effective volume and the fixed targetv‑s(ii) a S3, dispersing the fixed target involved by the blast shock wave, and determining the intersecting volume V of the fixed target and the effective volumes‑v(ii) a S4, according to the intersecting volume Vv‑sIntersecting volume Vs‑vAnd determining a final intersecting volume V, and evaluating the damage capability E of the explosive shock wave to the fixed target according to the final intersecting volume V. The damage capability of the explosion shock wave obtained by calculation to the fixed target can meet the requirement of any precision, and the arbitrarily intersected volume of the sphere and the hexahedron can be quickly given.
Description
Technical Field
The invention belongs to the field of blast shock wave damage effect evaluation, and particularly relates to a rapid evaluation method for the damage capability of blast shock waves to a fixed target.
Background
In the field of evaluation of the damage effect of the explosive shock wave, the damage capability of the explosive shock wave to a target is scientifically evaluated, and the method has important engineering application value for development of various blasting warheads, power evaluation and performance evaluation of equipped killer weapons. In addition, research on blast shock wave damage effect evaluation provides scientific reference for structural design and safety protection of key buildings, bridges and key targets.
The explosion shock wave damage capability is mainly characterized in that an experiment is carried out to obtain a shock wave overpressure value, then theoretical modeling is carried out, finally a damage criterion is obtained, and further the damage capability is calculated. However, this method is very labor and material intensive, and the experimental conditions are complicated.
Disclosure of Invention
In order to calculate the damage capability of the explosive shock wave to the fixed target, the invention aims to provide a method for rapidly evaluating the damage capability of the explosive shock wave to the fixed target, and aims to solve the problems that the existing method consumes huge manpower and material resources, is complex in experimental conditions and the like.
The invention is realized in such a way that a rapid evaluation method of the damage capability of the explosive shock wave to a fixed target comprises the following steps:
s1, initializing parameters including the effective damage radius R of the damage of the given explosion shock wave, the length L, the width W, the height H and the calculation precision epsilon of the target parameters;
s2, dispersing the effective volume formed by the blast shock wave, and determining the intersecting volume V of the effective volume and the fixed targetv-s;
S3, dispersing the fixed target involved by the blast shock wave, and determining the intersecting volume V of the fixed target and the effective volumes-v;
S4, according to the intersecting volume Vv-sIntersecting volume Vs-vAnd determining a final intersecting volume V, and evaluating the damage capability E of the explosive shock wave to the fixed target according to the final intersecting volume V.
Preferably, in step S2, M ═ 1/epsilon points are randomly generated by the monte carlo methodDispersing the effective volume formed by the blast shock wave; where v denotes a random point at which an effective volume formed by an explosive shock wave is discretely generated, and i denotes an ith random point i ═ 1, 2.
Preferably, in step S2, the intersecting volume Vv-sComprises the following steps:
wherein the content of the first and second substances,to calculate M pointsBoth in the effective volume of the blast shock wave and in the number of fixed target points.
Preferably, in step S3, M random points are randomly generated by the monte carlo methodDiscretizing the fixed target; where s is a random point generated by discretizing a fixed target, and j is a random point generated by discretizing a fixed target.
Preferably, in step S3, the intersecting volume Vs-vComprises the following steps:
wherein the content of the first and second substances,to calculate M pointsIn the hexahedron and in the middle of the sphere.
Preferably, in step S4, the final intersecting volume V is the intersecting volume Vv-sAnd the intersecting volume Vs-vThe average of the sums; the damage capability of the explosive shock wave to the fixed target is as follows:
wherein, VtIs the volume of the fixation target.
The invention overcomes the defects of the prior art and provides a method for rapidly evaluating the damage capability of an explosive shock wave to a fixed target. It is clear that the calculated volume is closer to the true value as the monte carlo method generates points faster and more evenly distributed. The intersection volume obtained by calculation of the invention not only can meet the requirement of any precision, but also can quickly give the intersection volume of the sphere and the hexahedron at will.
Compared with the defects and shortcomings of the prior art, the invention has the following beneficial effects:
(1) the damage capability of the explosion shock wave obtained by calculation on the fixed target can meet the requirement of any precision;
(2) the invention can quickly give the arbitrarily intersected volume of the sphere and the hexahedron.
Drawings
FIG. 1 is a flow chart of an evaluation method of the present invention;
FIG. 2 illustrates the intersection of a sphere and a hexahedron in an embodiment of the present invention;
FIG. 3 illustrates the variation of the intersection volume calculation with respect to discrete points in an embodiment of the present invention;
FIG. 4 shows the absolute error and the relative error of the calculation result of the intersected volume in the embodiment of the invention and the variation relation of discrete points.
Detailed Description
In order to make the objects, technical solutions and advantages of the present invention more apparent, the present invention is described in further detail below with reference to the accompanying drawings and embodiments. It should be understood that the specific embodiments described herein are merely illustrative of the invention and are not intended to limit the invention.
A method for rapidly evaluating the damage capability of an explosive shock wave to a fixed target, as shown in fig. 1, the method comprises the following steps:
s1, parameter initialization
And giving an effective damage radius R of explosion shock wave damage, the length L, the width W and the height H of fixed target parameters and the calculation precision epsilon.
S2, dispersing the effective volume formed by the blast shock wave, and determining the intersecting volume V of the effective volume and the fixed targetv-s
whereinCenter coordinates of an effective volume formed for the blast shock wave; v denotes a random point at which an effective volume formed by the blast shock wave is discretely generated, and i denotes an ith random point i ═ 1, 2.
According to step S2, M points are calculatedIn both the effective volume of the blast shock wave and the number of fixed target midpointsI.e. calculating M pointsA number satisfying the following formula:
then the intersecting volume V at this timev-sComprises the following steps:
s3, dispersing the fixed target involved by the blast shock wave, and determining the intersecting volume V of the fixed target and the effective volumes-v
In step S3, M points are randomly generated using the Monte Carlo methodThe fixed target is discretized and these points satisfy:
wherein (x)o,yo,zo) Is the center coordinate of the fixed target; s is the random point generated for the fixed target discretization, and j is the random point generated for the fixed target discretization.
According to step S3, M points are calculatedIn the hexahedron and the middle point of the sphereI.e. calculating M pointsA number satisfying the following formula:
then the intersecting volume V at this times-vComprises the following steps:
s4, according to the intersecting volume Vv-sIntersecting volume Vs-vDetermining the final intersection volume VtAnd evaluating the damage capability E of the explosive shock wave to the fixed target according to the final intersecting volume V
According to steps S2 and S3, the final intersection volume V of the sphere and the hexahedron is calculated:
at this time, the damage capability of the explosive shock wave to the fixed target is as follows:
E=V/Vt=(Vv-s+Vs-v)/2Vt(8)
wherein, VtIs the volume of the fixation target.
Example 2
Assuming that the effective acting volume of the blast shock wave is a sphere with a radius R equal to 2 and a center coordinate of the circleThe function of the sphere is then:
further assume that the fixed target is a hexahedron, whose parameters are:
and drawing the intersection graph of the sphere and the hexahedron by using Matlab software, as shown in FIG. 2. It can be seen that the sphere is completely enclosed in the hexahedron, and the intersecting volume is the volume of the sphere, and the calculation formula according to the volume of the sphere is:
the method of the invention is used to calculate the intersection volume, and the absolute error and the relative error are respectively used to calculate the precision of the method provided by the invention:
here, when the intersecting volume is calculated by the present invention, the accuracy range is set to [1/80000,1/10000], the number of discrete points M is [10000,80000], the sphere and the hexahedron are discretized, and finally the relationship between the calculation result of the intersecting volume and the number of discrete points is shown in fig. 3, and the corresponding absolute error and relative error are shown in fig. 4.
As can be seen from fig. 3 and 4, when the number of discrete points M is 50,000, the intersecting volume is Ve33.4780, the corresponding error is
At this time, the estimation result has a small error and is closer to the real intersection volume.
Further calculating the damage capability of the explosive shock wave to the fixed target when M is 50,000
Obviously, the invention can rapidly calculate the damage capability of the explosive shock wave to the fixed target.
The above description is only for the purpose of illustrating the preferred embodiments of the present invention and is not to be construed as limiting the invention, and any modifications, equivalents and improvements made within the spirit and principle of the present invention are intended to be included within the scope of the present invention.
Claims (6)
1. A method for rapidly evaluating the damage capability of an explosive shock wave to a fixed target is characterized by comprising the following steps:
s1, initializing parameters including the effective damage radius R of the damage of the given explosion shock wave, the length L, the width W, the height H and the calculation precision epsilon of the target parameters;
s2, dispersing the effective volume formed by the blast shock wave, and determining the intersecting volume V of the effective volume and the fixed targetv-s;
S3, dispersing the fixed target involved by the blast shock wave, and determining the intersecting volume V of the fixed target and the effective volumes-v;
S4, according to the intersecting volume Vv-sIntersecting volume Vs-vAnd determining a final intersecting volume V, and evaluating the damage capability E of the explosive shock wave to the fixed target according to the final intersecting volume V.
2. The method for rapidly evaluating the destructibility of a fixed object by an explosive blast as defined in claim 1, wherein M ═ 1/epsilon points are randomly generated by the monte carlo method in step S2Dispersing the effective volume formed by the blast shock wave; where v denotes a random point at which an effective volume formed by an explosive shock wave is discretely generated, and i denotes an ith random point i ═ 1, 2.
3. The method for rapid assessment of the ability of an explosive blast to damage a fixed target according to claim 2, wherein in step S2, said intersecting volume Vv-sComprises the following steps:
4. The method of claim 1 for rapid evaluation of the ability of an explosive blast to damage a fixed targetWherein, in step S3, M random points are randomly generated by the Monte Carlo methodDiscretizing the fixed target; where s is a random point generated by discretizing a fixed target, and j is a random point generated by discretizing a fixed target.
5. The method for rapid assessment of the ability of an explosive blast to damage a fixed target according to claim 4, wherein in step S3, said intersecting volume Vs-vComprises the following steps:
6. The method for rapid evaluation of the ability of an explosive blast to damage a fixed target according to claim 1, wherein in step S4, said final intersecting volume V is an intersecting volume Vv-sAnd the intersecting volume Vs-vThe average of the sums; the damage capability of the explosive shock wave to the fixed target is as follows:
wherein, VtIs the volume of the fixation target.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910929225.8A CN110866296A (en) | 2019-09-29 | 2019-09-29 | Method for rapidly evaluating damage capability of explosive shock waves to fixed target |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910929225.8A CN110866296A (en) | 2019-09-29 | 2019-09-29 | Method for rapidly evaluating damage capability of explosive shock waves to fixed target |
Publications (1)
Publication Number | Publication Date |
---|---|
CN110866296A true CN110866296A (en) | 2020-03-06 |
Family
ID=69652862
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910929225.8A Pending CN110866296A (en) | 2019-09-29 | 2019-09-29 | Method for rapidly evaluating damage capability of explosive shock waves to fixed target |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110866296A (en) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN116861712A (en) * | 2023-09-05 | 2023-10-10 | 中国兵器装备集团兵器装备研究所 | Efficient damage modeling method and system based on few public parameters |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103177191A (en) * | 2013-04-11 | 2013-06-26 | 北京理工大学 | Warhead damage performance evaluation method |
KR20160114963A (en) * | 2015-03-25 | 2016-10-06 | 삼성중공업 주식회사 | Methods for modelling explosion pressure time history model |
US20180239848A1 (en) * | 2017-02-21 | 2018-08-23 | Livermore Software Technology Corporation | Numerical Blast Simulation Methods and Systems Thereof |
CN110119580A (en) * | 2019-05-17 | 2019-08-13 | 中国人民解放军军事科学院国防工程研究院 | A kind of ground surface works target Damage assessment system and method |
-
2019
- 2019-09-29 CN CN201910929225.8A patent/CN110866296A/en active Pending
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103177191A (en) * | 2013-04-11 | 2013-06-26 | 北京理工大学 | Warhead damage performance evaluation method |
KR20160114963A (en) * | 2015-03-25 | 2016-10-06 | 삼성중공업 주식회사 | Methods for modelling explosion pressure time history model |
US20180239848A1 (en) * | 2017-02-21 | 2018-08-23 | Livermore Software Technology Corporation | Numerical Blast Simulation Methods and Systems Thereof |
CN110119580A (en) * | 2019-05-17 | 2019-08-13 | 中国人民解放军军事科学院国防工程研究院 | A kind of ground surface works target Damage assessment system and method |
Non-Patent Citations (1)
Title |
---|
李青等: "催泪子母弹对体目标作用效能研究", 《火力与指挥控制》 * |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN116861712A (en) * | 2023-09-05 | 2023-10-10 | 中国兵器装备集团兵器装备研究所 | Efficient damage modeling method and system based on few public parameters |
CN116861712B (en) * | 2023-09-05 | 2023-12-22 | 中国兵器装备集团兵器装备研究所 | Efficient damage modeling method and system based on few public parameters |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Randles et al. | Smoothed particle hydrodynamics: some recent improvements and applications | |
Owen et al. | The modelling of multi‐fracturing solids and particulate media | |
Randles et al. | Normalized SPH with stress points | |
Wang et al. | A computational framework for the simulation of high‐speed multi‐material fluid–structure interaction problems with dynamic fracture | |
CN106599382A (en) | Stress solution method based on force boundary and balance conditions | |
AU2007208111B2 (en) | Sector meshing and neighbor searching for object interaction simulation | |
Benedetti et al. | A fast dual boundary element method for 3D anisotropic crack problems | |
CN110866296A (en) | Method for rapidly evaluating damage capability of explosive shock waves to fixed target | |
AU2007208110B2 (en) | Object discretization to particles for computer simulation and analysis | |
Warren | Simulations of the penetration of limestone targets by ogive-nose 4340 steel projectiles | |
Wu et al. | Debris cloud structure and hazardous fragments distribution under hypervelocity yaw impact | |
Zukas et al. | Ricochet of deforming projectiles from deforming plates | |
Banerjee et al. | Simulation of impact and fragmentation with the material point method | |
Sun et al. | Numerical analysis of the trajectory stability and penetration ability of different lateral-abnormal projectiles for non-normal penetration into soil based on modified integrated force law method | |
Zhao et al. | Discrete element modelling of dynamic behaviour of rockfills for resisting high speed projectile penetration | |
Zhdan et al. | Drag coefficients for bodies of meteorite-like shapes | |
CN113408168B (en) | High-precision numerical simulation method based on Riemann problem accurate solution | |
Zanotti | Hypervelocity impacts on planetary bodies: modelling craters formation and ejecta plume evolution | |
Ghajar et al. | Numerical analysis of hydrodynamic interaction between the linear waves of Caspian Sea and Amirkabir semi-submersible drilling platform | |
KR101456145B1 (en) | A damage assessment method for personnel and materiel in or around tunnel subjected to blast impact | |
Wilde et al. | A 3-D Dual BEM formulation for the analysis of crack growth | |
Nordendale et al. | Modeling of high rate ballistic impact of brittle armors with abaqus explicit | |
Saverkin | Comparative analysis of digital radar data processing algorithms | |
Yao | SDOT: An explicit mesh-free detonation tracking package in 2D and 3D | |
RU2611892C1 (en) | Method of three-dimensional simulation of specified hydrogeological feature implemented in computer system |
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 | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20200306 |
|
RJ01 | Rejection of invention patent application after publication |