Explosion center position determination method and system combined with dynamic pressure

Abstract

The invention relates to the technical field of blasting, and discloses a method and a system for determining the position of a blasting center combined with dynamic pressure so as to ensure that the substantial damage effect meets the expectation. The method comprises the following steps: determining an A damage grade and a first initial explosion center position expected for an A target; the specific impulse, the dynamic pressure and the positive pressure acting time of the ammunition at the target position are calculated according to the equivalent TNT equivalent of the current ammunition and the distance between the first initial explosion center position and the A target; a damage coefficient is calculated, and the substantial damage level of a current damage coefficient calculation value is judged; and judging whether the substantial damage level is consistent with the A damage level, if so, determining the first initial blasting center position as a reasonable position, and if not, adjusting the first initial blasting center position according to a difference value between the substantial damage level and the expected damage level until the substantial damage level is consistent with the expected damage level.

Description

Technical Field

[0001] The present invention relates to the technical field of blasting, and in particular, to a method and system for determining the position of the explosion center with combined dynamic pressure. Background Art

[0002] As the most important damage element in explosion damage, the shock wave is a commonly used means for target damage. In the assessment of shock wave damage, the shock wave damage criteria based on the shock wave overpressure are mainly adopted, such as the shock wave overpressure criterion, the impulse criterion, or the overpressure-impulse criterion, without considering the influence of dynamic pressure. In the actual process, dynamic pressure may also cause structural damage, or cause damage such as the target toppling, translating, rolling, and throwing.

[0003] For example: When conducting a shock wave test on a certain oil tanker carriage, a shock wave overpressure of 0.09 MPa cannot damage the structure of the oil tanker carriage, but at this time, the rolling caused by the dynamic pressure makes the whole vehicle basically scrapped.

[0004] With the continuous emergence of high-energy explosives and the increase in the warhead equivalent, the energy released instantaneously by its explosion shock wave is higher and the dynamic pressure is greater, becoming an important damage factor for targets such as armored vehicles and buildings. If only the criterion based on the shock wave overpressure is used, the shock wave damage effect cannot be accurately evaluated. Therefore, it is of great practical significance to establish a combined damage assessment method for the peak shock wave overpressure, specific impulse, and dynamic pressure considering dynamic pressure and use it to guide the reasonable determination of the explosion center position. Summary of the Invention

[0005] The purpose of the present invention is to disclose a method and system for determining the position of the explosion center with combined dynamic pressure to ensure that the actual damage effect meets the expectations.

[0006] To achieve the above purpose, the method for determining the position of the explosion center with combined dynamic pressure disclosed by the present invention includes: Step S1, determining the expected damage level of Target A and the first initial explosion center position; Step S2, calculating according to the equivalent TNT equivalent of the current ammunition and the distance between the first initial explosion center position and Target A and the positive pressure action time of the ammunition at the target position; wherein, are respectively the peak shock wave overpressure, specific impulse, and dynamic pressure of the target; Step S3, calculating the damage coefficient , and the calculation formula is: ; wherein, are respectively the damage thresholds corresponding to the peak shock wave overpressure, specific impulse, and dynamic pressure of the target, They are the peak overpressure, specific impulse, and dynamic pressure damage coefficient of the shock wave related to the target and the positive pressure action time respectively; then determine the actual damage level where the current damage coefficient calculation value is located. Indicates no damage. Indicates mild damage. Indicates moderate damage. Is severe damage. Is a pre-calibrated constant. Step S4: Determine whether the actual damage level is consistent with the A damage level. If it is consistent, determine the first initial burst center position as the reasonable position. If it is not consistent, adjust the first initial burst center position according to the difference between the actual damage level and the expected damage level until the actual damage level is consistent with the expected damage level.

[0007] Preferably, the present invention further includes: During the process of iteratively adjusting the distance between the first initial burst center position and the target A based on step S4, determine the area range A that meets the A damage level. Then obtain the expected B damage level for the B target. Randomly select a position within the area range A as the second initial burst center position, and then determine the area range B that meets the B damage level in the same way as determining the area range A that meets the A damage level. Take the intersection of the area range A and the area range B to obtain the target area that simultaneously meets the A damage level expected for the A target and the B damage level expected for the B target.

[0008] Preferably, the present invention The calculation formula is: ; where Is the equivalent TNT equivalent. Is the distance from the burst center to the target.

[0009] Preferably, the present invention The calculation formula is: ; ; where Is the scaled distance. Is the impulse coefficient.

[0010] Preferably, the present invention The calculation formula is: ; ; ; ; ; Where , Are the air pressures in front of and behind the wave front respectively. , are the air densities before and after the wavefront respectively; is the air adiabatic index, is the dynamic pressure velocity, is the air velocity before the wavefront.

[0011] Preferably, the positive pressure time of the present invention has the following calculation formula: .

[0012] Preferably, during the calibration of the present invention according to the target and the positive pressure action time, the distribution of the value range of the constraint is as follows: I. If the target is a metal material target: When the relationship between the positive pressure action time and the vibration period T of the target structure itself satisfies / T≥10, , , The value range is: 0.8≤ ≤1.2, 0.8≤ ≤1.1, 0.8≤ ≤1.3, and it satisfies < < ; When the relationship between the positive pressure action time and the vibration period T of the target structure itself satisfies / T≤0.25, , , The value range is: 0.8≤ ≤1.1, 0.8≤ ≤1.2, 0.8≤ ≤1.3, and it satisfies < < ; When the ratio of the positive pressure action time to the vibration period T of the target structure itself is in the range of: 10≤ / T≤0.25, , , The value range is: 0.8≤ = ≤1.2, 0.8≤ ≤1.3, and it satisfies = < ; II. If the target is a reinforced concrete target: When the positive pressure acting time and the vibration period T of the target structure itself satisfy the relationship / T≥10, 、 、 the value range is: 0.8≤ ≤1.3, 0.8≤ ≤1.2, 0.8≤ ≤1.1, and satisfy < < ; When the positive pressure acting time and the vibration period T of the target structure itself satisfy the relationship / T≤0.25, 、 、 the value range is: 0.8≤ ≤1.2, 0.8≤ ≤1.3, 0.8≤ ≤1.1, and satisfy < < ; When the ratio of the positive pressure acting time to the vibration period T of the target structure itself is in the range: 10≤ / T≤0.25, 、 、 the value range is: 0.8≤ = ≤1.3, 0.8≤ ≤1.1, and satisfy < = 。

[0013] To achieve the above object, the present invention also discloses a method for determining the burst center position in combination with dynamic pressure, including a memory, a processor, and a computer program stored on the memory and executable on the processor. When the processor executes the computer program, the above method is implemented.

[0014] The present invention has the following beneficial effects: Simple and practical, comprehensively considering the influence of the peak overpressure of the shock wave, specific impulse, and dynamic pressure on the damage of the target, and can effectively ensure that the actual damage effect meets the expectations. Among them, in the present invention, 、 、 The determination of parameters such as etc. can be carried out based on the measured test data combined with the finite element simulation. In the specific implementation process, each ammunition enterprise can formulate an exclusive enterprise standard according to the simulation and measured results to guide downstream customers to accurately place the explosion center position during use.

[0015] The present invention will be further described in detail below with reference to the accompanying drawings. BRIEF DESCRIPTION OF THE DRAWINGS

[0016] The accompanying drawings that form a part of this application are used to provide a further understanding of the present invention. The schematic embodiments of the present invention and their descriptions are used to explain the present invention and do not constitute an improper limitation of the present invention. In the drawings: Figure 1 It is a schematic flow chart of a method for determining the explosion center position of combined dynamic pressure disclosed in an embodiment of the present invention. DETAILED DESCRIPTION OF THE EMBODIMENTS

[0017] The embodiments of the present invention will be described in detail below with reference to the accompanying drawings. However, the present invention can be implemented in many different ways defined and covered by the claims.

[0018] Embodiment 1 This embodiment discloses a method for determining the explosion center position of combined dynamic pressure, as Figure 1 shown, including: Step S1: Determine the expected damage level A and the first initial explosion center position for target A.

[0019] Step S2: Calculate the equivalent TNT equivalent of the current ammunition, the distance between the first initial explosion center position and target A, specific impulse, dynamic pressure, and the positive pressure action time of the ammunition at the target position.

[0020] In this step, the peak overpressure of the target shock wave The calculation formula can be: ; where is the equivalent TNT equivalent; is the distance from the explosion center to the target.

[0021] The specific impulse The calculation formula can be: ; ; where is the scaled distance; is the impulse coefficient.

[0022] The dynamic pressure The calculation formula can be: ; ; ; ; .

[0023] Among them, and are the air pressures before and after the wavefront respectively, and are the air densities before and after the wavefront respectively; is the air adiabatic index, is the dynamic pressure velocity, is the air velocity before the wavefront.

[0024] In this step, the positive pressure time can be calculated by the formula: .

[0025] Step S3: Calculate the damage coefficient and determine the actual damage level where the current damage coefficient calculation value is located.

[0026] In this step, the damage coefficient is calculated by the formula: .

[0027] Among them, are the damage thresholds corresponding to the peak overpressure, specific impulse, and dynamic pressure of the target shock wave respectively, are the damage coefficients of the peak overpressure, specific impulse, and dynamic pressure of the shock wave related to the target and the positive pressure action time respectively.

[0028] In this embodiment, is pre-calibrated according to the target and the positive pressure action time. Generally, during the calibration process, the distribution of the value range of the constraint is: I. If the target is a metal material target such as a weapon or technical equipment: When the positive pressure action time is much greater than the vibration period T of the target structure itself, that is, / T≥10, , , The value range is: 0.8≤ ≤1.2, 0.8≤ ≤1.1, 0.8≤ ≤1.3, and it satisfies < < .

[0029] When the positive pressure action time is much less than the vibration period T of the target structure itself, that is, / T≤0.25, , , The value range is: 0.8≤ ≤1.1, 0.8≤ ≤1.2, 0.8≤ ≤1.3, and satisfy < < 。

[0030] When the positive pressure acting time The ratio to the vibration period T of the target structure itself is in the range: 10≤ / T≤0.25, 、 、 The value range is: 0.8≤ = ≤1.2, 0.8≤ ≤1.3, and satisfy = < 。

[0031] II. If the target is a reinforced concrete target such as an open works, a shelter, a command center, etc.: When the positive pressure acting time Is much greater than the vibration period T of the target structure itself, that is / T≥10, 、 、 The value range is: 0.8≤ ≤1.3, 0.8≤ ≤1.2, 0.8≤ ≤1.1, and satisfy < < 。

[0032] When the positive pressure acting time Is much less than the vibration period T of the target structure itself, that is / T≤0.25, 、 、 The value range is: 0.8≤ ≤1.2, 0.8≤ ≤1.3, 0.8≤ ≤1.1, and satisfy < < 。

[0033] When the positive pressure acting time The ratio to the vibration period T of the target structure itself is in the range: 10≤ / T≤0.25, 、 、 The value range is: 0.8 ≤ = ≤1.3,0.8≤ ≤1.1, and meets < = .

[0034] In this step, preferably, the actual damage level of the current damage coefficient calculation value can be determined as follows: Indicates no damage, Indicates minor injury; Indicates moderate damage. For severe injuries; is a pre-calibrated constant.

[0035] Step S4, determine whether the actual damage level is consistent with the A damage level. If consistent, determine the first initial explosion center position as a reasonable position. If inconsistent, adjust the first initial explosion center position according to the difference between the actual damage level and the expected damage level until the actual damage level is consistent with the expected damage level.

[0036] In this step, usually, for the same ammunition, the closer the distance between the explosion center and the target, the greater the damage coefficient value, so that the area range A that meets the A damage level can be determined in the process of iteratively adjusting the distance between the first initial explosion center position and the target A. Furthermore, the method of this embodiment can also obtain the B damage level expected for the B target, randomly select a position in the area range A as the second initial explosion center position, and then determine the area range B that meets the B damage level in the same way as determining the area range A that meets the A damage level. Finally, the intersection of the area range A and the area range B is taken to obtain a target area that meets both the A damage level expected for the A target and the B damage level expected for the B target.

[0037] Similarly, if there are more damage levels specified for different targets, the above method can be used and so on, and no further explanation is given. When the intersection is empty, an error message can be output for the user to adjust.

[0038] Example 2 This embodiment discloses a method for determining the explosion center position by combining dynamic pressure, including a memory, a processor, and a computer program stored in the memory and executable on the processor. When the processor executes the computer program, the above method is implemented.

[0039] In summary, the method and system disclosed in the embodiments of the present invention have the following beneficial effects: It is simple and practical, and comprehensively considers the impact of shock wave overpressure peak value, specific impulse and dynamic pressure on target damage, which can effectively ensure that the actual damage effect meets expectations. , , The determination of parameters such as can be based on the combination of measured test data and finite element simulation. In the specific implementation process, each ammunition enterprise can formulate exclusive enterprise standards according to the simulation and measurement results to guide downstream customers to accurately place the explosion center during use.

[0040] The above are only the preferred embodiments of the present invention and are not intended to limit the present invention. For those skilled in the art, the present invention can have various changes and modifications. Any modification, equivalent replacement, improvement, etc. made within the spirit and principle of the present invention shall be included within the protection scope of the present invention.

Claims

1. A method for determining the position of the explosion center by combining dynamic pressure, characterized in that, Including: Step S1: Determine the A damage level expected for target A and the first initial burst center position; Step S2: Calculate based on the equivalent TNT equivalent of the current ammunition and the distance between the first initial explosion center position and the A target and the positive pressure action time of the ammunition at the target position; where are the peak overpressure, specific impulse, and dynamic pressure of the target shock wave respectively; Step S3, calculate the damage coefficient , and the calculation formula is: ； wherein, are the damage thresholds corresponding to the peak overpressure, specific impulse, and dynamic pressure of the target shock wave respectively, are the overpressure peak, specific impulse, and dynamic pressure damage coefficients of the shock wave related to the target and the positive pressure action time respectively; then determine the actual damage level where the current calculated value of the damage coefficient is located, indicates no damage, indicates minor damage; indicates moderate damage, is severe damage; is a pre-calibrated constant; Step S4: Judge whether the actual damage level is consistent with the A damage level. If they are consistent, determine the first initial burst center position as the reasonable position. If they are not consistent, adjust the first initial burst center position according to the difference between the actual damage level and the expected damage level until the actual damage level is consistent with the expected damage level.

2. The method for determining the position of the explosion center of combined dynamic pressure according to claim 1, characterized in that It also includes: During the process of iteratively adjusting the distance between the first initial burst center position and target A based on Step S4, determine the area range A that meets the A damage level; Then obtain the B damage level expected for target B. Randomly select a position within the area range A as the second initial burst center position, and then determine the area range B that meets the B damage level in the same way as determining the area range A that meets the A damage level; Take the intersection of the area range A and the area range B to obtain the target area that simultaneously meets the A damage level expected for target A and the B damage level expected for target B.

3. The method for determining the position of the explosion center of combined dynamic pressure according to claim 1 or 2, characterized in that, The calculation formula is as follows: wherein, is the equivalent TNT equivalent; is the distance from the explosion center to the target.

4. The method for determining the position of the explosion center of combined dynamic pressure according to claim 3, characterized in that, The calculation formula is: ; ; where is the proportional distance; is the impulse coefficient.

5. The method for determining the position of the explosion center of combined dynamic pressure according to claim 4, wherein The calculation formula is as follows: ； ； ； ； ； Among them, and are the air pressures before and after the wave front respectively, and are the air densities before and after the wave front respectively; is the air adiabatic index, is the dynamic pressure velocity, is the air velocity before the wave front.

6. The method for determining the position of the explosion center of combined dynamic pressure according to claim 5, characterized in that Positive pressure time The calculation formula is as follows: .

7. The method for determining the position of the explosion center of combined dynamic pressure according to claim 6, characterized in that During the calibration according to the target and the positive pressure application time, the distribution of the value range of the constraint is as follows:

1. If the target is a metal material target: When the positive pressure acting time and the relationship with the vibration period T of the target structure itself satisfies / T≥10, 、 、 The value range is: 0.8≤ ≤1.2, 0.8≤ ≤1.1, 0.8≤ ≤1.3, and satisfies < < ; When the positive pressure acting time and the relationship with the vibration period T of the target structure itself satisfies / T ≤ 0.25, 、 、 The value range is: 0.8 ≤ ≤ 1.1, 0.8 ≤ ≤ 1.2, 0.8 ≤ ≤ 1.3, and satisfies < < ; When the ratio of the positive pressure acting time to the vibration period T of the target structure itself is in the range of: 10 ≤ / T ≤ 0.25, , , The value range is: 0.8 ≤ = ≤ 1.2, 0.8 ≤ ≤ 1.3, and it satisfies = < ; 2. If the target is a reinforced concrete target: When the positive pressure acting time and the relationship between the vibration period T of the target structure itself satisfies / T≥10, 、 、 The value range is: 0.8≤ ≤1.3, 0.8≤ ≤1.2, 0.8≤ ≤1.1, and satisfies < < ; When the positive pressure acting time and the vibration period T of the target structure itself satisfy / T ≤ 0.25, 、 、 The value range is: 0.8 ≤ ≤ 1.2, 0.8 ≤ ≤ 1.3, 0.8 ≤ ≤ 1.1, and satisfy < < ; When the ratio of the positive pressure acting time to the vibration period T of the target structure itself is in the range: 10 ≤ / T ≤ 0.25, 、 、 The value range is: 0.8 ≤ = ≤ 1.3, 0.8 ≤ ≤ 1.1, and satisfy < = .

8. A method for determining the position of the explosion center of combined dynamic pressure, including a memory, a processor, and a computer program stored on the memory and operable on the processor, characterized in that, When the processor executes the computer program, it implements the method described in any one of claims 1 to 7 above.

Images