Deduction method and device for explosive loading DDT position and storage medium

Abstract

The invention discloses an explosive loading DDT position deduction method. The method comprises the following steps: S1, determining take-off time and minimum take-off time according to acquired data; s2, charging explosive, shell parameters and PDV probe position information are input, a detonation wave velocity range is preset, and a resolution ratio PP I value is set; s3, the position information is converted into a Cartesian coordinate system, a medicine loading physical model is drawn, random points are generated according to a PP I set value, and space distribution is completed; s4, the maximum distance difference and the minimum distance difference are determined according to the PDV take-off time difference, an outer ball and an inner ball are drawn, and a point set range is obtained; s5, generating random points, screening in a point set range, and searching and selecting points meeting conditions; and S6, outputting the coordinates of the selected point location, converting the coordinates into longitude and latitude forms, and remarking the elevation of the coordinates to determine the area position where the DDT occurs. The method can deduce and determine the position of a region where DDT possibly occurs in the accident response evolution process of charge.

Description

Technical Field

[0001] The invention relates to the technical field of explosive charge safety, and in particular to a deduction method, equipment and storage medium for a DDT position of an explosive charge. Background Art

[0002] The deflagration-detonation process of explosive charges (also known as the DDT process) refers to the process in which the deflagration propagating in the explosive gradually develops and transforms into a detonation. The DDT process involves complex physical and chemical changes and is widely present in the detonation of explosives and the action of weapons and ammunition. It is also one of the safety factors that must be considered in the production and storage of explosives, and is one of the important performance parameters for the safety design and evaluation of explosives throughout their life cycle.

[0003] The deflagration-detonation process (DDT process) of explosive charges undergoes a very complex multi-stage chemical and physical reaction, and its combustion rate, pressure and temperature in the combustion reaction zone all undergo changes of orders of magnitude, with the time varying from tens of microseconds to several milliseconds.

[0004] In order to explore the reaction evolution process of the charge under accident stimulation and clarify the conditions for the transition of the charge deflagration-detonation reaction intensity, it is necessary to carry out artificial ignition engineering tests of the charge to observe and analyze its accident reaction evolution mechanism. The determination of the location of the area where the deflagration-detonation process of the explosive charge occurs is particularly important for the analysis of the charge accident reaction evolution process.

[0005] At present, the existing technology mainly uses PDV (shell expansion velocity-time information) to capture the charge shell expansion velocity-time development and analyze the evolution process of charge reaction intensity. Generally, when DDT occurs, the velocity-time curve will show a microsecond jump feature (jumping a kilometer within 0.5-1 microsecond). In the research process, researchers generally determine the DDT position based on the velocity jump characteristics collected by PDV measuring points at different locations, the assumption of consistent detonation wave velocity and the velocity jump time difference, combined with engineering experience for comprehensive analysis and deduction.

[0006] However, this traditional method is time-consuming, fails to quantify the DDT area, and the reliability of the deduced DDT location cannot be verified. In addition, no alternative technical solutions have been found in the prior art.

[0007] To this end, the present application proposes a method for deducing the position of explosive charge DDT to solve the above technical problems. Summary of the invention

[0008] The main purpose of the present invention is to provide a method for deducing the DDT position of explosive charge, and to deduce the DDT area of ​​explosive based on the charge shell expansion velocity-time curve to solve the technical problems raised in the background technology.

[0009] The present invention adopts the following technical solutions to solve the above technical problems:

[0010] A method for deducing the position of explosive charge DDT includes the following specific operation steps:

[0011] S1. Collect speed-time curve data at different PDV measuring points to determine the take-off time T i and minimum take-off time T min ;

[0012] S2 inputs the charge explosive, basic shell size parameters, and PDV probe location longitude and latitude information, and presets the detonation wave velocity range V based on the detonation wave velocity. min ~V max , set the resolution PPI value;

[0013] S3. Convert the longitude and latitude positions into a Cartesian coordinate system, draw a charge physical model based on the actual input charge explosive and shell basic size parameters, complete coloring and transparency settings to distinguish between the explosive and the shell, and then generate random points based on the PPI setting value to complete the space allocation;

[0014] S4. Determine the maximum distance difference D based on the PDV take-off time difference at different positions max and the minimum distance difference D min As the outer sphere range and inner sphere range respectively, draw the outer sphere and inner sphere respectively:

[0015] D max =V max *T i -V min *T min

[0016] D min =V min *T i -V min *T min

[0017] Get the point set range D min ~D max ;

[0018] S5. Generate random points based on the point set range as the determination condition, and then find and select points that meet the preset screening conditions within the point set range according to the point set range;

[0019] S6. Output the coordinates of the selected points, convert the obtained Cartesian coordinates into the form of longitude and latitude, and note the elevation to determine the regional location where DDT occurs.

[0020] Preferably, the specific operation process of step S1 includes:

[0021] S11. Input the lowest recorded detonation velocity, and record the time when the lowest detonation velocity is reached in the velocity-time curve data collected at each measuring point as the take-off time T i ;

[0022] S12. Iteratively compare the take-off time T of each measuring point i , select the earliest time as the minimum take-off time T min .

[0023] Preferably, the basic size parameters of the charge explosive and shell in step S2 include: the explosive radius of the spherical explosive or the hemispherical explosive and the shell surface thickness.

[0024] Preferably, the specific operation steps of searching for points within the point set that meet the screening conditions in step S5 include:

[0025] S51. Preset the appropriate point as (x, y, z), and the DDT position is O (x0, y0, z0), and the point with the minimum take-off time is Min (x min ,y min , z min ), preset Function condition function:

[0026]

[0027] Use the Function conditional function to filter points;

[0028] S52. Preset the actual medicine ball screening formula to perform point screening again. The medicine ball screening formula is:

[0029] a 2 ≤x 2 +y 2 +z 2 ≤ b 2

[0030] Where a represents the inner shell diameter of the explosive ball, and b represents the inner surface of the explosive shell.

[0031] On the other hand, the present invention further discloses a computer-readable storage medium storing a computer program, wherein when the computer program is executed by a processor, the processor executes the steps of the above method.

[0032] On the other hand, the present invention further discloses a computer device, including a memory and a processor, wherein the memory stores a computer program, and when the computer program is executed by the processor, the processor executes the steps of the above method.

[0033] It can be seen from the above technical solution that the present invention provides a method for deducing the position of explosive charge DDT. Compared with the prior art, the present invention has the following advantages:

[0034] The method of the present invention deduces the DDT occurrence area based on the shell expansion velocity-time information (PDV), and can deduce and determine the regional location where the charge may undergo deflagration-detonation (DDT) during the evolution of the accident reaction, thereby helping relevant researchers to quickly determine the location where the charge may undergo deflagration-detonation (DDT) based on PDV data, and providing assistance and reference for the analysis of the charge reaction evolution mechanism and engineering test analysis.

[0035] It should be understood that the content described in this section is not intended to identify the key or important features of the embodiments of the present invention, nor is it intended to limit the scope of the present invention. Other features of the present invention will become easy to understand through the following description. Of course, it is not necessary to achieve all of the advantages described above simultaneously for any product implementing the present invention. BRIEF DESCRIPTION OF THE DRAWINGS

[0036] The drawings constituting a part of the present application are used to provide a further understanding of the present invention. The exemplary 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:

[0037] Figure 1 It is a schematic diagram of the overall operation flow of the method of the present invention;

[0038] Figure 2 A schematic diagram of generating a random point allocation space for the present invention;

[0039] Figure 3 It is a schematic diagram of the speed-time curve obtained by PDV at different positions of the present invention;

[0040] Figure 4 A schematic diagram of setting up a 2D hemispherical physical model of the present invention;

[0041] Figure 5 It is a linear schematic diagram of simulating PDV data for the 2D hemispherical physical model of the present invention;

[0042] Figure 6 It is a schematic diagram of the output result of the DDT program of the 2D hemispherical physical model of the present invention;

[0043] Figure 7 A schematic diagram of setting up a 3D spherical physical model of the present invention;

[0044] Figure 8 It is a linear schematic diagram of simulation data of the 3D spherical physical model of the present invention;

[0045] Fig. 9It is a schematic diagram for comparing the output results of the DDT program of the present invention. DETAILED DESCRIPTION

[0046] The technical solutions in the embodiments of the present invention will be clearly and completely described below in conjunction with the drawings in the embodiments of the present invention. Obviously, the described embodiments are only part of the embodiments of the present invention, rather than all of the embodiments. In the absence of conflict, the embodiments in this application and the features in the embodiments can be combined with each other. Based on the embodiments in the present invention, all other embodiments obtained by ordinary technicians in this field without making creative work are within the scope of protection of the present invention.

[0047] In the embodiment, see Figures 1 to 9 .

[0048] like Figure 1 As shown, a method for deducing the DDT position of explosive charge proposed in an embodiment of the present invention deduces the DDT occurrence area based on the shell expansion velocity-time information (PDV), and can find the explosion point position verified by 2D and 3D simulation, including the following specific operation steps:

[0049] S1. Collect speed-time curve data at different PDV measuring points and determine the take-off time T i and minimum take-off time T min , which specifically include:

[0050] S11. Input the lowest recorded detonation velocity, and record the time when the lowest detonation velocity is reached in the velocity-time curve data collected at each measuring point as the take-off time T i At this time, the minimum detonation velocity as a screening threshold can accurately capture the starting point of the detonation wave propagation, thereby facilitating the accurate determination of the initial moment of the explosion reaction and providing a reliable time reference for subsequent analysis;

[0051] S12. Iteratively compare the take-off time T of each measuring point i , select the earliest time as the minimum take-off time T min ;

[0052] S2 inputs the charge explosive, basic shell size parameters, and PDV probe location longitude and latitude information, and presets the detonation wave velocity range V based on the detonation wave velocity. min ~V max , set the resolution PPI value, where the basic size parameters of the charge explosive and shell include: the explosive radius and shell surface thickness of spherical explosive or hemispherical explosive;

[0053] S3. Convert the longitude and latitude positions into Cartesian coordinate system, draw the charge physical model according to the actual input charge explosive and shell basic size parameters, complete the coloring and transparency settings to distinguish the explosive and shell, generate random points according to the PPI setting value, and complete the space allocation. Figure 2 ,shown, shows the point set for drawing the generated random points within the medicine ball;

[0054] The charge physics model constructed at this time has high simulation accuracy, which can enhance the model's understanding and prediction capabilities of the DDT occurrence mechanism. At the same time, the drawn charge physics model can also play a role in clearly displaying the structural characteristics, thereby achieving the effect of optimizing the visualization of spatial distribution, and finally making it easier for researchers to intuitively understand the experimental layout;

[0055] S4. By analyzing the actual physical process of the detonation wave acting on the shell surface after the charge has DDT, it is believed that DDT occurs before the minimum take-off time, such as Figure 3 As shown, assuming that the DDT occurs at T D , and the location where DDT occurs is assumed to be O. Due to the limitations of experimental diagnostic technology, the experiment cannot obtain the actual T D and O, but an equation without two parameters can be established according to the PDV take-off time difference at different positions and based on this equation, the maximum distance difference and the minimum distance difference can be determined. Therefore, according to the PDV take-off time difference at different positions, the maximum distance difference D is determined. max and the minimum distance difference D min As the outer sphere range and inner sphere range respectively, write a program in Matlab to draw the outer sphere and inner sphere with the maximum distance difference and the minimum distance difference respectively:

[0056]

[0057] It is assumed that T i All units of D are in us, so max The unit is m. Here T i represents all the sampled take-off times except the minimum take-off time, T min is the minimum time obtained;

[0058] In summary, we can get the point set range D min ~D max ;

[0059] S5. Using the point set range as a determination condition to determine the generation of random points, and then finding and selecting points within the point set range that meet the preset screening conditions based on the point set range, wherein the specific operation steps of finding points within the point set range that meet the screening conditions include:

[0060] S51. Preset the appropriate point as (x, y, z), and the DDT position is O (x0, y0, z0), and the point with the minimum take-off time is Min (x min ,y min , z min ), preset Function condition function:

[0061]

[0062] Use the Function conditional function to filter points;

[0063] S52. Preset the actual medicine ball screening formula to perform point screening again. The medicine ball screening formula is:

[0064] a 2 ≤x 2 +y 2 +z 2 ≤ b 2

[0065] Where a represents the inner shell diameter of the explosive ball, and b represents the inner surface of the explosive shell;

[0066] S6. Output the coordinates of the selected points, convert the obtained Cartesian coordinates into the form of longitude and latitude, and note the elevation to determine the regional location where DDT occurs.

[0067] In summary, the method of the present invention is based on the assumption of uniform detonation wave velocity and can deduce and determine the regional location where deflagration-detonation (DDT) may occur in the charge during the evolution of the accident reaction according to the take-off time difference of PDV data, thereby helping relevant researchers to quickly determine the location where deflagration-detonation (DDT) may occur in the charge based on PDV data, and further provide help and reference for the analysis of charge reaction evolution mechanism and engineering test analysis.

[0068] Test example:

[0069] In a specific embodiment, the present invention uses a practical example to apply the deduction method of the explosive charge DDT position provided in the present application, thereby accurately deducing and calculating the explosive deflagration-detonation (DDT) area based on the charge shell expansion velocity-time curve. The method specifically includes:

[0070] The 2D hemispherical surface model is drawn using Autodyn numerical simulation software, such as Figure 4 As shown, the green area is the explosive PBX9404 with a radius of 180 mm, and the blue area is the surface of the steel shell with a thickness of 35 mm;

[0071] In order to ensure that the detonation wave propagates in a straight line inside the sphere and meet the DDT program design assumptions, no inner spherical shell is set inside the explosive;

[0072] at this time Figure 4 The light blue dots are the PDV speed measurement setting points. Since the simulation is a 2D sphere, the specific locations are distinguished by longitude and latitude. The longitude is considered to be 0°, and the hemisphere is distinguished by north and south latitude. The red is the setting burst point position, which is (90 north latitude, 0 longitude, R150mm), converted into xyz coordinates (0,0,150);

[0073] The detonation wave velocity is set to 8.8 km / s. According to the numerical simulation settings, the simulation is started to obtain the shell expansion distance-time curves of the above 9 measurement points as shown in the figure. Figure 5 As shown;

[0074] At this time, based on the data results, the take-off time of each point is obtained: (90 South latitude, longitude 0, 42.2us), (67.5 South latitude, longitude 0, 41.5us), (45 South latitude, longitude 0, 39.1us), (22.5 South latitude, longitude 0, 35.6us), (latitude 0, longitude 0, 30.48us), (22.5 North latitude, longitude 0, 24.6us), (45 North latitude, longitude 0, 18.1us), (67.5 North latitude, longitude 0, 12.2).

[0075] Input the above data into the process of this method, set PPI to 100, and get the output result as follows Figure 6 As shown:

[0076] Z=135.7672, 135.8832, 133.9374, 138.4290, 135.1927, 135.3652, 136.3969, 135.8453, 137.6664;

[0077] At this time, the prediction accuracy is 91.8%, and the deviation from the actual explosion point location is 8.2%.

[0078] Furthermore, in order to improve the accuracy of the present application in actual use, the method may further continue to perform the following operation steps during specific use:

[0079] The 3D sphere model was drawn using Autodyn numerical simulation software for 3D numerical simulation experimental verification, such as Figure 7 As shown, the sphere radius of the explosive PBX9404 is still 180 mm, and the shell surface thickness is 35 mm, which is consistent with the 2D spherical drawing size.

[0080] In order to ensure that the detonation wave propagates in a straight line inside the sphere and meet the DDT program design assumptions, an inner spherical shell is not set inside the explosive. Figure 7The light blue points are PDV velocity measurement setting points, totaling 20 measurement points. The measurement point locations are still set according to the engineering longitude and latitude. The xyz coordinates of the explosion point are (147.8, 61.2, 0), 160 mm away from the center of the ball, and the detonation wave velocity is still set to 8.8 km / s.

[0081] According to the numerical simulation settings, the simulation is started to obtain the shell expansion distance-time curves of the above 20 measuring points as shown in Figure 8 As shown;

[0082] According to the data results, the take-off time of each point is obtained: (45 North, longitude 0, 9.5us), (60 North, longitude 0, 13.9us), (45 South, longitude 0, 22.61us), (60 South, longitude 0, 26.2us), (45 North, 45 East, 15.1us), (60 North, 45 East, 16.5us), (45 South, 45 East, 25.3us), (60 South, 45 East, 28.3us), (45 North, 135 East, 31.4us), (60 North, 135 East, 28us), (45 South, 135 East, 37.7us), (60 South 0, East longitude 135, 36.6us), (N45, W45, 14.21us), (N60, W45, 16.1us), (S45, W45, 25.1us), (S60, W45, 28.1us), (N45, W135, 31.5us), (N60, W135, 28us), (S45, W135, 38.1us), (S60, W135, 36.7us), input the above data into the process of this method, set PPI to 100, and input the detonation wave velocity range from 6 to 10 km / s. The output result is as follows Fig. 9 As shown;

[0083] List the main coordinate point results:

[0084] X=153.6292, 132.4214, 134.2175, 136.6336;

[0085] Y=41.3746;49.2331, 47.1006, 50.1317;

[0086] Z=8.6386, -4.2508, -0.1095, 10.6757;

[0087] Distance from the center of the ball H = 159.3375, 141.3414142.2421, 145.9312;

[0088] Compared with the simulated burst point position (147.8, 61.2, 0, R159.97mm), the burst point position closest calculated by MATLAB is (153.63, 41.37, 8.64, H159.34mm), with a prediction deviation of 0.4%.

[0089] The above results show that the calculation method provided by the present application can well find the explosion point location verified by 2D and 3D simulations, and the calculation program constructed using MATLAB based on the calculation method of the present application can, based on the assumption of uniform detonation wave velocity and the take-off time difference of PDV data, deduce and determine the location of the area where deflagration-detonation may occur in the evolution of the accident reaction of the charge, which can help relevant researchers to quickly determine the location where DDT may occur in the charge based on PDV data, and provide help and reference for the analysis of charge reaction evolution mechanism and engineering test analysis.

[0090] On the other hand, the present invention further discloses a computer-readable storage medium storing a computer program, wherein when the computer program is executed by a processor, the processor executes the steps of the above method.

[0091] On the other hand, the present invention further discloses a computer device, including a memory and a processor, wherein the memory stores a computer program, and when the computer program is executed by the processor, the processor executes the steps of the above method.

[0092] In another embodiment provided by the present application, a computer program product including instructions is also provided, which, when executed on a computer, enables the computer to execute the method for deducing the DDT position of any of the above-mentioned embodiments.

[0093] It is understandable that the system provided by the embodiment of the present invention corresponds to the method provided by the embodiment of the present invention, and the explanation, examples and beneficial effects of the relevant contents can refer to the corresponding parts in the above method.

[0094] The embodiment of the present application also provides an electronic device, including a processor, a communication interface, a memory and a communication bus, wherein the processor, the communication interface and the memory communicate with each other through the communication bus.

[0095] Memory, used to store computer programs;

[0096] The processor is used to implement the above-mentioned method of deducing the DDT position of the explosive charge when executing the program stored in the memory.

[0097] The communication bus mentioned in the above electronic device can be a Peripheral Component Interconnect (PCI) bus or an Extended Industry Standard Architecture (EISA) bus, etc. The communication bus can be divided into an address bus, a data bus, a control bus, etc.

[0098] The communication interface is used for communication between the above electronic device and other devices.

[0099] The memory may include a random access memory (RAM) or a non-volatile memory (NVM), such as at least one disk memory. Optionally, the memory may also be at least one storage device located away from the aforementioned processor.

[0100] The above-mentioned processor can be a general-purpose processor, including a central processing unit (CPU), a network processor (NP), etc.; it can also be a digital signal processor (DSP), an application-specific integrated circuit (ASIC), a field-programmable gate array (FPGA) or other programmable logic devices, discrete gates or transistor logic devices, and discrete hardware components.

[0101] It should also be noted that electronic devices also include terminal devices, which can also be called terminals, user equipment (UE), mobile stations (MS), mobile terminals (MT), etc. Terminal devices can be mobile phones, smart TVs, wearable devices, tablet computers (Pad), computers with wireless transceiver functions, virtual reality (VR) terminal devices, augmented reality (AR) terminal devices, wireless terminals in industrial control, wireless terminals in self-driving, wireless terminals in remote medical surgery, wireless terminals in smart grids, wireless terminals in transportation safety, wireless terminals in smart cities, wireless terminals in smart homes, etc. The embodiments of the present application do not limit the specific technology and specific device form adopted by the terminal devices.

[0102] In the above embodiments, it can be implemented in whole or in part by software, hardware, firmware or any combination thereof. When implemented using software, it can be implemented in whole or in part in the form of a computer program product. The computer program product includes one or more computer instructions. When the computer program instructions are loaded and executed on a computer, the process or function described in the embodiment of the present application is generated in whole or in part. The computer may be a general-purpose computer, a special-purpose computer, a computer network, or other programmable device. The computer instructions may be stored in a computer-readable storage medium, or transmitted from one computer-readable storage medium to another computer-readable storage medium, for example, the computer instructions may be transmitted from a website site, computer, server or data center by wired (e.g., coaxial cable, optical fiber, digital subscriber line (DSL)) or wireless (e.g., infrared, wireless, microwave, etc.) mode to another website site, computer, server or data center. The computer-readable storage medium may be any available medium that a computer can access or a data storage device such as a server or data center that includes one or more available media integrated. The available medium may be a magnetic medium, (e.g., a floppy disk, a hard disk, a tape), an optical medium (e.g., a DVD), or a semiconductor medium (e.g., a solid-state drive Solid State Disk (SSD)), etc.

[0103] The above description is only a preferred embodiment of the present invention and is not intended to limit the present invention. Any modifications, equivalent substitutions, improvements, etc. made within the spirit and principle of the present invention should be included in the protection scope of the present invention.

[0104] In addition, it should be noted that if the embodiments of the present invention involve directional indications (such as up, down, left, right, front, back, etc.), the directional indications are only used to explain the relative position relationship, movement status, etc. between the components in a certain specific posture. If the specific posture changes, the directional indication will also change accordingly.

[0105] In addition, if there are descriptions involving "first", "second", etc. in the embodiments of the present invention, the descriptions of "first", "second", etc. are only used for descriptive purposes and cannot be understood as indicating or implying their relative importance or implicitly indicating the number of technical features indicated. Therefore, the features defined as "first" and "second" may explicitly or implicitly include at least one of the features. In addition, the meaning of "and / or" appearing in the full text includes three parallel schemes. Taking "A and / or B" as an example, it includes scheme A, or scheme B, or schemes that A and B meet at the same time. In addition, in the embodiments of the present invention, "multiple" refers to more than two. In addition, the technical solutions between the various embodiments can be combined with each other, but it must be based on the ability of ordinary technicians in the field to implement. When the combination of technical solutions is contradictory or cannot be implemented, it should be considered that such a combination of technical solutions does not exist and is not within the scope of protection required by the present invention.

Claims

1. A method for deducing the position of explosive charge DDT, characterized in that: The specific steps are as follows: S1. Collect speed-time curve data at different PDV measuring points to determine the take-off time T i and minimum take-off time T min ; S2 inputs the charge explosive, basic shell size parameters, and PDV probe location longitude and latitude information, and presets the detonation wave velocity range V based on the detonation wave velocity. min ~V max , set the resolution PPI value; S3. Convert the longitude and latitude positions into a Cartesian coordinate system, draw a charge physical model based on the actual input charge explosive and shell basic size parameters, complete coloring and transparency settings to distinguish between the explosive and the shell, and then generate random points based on the PPI setting value to complete the space allocation; S4. Determine the maximum distance difference D based on the PDV take-off time difference at different positions max and the minimum distance difference D min As the outer sphere range and inner sphere range respectively, draw the outer sphere and inner sphere respectively: Get the point set range D min ~D max ; S5. Generate random points based on the point set range as the determination condition, and then find and select points that meet the preset screening conditions within the point set range according to the point set range; S6. Output the coordinates of the selected points, convert the obtained Cartesian coordinates into the form of longitude and latitude, and note the elevation to determine the regional location where DDT occurs.

2. The method for deducing the position of explosive charge DDT as claimed in claim 1, characterized in that: The specific operation process of step S1 includes: S11. Input the lowest recorded detonation velocity, and record the time when the lowest detonation velocity is reached in the velocity-time curve data collected at each measuring point as the take-off time T i ; S12. Iteratively compare the take-off time T of each measuring point i , select the earliest time as the minimum take-off time T min .

3. The method for deducing the position of explosive charge DDT as claimed in claim 1, characterized in that: The basic size parameters of the charge explosive and shell in step S2 include: the explosive radius of the spherical explosive or the hemispherical explosive and the shell surface thickness.

4. The method for deducing the position of explosive charge DDT as claimed in claim 1, characterized in that: The specific operation steps of searching for points within the point set that meet the screening conditions in step S5 include: S51. Preset the appropriate point as (x, y, z), and the DDT position is O(x0, y0, z0), and the point with the minimum take-off time Min(x min ,y min , z min ), preset Function condition function: Use the Function conditional function to filter points; S52. Preset the actual medicine ball screening formula to perform point screening again. The medicine ball screening formula is: to 2 ≤x 2 +y 2 +z 2 ≤b 2 Where a represents the inner shell diameter of the explosive ball, and b represents the inner surface of the explosive shell.

5. A computer-readable storage medium, characterized in that: A computer program is stored, and when the computer program is executed by a processor, the processor is caused to perform the steps of the method according to any one of claims 1 to 4.

6. A computer device, characterized in that: The method comprises a memory and a processor, wherein the memory stores a computer program, and when the computer program is executed by the processor, the processor executes the steps of the method according to any one of claims 1 to 4.

Images