Method and system for determining explosion center location using combined dynamic pressure

Abstract

The present invention relates to the field of blasting technology, and discloses a method and system for determining the center of gravity position in combination with dynamic pressure to ensure that the actual damage effect meets expectations. The method includes: determining the expected damage level A and the first initial center of gravity position facing target A; calculating the damage coefficient based on the equivalent TNT equivalent of the current ammunition, the distance between the first initial center of gravity position and target A, the specific impulse, the dynamic pressure, and the positive pressure action time of the ammunition at the target position; calculating the damage coefficient and determining the actual damage level of the current damage coefficient calculation value; judging whether the actual damage level is consistent with the damage level A; if so, determining the first initial center of gravity position as a reasonable position; if not, adjusting the first initial center of gravity 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.

Description

Technical Field

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

[0002] Shock waves, the primary damaging element in blast damage, are a common method of target destruction. Shock wave damage assessments primarily use shock wave overpressure-based damage criteria, such as the shock wave overpressure criterion, impulse criterion, or overpressure-impulse criterion, without considering the impact of dynamic pressure. In practice, dynamic pressure can also cause structural damage or lead to target toppling, translation, tumbling, and other damage.

[0003] For example, when a tank car is subjected to a shock wave test, the shock wave overpressure of 0.09 MPa cannot damage the tank car structure, but the rolling caused by the dynamic pressure makes the entire car basically scrapped.

[0004] With the increasing emergence of high-energy explosives and the increase in warhead yield, the instantaneous release of shock waves from explosives increases energy and dynamic pressure, making them a significant destructive factor for targets such as armored vehicles and buildings. Using only shock wave overpressure as a benchmark cannot accurately assess shock wave damage. Therefore, establishing a combined damage assessment method for shock wave overpressure peak, specific impulse, and dynamic pressure that considers dynamic pressure and uses it to guide the rational determination of the explosion epicenter is of great practical significance. Summary of the Invention

[0005] The present invention aims to disclose a method and system for determining the explosion center position by combining dynamic pressure to ensure that the actual damage effect meets expectations.

[0006] To achieve the above-mentioned purpose, the method for determining the explosion center position by combining dynamic pressure disclosed in the present invention includes:

[0007] Step S1: Determine the expected damage level A and the first initial explosion center position toward target A;

[0008] Step S2: Calculate the distance between the first initial explosion center position and the target A according to the equivalent TNT equivalent of the current ammunition and the positive pressure action time of the ammunition at the target location; among them, are the target shock wave overpressure peak, specific impulse and dynamic pressure respectively;

[0009] Step S3: Calculate the damage coefficient , the calculation formula is:

[0010] ;

[0011] in, are the damage thresholds corresponding to the target shock wave overpressure peak, specific impulse and dynamic pressure, respectively. The shock wave overpressure peak value, specific impulse and dynamic pressure damage coefficient related to the target and the positive pressure action time are respectively; then the actual damage level of the current damage coefficient calculation value is determined, Indicates no damage. Indicates minor injury; Indicates moderate damage. Severe injury; is a pre-calibrated constant;

[0012] Step S4: Determine whether the actual damage level is consistent with the damage level A. If so, determine the first initial ground center position as a reasonable position. If not, adjust the first initial ground 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.

[0013] Preferably, the present invention further comprises:

[0014] In the process of iteratively adjusting the distance between the first initial ground zero position and the target A based on step S4, determining an area range A that meets the damage level A;

[0015] Then, the expected damage level B for target B is obtained, a position is randomly selected within the area range A as the second initial ground zero position, and then the area range B meeting the damage level B is determined in the same manner as the area range A meeting the damage level A.

[0016] The intersection of the area range A and the area range B is taken to obtain a target area that simultaneously satisfies the damage level A expected for the A target and the damage level B expected for the B target.

[0017] Preferably, the present invention The calculation formula is: ;in, is the equivalent TNT equivalent; is the distance from the explosion center to the target.

[0018] Preferably, the present invention The calculation formula is: ; ;in, is the proportional distance; is the impulse coefficient.

[0019] Preferably, the present invention The calculation formula is:

[0020] ; ;

[0021] ; ; ;

[0022] in, 、 are the air pressures before and after the wave front, 、 are the air densities before and after the wave front, respectively; is the air adiabatic index, is the dynamic pressure velocity, is the air speed in front of the wave front.

[0023] Preferably, the positive pressure time of the present invention is The calculation formula is: .

[0024] Preferably, the present invention During the calibration process based on the target and the positive pressure duration, the distribution of the constraint value range is as follows:

[0025] 1. If the target is a metal material target:

[0026] When positive pressure acts for a certain time 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 meet < < ;

[0027] When positive pressure acts for a certain time 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 meet < < ;

[0028] When positive pressure acts for a certain time The ratio of the target structure's vibration period T is: 10≤ / T≤0.25, 、 、 The value range is: 0.8≤ = ≤1.2,0.8≤ ≤1.3, and meet = < ;

[0029] 2. If the target is a reinforced concrete target:

[0030] When positive pressure acts for a certain time The relationship with 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 meet < < ;

[0031] When positive pressure acts for a certain time The relationship with the vibration period T of the target structure itself satisfies / T≤0.25, 、 、 The value range is: 0.8≤ ≤1.2,0.8≤ ≤1.3,0.8≤ ≤1.1, and meet < < ;

[0032] When positive pressure acts for a certain time The ratio of the target structure's vibration period T is: 10≤ / T≤0.25, 、 、 The value range is: 0.8≤ = ≤1.3,0.8≤ ≤1.1, and meet < = .

[0033] To achieve the above objectives, the present invention also discloses a combined dynamic pressure explosion center location determination system, comprising a memory, a processor, and a computer program stored in the memory and executable on the processor, wherein the processor implements the above method when executing the computer program.

[0034] The present invention has the following beneficial effects:

[0035] It is simple and practical, and comprehensively considers the impact of shock wave overpressure peak, 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 the actual test data can be carried out based on the finite element simulation method. In the specific implementation process, each ammunition company can formulate its own corporate standards based on the simulation and actual test results to guide downstream customers to accurately deliver the explosive center during use.

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

[0037] The accompanying drawings, which constitute part of this application, are intended to provide a further understanding of the present invention. The exemplary embodiments of the present invention and their descriptions are intended to explain the present invention and do not constitute an undue limitation of the present invention. In the accompanying drawings:

[0038] Figure 1 It is a flow chart of a method for determining the explosion center position by combining dynamic pressure according to an embodiment of the present invention. DETAILED DESCRIPTION

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

[0040] Example 1

[0041] This embodiment discloses a method for determining the explosion center position by combining dynamic pressure. Figure 1 Shown, including:

[0042] Step S1: Determine the expected damage level A and the first initial explosion center position for target A.

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

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

[0045] Specific impulse The calculation formula can be: ; ;in, is the proportional distance; is the impulse coefficient.

[0046] Dynamic pressure The calculation formula can be:

[0047] ; ;

[0048] ; ; .

[0049] in, 、 are the air pressures before and after the wave front, 、 are the air densities before and after the wave front, respectively; is the air adiabatic index, is the dynamic pressure velocity, is the air speed in front of the wave front.

[0050] In this step, the positive pressure time The calculation formula can be: .

[0051] Step S3: Calculate the damage coefficient and determine the actual damage level of the current calculated damage coefficient value.

[0052] In this step, the damage factor The calculation formula is:

[0053] .

[0054] in, are the damage thresholds corresponding to the target shock wave overpressure peak, specific impulse and dynamic pressure, respectively. They are respectively the shock wave overpressure peak value, specific impulse and dynamic pressure damage coefficient related to the target and positive pressure action time.

[0055] In this embodiment, Calibration is performed in advance according to the target and the duration of positive pressure. Usually, during the calibration process, the distribution of the constraint value range is as follows:

[0056] 1. If the target is a weapon, technical equipment or other metal material target:

[0057] When positive pressure acts for a certain time is much larger 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 meet < < .

[0058] When positive pressure acts for a certain time is much smaller 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 meet < < .

[0059] When positive pressure acts for a certain time The ratio of the target structure's vibration period T is: 10≤ / T≤0.25, 、 、 The value range is: 0.8≤ = ≤1.2,0.8≤ ≤1.3, and meet = < .

[0060] 2. If the target is a reinforced concrete structure such as an open-air fortification, a sheltered fortification, or a command center:

[0061] When positive pressure acts for a certain time is much larger 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 meet < < .

[0062] When positive pressure acts for a certain time is much smaller 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 meet < < .

[0063] When positive pressure acts for a certain time The ratio of the target structure's vibration period T is: 10≤ / T≤0.25, 、 、 The value range is: 0.8≤ = ≤1.3,0.8≤ ≤1.1, and meet < = .

[0064] 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. Severe injury; is a pre-calibrated constant.

[0065] Step S4: Determine whether the actual damage level is consistent with the damage level A. If so, determine the first initial ground zero position as a reasonable position. If not, adjust the first initial ground zero 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.

[0066] In this step, for the same ammunition, the closer the distance between the ground zero and the target, the greater the damage coefficient. Therefore, during the iterative adjustment of the distance between the first initial ground zero position and target A, area range A meeting damage level A can be determined. Furthermore, the method of this embodiment can also obtain damage level B expected for target B, randomly select a location within area range A as the second initial ground zero position, and then determine area range B meeting damage level B using the same method used to determine area range A meeting damage level A. Finally, the intersection of area ranges A and B is taken to obtain a target area that meets both damage level A expected for target A and damage level B expected for target B.

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

[0068] Example 2

[0069] This embodiment discloses a system 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.

[0070] In summary, the method and system disclosed in the embodiments of the present invention have the following beneficial effects:

[0071] It is simple and practical, and comprehensively considers the impact of shock wave overpressure peak, 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 the actual test data can be carried out based on the finite element simulation method. In the specific implementation process, each ammunition company can formulate its own corporate standards based on the simulation and actual test results to guide downstream customers to accurately deliver the explosive center during use.

[0072] The foregoing description is merely a preferred embodiment of the present invention and is not intended to limit the present invention. Those skilled in the art will readily appreciate that various modifications and variations of the present invention are possible. Any modifications, equivalent substitutions, or improvements made within the spirit and principles of the present invention are intended to be within the scope of protection of the present invention.

Claims

1. A method for determining the explosion center location by combining dynamic pressure, characterized in that: include: Step S1: Determine the expected damage level A and the first initial explosion center position toward target A; Step S2: Calculate the distance between the first initial explosion center position and the target A according to the equivalent TNT equivalent of the current ammunition and the positive pressure action time of the ammunition at the target location; among them, are the target shock wave overpressure peak, specific impulse and dynamic pressure respectively; Step S3: Calculate the damage coefficient , the calculation formula is: ； in, are the damage thresholds corresponding to the target shock wave overpressure peak, specific impulse and dynamic pressure, respectively. The shock wave overpressure peak value, specific impulse and dynamic pressure damage coefficient related to the target and the positive pressure action time are respectively; then the actual damage level of the current damage coefficient calculation value is determined, Indicates no damage. Indicates minor injury; Indicates moderate damage. Severe injury; is a pre-calibrated constant; Step S4: Determine whether the substantial damage level is consistent with the A damage level. If so, determine the first initial ground center position as a reasonable position. If not, adjust the first initial ground center position according to the difference between the substantial damage level and the A damage level until the substantial damage level is consistent with the A damage level.

2. The method for determining the explosion center position by combining dynamic pressure according to claim 1, characterized in that: Also includes: In the process of iteratively adjusting the distance between the first initial ground center position and the A target based on step S4, determining the area range A that meets the A damage level; Then, the expected damage level B for target B is obtained, a position is randomly selected within the area range A as the second initial ground zero position, and then the area range B meeting the damage level B is determined in the same manner as the area range A meeting the damage level A. The intersection of the area range A and the area range B is taken to obtain a target area that simultaneously satisfies the damage level A expected for the A target and the damage level B expected for the B target.

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

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

5. The method for determining the explosion center position by combining dynamic pressure according to claim 4, characterized in that: The calculation formula is: ； ； ； ； ； in, 、 are the air pressures before and after the wave front, 、 are the air densities before and after the wave front, respectively; is the air adiabatic index, is the dynamic pressure velocity, is the air speed in front of the wave front.

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

7. A system for determining the center of explosion using combined dynamic pressure, comprising a memory, a processor, and a computer program stored in the memory and executable on the processor, characterized in that: When the processor executes the computer program, the method according to any one of claims 1 to 6 is implemented.

Images