Method for calculating the intensity of an explosive particle jet and related apparatus

Abstract

This application belongs to the field of explosion mechanics and provides a method and related equipment for calculating the intensity of an explosive particle jet. The method, under the condition that the explosive particles are driven by a high-pressure gas cloud, determines a first cloud region based on the mass distribution of the explosive particles, and determines a first jet region based on the packing density gradient of the explosive particles and the first cloud region; based on the mass of the explosive particles in the first jet region and the mass of the explosive particles in the first cloud region, it calculates the first mass specific gravity of the explosive particles in the first jet region; it determines the local packing density of the explosive particles at each jet location within the first jet region, and the minimum packing density of the explosive particles within the first jet region; based on the first mass specific gravity, local packing density, and minimum packing density, it obtains the first jet intensity. The above method can accurately determine the intensity of an explosive particle jet.

Description

Technical Field

[0001] This application relates to the field of explosion mechanics, specifically to a method for calculating the intensity of an explosive particle jet and related equipment. Background Technology

[0002] Explosive dispersion refers to the outward dispersion of particulate matter driven by an explosion. Numerous experiments have observed that during explosive dispersion, the particulate matter does not disperse uniformly outward, but rather in the form of a "finger-like jet," exhibiting strong localized aggregation. The formation of this jet structure directly affects the energy transfer and distribution of the explosive source, ultimately impacting the performance of the device. For example, in warheads (such as armor-piercing projectiles and thermobaric weapons), the intensity of the particulate jet directly affects the penetration depth, damage range, and target penetration capability.

[0003] Although the formation mechanism and influencing factors of jet structures have been studied in related technologies, there is a lack of reasonable methods for quantitative characterizing jet intensity. This limits the establishment of structure-property relationship between jet intensity and system structural parameters, and hinders the design optimization of charge structure. Summary of the Invention

[0004] This application provides a method and related equipment for calculating the intensity of explosive particle jets, in order to solve the technical problem of difficulty in accurately determining the intensity of jets due to the lack of quantitative characterization of jet intensity.

[0005] The first aspect of this application provides a method for calculating the intensity of an explosive particle jet. The method includes: determining a first cloud region based on the mass distribution of the explosive particles when the explosive particles are driven by a high-pressure gas cloud; determining a first jet region based on the packing density gradient of the explosive particles and the first cloud region; calculating a first mass proportion of the explosive particles in the first jet region based on the mass of the explosive particles in the first jet region and the mass of the explosive particles in the first cloud region; determining the local packing density of the explosive particles at each jet location in the first jet region and the minimum packing density of the explosive particles in the first jet region; and obtaining a first jet intensity based on the first mass proportion, the local packing density, and the minimum packing density.

[0006] In some embodiments, determining the first cloud region based on the mass distribution of the explosive particles includes: dividing the distribution range of the explosive particles into multiple concentric circle regions based on multiple virtual concentric circles, wherein the center of the multiple virtual concentric circles is the center of the high-pressure air mass; determining a second mass proportion of the explosive particles in each concentric circle region based on the ratio of the mass of the explosive particles in each concentric circle region to the total mass of all explosive particles; determining the inner boundary of the first cloud region based on the virtual concentric circles corresponding to the concentric circle regions where the second mass proportion is within a first preset range, and determining the outer boundary of the first cloud region based on the virtual concentric circles corresponding to the concentric circle regions where the second mass proportion is within a second preset range, wherein the second preset range is larger than the first preset range; and determining the first cloud region based on the area enclosed by the inner and outer boundaries of the first cloud region.

[0007] In some embodiments, determining the first jet region based on the packing density gradient of the explosive particles and the first cloud region includes: determining the inner boundary of the first jet region based on a virtual concentric circle where the packing density gradient of the explosive particles is at its minimum; determining the outer boundary of the first jet region based on the outer boundary of the first cloud region; and determining the first jet region based on the area enclosed by the inner and outer boundaries of the first jet region.

[0008] In some embodiments, calculating the first mass proportion of explosive particles in the first jet region based on the mass of explosive particles in the first jet region and the mass of explosive particles in the first cloud region includes: determining the first mass proportion based on the ratio of the mass of explosive particles in the first jet region to the mass of explosive particles in the first cloud region.

[0009] In some embodiments, the intensity of the first jet is expressed by the following formula: in, Indicates the intensity of the first jet. This represents the local packing density of the i-th explosive particle at the location of the jet within the first jet region. This represents the minimum packing density of explosive particles within the first jet region. This indicates the total number of explosive particles in the first jet region. This indicates the first mass ratio of the explosive particles within the first jet region. In some embodiments, the method further includes: acquiring a grayscale stained image of the dispersed explosive particles; determining a second cloud region and a second jet region based on the grayscale values ​​in the grayscale stained image; dividing the grayscale stained image into multiple grid regions; determining a first grayscale value for each grid region within the second jet region and a second grayscale value for each grid region within the second cloud region; determining a third grayscale value based on the minimum grayscale value within the second jet region; and calculating a second jet intensity based on the first grayscale value, the second grayscale value, and the third grayscale value.

[0010] In some embodiments, the second jet intensity is expressed by the following formula: in, Indicates the intensity of the second jet. This indicates the total number of grid regions within the second jet region. This represents the first grayscale value of the i-th grid region within the second jet region. This represents the third grayscale value within the second jet region. The radius representing the location of the explosive particle. This indicates the inner radius of the second jet region. This indicates the total number of grid regions within the second cloud / fog area. This represents the second grayscale value corresponding to the j-th grid region within the second cloud / fog area. This represents the inner radius of the second cloud / fog region. A second aspect of this application provides a device for calculating the intensity of an explosive particle jet. The device includes: a determining unit, configured to determine a first cloud region based on the mass distribution of the explosive particles and a first jet region based on the packing density gradient of the explosive particles and the first cloud region when the explosive particles are driven to move by a high-pressure gas cloud; a calculating unit, configured to calculate a first mass ratio of the explosive particles in the first jet region based on the mass of the explosive particles in the first jet region and the mass of the explosive particles in the first cloud region; the determining unit is further configured to determine the local packing density of the explosive particles at each jet location in the first jet region and the minimum packing density of the explosive particles in the first jet region; the calculating unit is further configured to obtain a first jet intensity based on the first mass ratio, the local packing density, and the minimum packing density.

[0011] A third aspect of this application provides an electronic device, which includes: a memory, a processor, and a computer program stored in the memory and executable on the processor. When the processor executes the computer program, it implements the steps in the method provided in the first aspect described above.

[0012] A fourth aspect of this application provides a computer-readable storage medium storing a computer program, which, when executed by a processor, implements the steps of the method provided in the first aspect.

[0013] In several embodiments of this application, by determining the first cloud zone and the first jet zone based on the distribution range of the explosive particles when the explosive particles are driven by a high-pressure gas cloud, a clear boundary is defined for evaluating the jet intensity of the explosive particles. Based on the mass of the explosive particles in the first jet zone and the mass of the explosive particles in the first cloud zone, the first mass proportion of the explosive particles in the first jet zone is calculated. The first mass proportion quantifies the proportion of explosive particles entering the first jet zone, revealing the density of explosive particles in the jet morphology, determining the local packing density of explosive particles at each jet location in the first jet zone, and the minimum packing density of explosive particles in the first jet zone. Then, based on the first mass proportion, local packing density, and minimum packing density, the first jet intensity is obtained. By integrating the distribution morphology and packing density of the explosive particles, the jet intensity value is directly quantified, which can accurately and efficiently quantitatively characterize the jet intensity of the explosive particles, providing a reliable basis for optimizing the design of explosive devices. Attached Figure Description

[0014] To more clearly illustrate the technical solutions of the embodiments of this application, the drawings used in the description of the embodiments of this application will be briefly introduced below. Obviously, the drawings described below are only some embodiments of this application. For those skilled in the art, other drawings can be obtained based on these drawings without creative effort.

[0015] Figure 1 This is a schematic diagram of the structure of the electronic device provided in the embodiments of this application.

[0016] Figure 2 This is a flowchart of the method for calculating the intensity of an explosive particle jet provided in the embodiments of this application.

[0017] Figure 3 This is a schematic diagram of a high-pressure air mass and a cloud / fog region provided in an embodiment of this application.

[0018] Figure 4 This is a schematic diagram of the first cloud and fog zone and the second jet zone provided in the embodiments of this application.

[0019] Figure 5 This is a flowchart of a method for calculating the intensity of an explosive particle jet according to another embodiment of this application.

[0020] Figure 6 This is a schematic diagram of a grayscale stained image of explosive particles provided in an embodiment of this application.

[0021] Figure 7 This is a result diagram provided in an embodiment of this application.

[0022] Figure 8 This is a structural diagram of the device for calculating the intensity of explosive particle jets provided in the embodiments of this application. Detailed Implementation

[0023] The technical solutions of the embodiments of this application will be clearly and completely described below with reference to the accompanying drawings. Obviously, the described embodiments are only some embodiments of this application, not all embodiments. Based on the embodiments of this application, all other embodiments obtained by those skilled in the art without creative effort are within the scope of protection of this application.

[0024] Hereinafter, the terms "first" and "second" are used for descriptive purposes only and should not be construed as indicating or implying relative importance or implicitly specifying the number of indicated technical features. Thus, a feature defined as "first" or "second" may explicitly or implicitly include one or more of that feature. In the description of embodiments of this application, words such as "exemplary" or "for example" are used to indicate examples, illustrations, or descriptions. Any embodiment or design described as "exemplary" or "for example" in the embodiments of this application should not be construed as being more preferred or advantageous than other embodiments or design solutions. Specifically, the use of words such as "exemplary" or "for example" is intended to present the relevant concepts in a specific manner.

[0025] Unless otherwise defined, all technical and scientific terms used herein have the same meaning as commonly understood by one of ordinary skill in the art to which this application pertains. The terminology used in this application's specification is for the purpose of describing particular embodiments only and is not intended to be limiting of the application. It should be understood that, unless otherwise stated, " / " in this application means "or". For example, A / B can mean A or B. "And / or" in this application is merely a description of the relationship between related objects, indicating that three relationships can exist. For example, A and / or B can represent: A alone, A and B simultaneously, and B alone. "At least one" refers to one or more. "More than one" refers to two or more. For example, at least one of a, b, or c can represent: a, b, c, a and b, a and c, b and c, and a, b, and c (seven cases).

[0026] In explosion dispersion experiments, the particulate medium was observed to exhibit a non-uniform dispersion characteristic of "finger-like jets," with strong local aggregation. Further research revealed that the formation mechanism of the jet is extremely complex: from a fluid dynamics perspective, the high-temperature and high-pressure gas generated by the explosion creates significant differences in the pressure gradient within the medium, while the interaction between particles drives the continuous evolution of the jet structure; regarding the physical properties of the particles themselves, differences in particle size and shape affect their trajectory in the airflow, thus altering the jet structure; in addition, different charge structures and properties also influence the release and propagation of explosion energy, ultimately affecting the jet.

[0027] Despite extensive experimental research into the formation mechanism of jet structures and the aforementioned influencing factors, current evaluations of jet intensity largely remain at the qualitative descriptive stage, lacking a reasonable and effective quantitative characterization method. For example, assessments are typically made through simple methods such as observing the jet's appearance and measuring its penetration depth. However, these qualitative methods lack accuracy and objectivity, failing to precisely reflect the actual magnitude and variation of jet intensity, and thus failing to meet the urgent needs of modern engineering fields for performance optimization and precise control of explosive devices.

[0028] To address the aforementioned issues, this application provides a method for calculating the intensity of an explosive particle jet. This method calculates the jet intensity using the mass and bulk density of the explosive particles, thus resolving the technical problem of difficulty in accurately determining the jet intensity due to the lack of quantitative characterization of the jet intensity.

[0029] Please see Figure 1 , Figure 1 This is a schematic diagram of the structure of an electronic device provided in an embodiment of this application. The electronic device 10 can be a mobile phone, tablet computer, smart wearable device, augmented reality (AR) / virtual reality (VR) device, laptop computer, netbook, or other electronic device with a display screen. This embodiment of the application does not impose any restrictions on the specific type of electronic device.

[0030] like Figure 1 As shown, the electronic device 10 may include a display screen 100, a communication module 101, a memory 102, a processor 103, an input / output (I / O) interface 104, and a bus 105. The processor 103 is coupled to the display screen 100, the communication module 101, the memory 102, and the I / O interface 104 via the bus 105.

[0031] It is understood that the structures illustrated in the embodiments of this application do not constitute a specific limitation on the electronic device 10. In other embodiments of this application, the electronic device 10 may include more or fewer components than illustrated, or combine some components, or split some components, or have different component arrangements. The illustrated components may be implemented in hardware, software, or a combination of software and hardware.

[0032] like Figure 2 The diagram shown is a flowchart of a method for calculating the intensity of an explosive particle jet according to an embodiment of this application. This method for calculating the intensity of an explosive particle jet can be applied to electronic devices (e.g.,...). Figure 1 The order of steps in the flowchart shown in the electronic device 10 can be changed, and some steps can be omitted.

[0033] S201, under the condition that the explosive particles are driven by a high-pressure gas mass, the first cloud zone is determined according to the mass distribution of the explosive particles, and the first jet zone is determined according to the packing density gradient of the explosive particles and the first cloud zone.

[0034] In some embodiments of this application, taking an underground explosive explosion scenario as an example, the explosive undergoes a violent chemical reaction in a very short time, generating a large amount of gas and releasing a large amount of heat. Within the explosion space, the gas volume increases sharply and the temperature rises significantly, causing the gas to be intensely compressed, forming a high-pressure gas cloud.

[0035] The initial pressure corresponds to the pressure difference between the interior of a high-pressure air mass and the external environment. Due to this pressure difference, the high-pressure air mass tends to expand, driving the surrounding particulate matter (i.e., explosive particles) composed of rocks, soil, etc., to move outward. When the energy of the high-pressure air mass is exhausted and the motion of the explosive particles becomes stable and no longer changes significantly over time, the flow field reaches equilibrium.

[0036] In the numerical simulation of the explosion system, a large number of explosive particles surround the outside of the high-pressure gas mass, for example, 1 million explosive particles. The number of explosive particles can be different in different explosion systems. Under the action of the initial pressure, the high-pressure gas mass is released as the explosion source. The released high-pressure gas mass drives a large number of explosive particles on the outside to move outward. When the flow field reaches equilibrium during the explosion, the first cloud and fog zone is determined according to the spatial dispersion range of the explosive particles. For example, a high-pressure air mass can be a spherical region, as shown in the reference. Figure 3 As shown, this is a quarter-segment region of a high-pressure air mass viewed from above. The area surrounding this quarter-segment region is the quarter-ring region of the first cloud / fog zone. The high-pressure air mass acts as the explosion source, driving a large number of explosive particles outwards. The first cloud / fog zone has its origin at the center of the high-pressure air mass, and its inner radius is [missing information]. The outer radius of the first cloud and fog zone is .

[0037] Electronic devices can acquire changes in the flow field pressure within the explosive system via pressure sensors connected to the explosive device. These pressure sensors can be miniature piezoresistive, piezoelectric, or similar types. When the difference between the flow field pressure and normal atmospheric pressure is within a preset range, and the flow field pressure no longer changes, it is determined that the flow field has reached equilibrium, at which point relevant data about the explosive particles are acquired.

[0038] For example, the preset range is set to (0, 2), and the normal atmospheric pressure is approximately 101.3 kPa. Sensor monitoring revealed that after the explosion began, the flow field pressure rose rapidly, reaching a peak of 500 kPa at 0.1 s. After a period of time, the pressure stabilized at 102 kPa at 1 s, and remained at 102 kPa thereafter. The difference between 102 kPa and 101.3 kPa was 0.7, which is within the preset range (0, 2). Since the flow field pressure no longer changed, it can be determined that the flow field reached equilibrium, indicating that the motion state of the explosive particles had stabilized and was no longer changing.

[0039] When the flow field reaches equilibrium, electronic devices can use numerical simulation software to acquire relevant data about the explosive particles. In this software, multiple monitoring points are pre-set in areas where the explosive particles are likely to appear. During the simulated explosion, the software automatically records relevant data about the explosive particles at each monitoring point, including particle velocity and mass. Numerical simulation software can include ANSYS Fluent, LS-DYNA, etc.

[0040] For example, when simulating the movement of explosive particles driven by the explosion of a high-pressure gas mass, 30 monitoring points are evenly arranged in a ring-shaped area 1m-5m away from the explosion center. During the explosion, the numerical simulation software automatically acquires relevant data such as the velocity, mass, packing density, packing density gradient, and position of the explosive particles at each monitoring point. In some embodiments of this application, the first cloud / fog zone can be determined based on relevant data, including: dividing the distribution range of the explosive particles into multiple concentric circle regions based on multiple virtual concentric circles, wherein the center of the multiple virtual concentric circles is the center of the high-pressure gas mass. The interval between each pair of adjacent virtual circles can be set and adjusted according to actual needs. In one example, the interval between each pair of adjacent virtual circles can be 1 mm; in another example, the interval between each pair of adjacent virtual circles can be 0.5 mm. Each concentric circle is sequentially numbered 0, 1, 2, 3…N. The radius of the circle is , It is a positive integer. (Reference) Figure 3 As shown, the virtual concentric circles are represented by the white dashed lines in the figure.

[0041] Calculate the ratio of the mass of the explosive particles within each concentric circle region to the total mass of all explosive particles to determine the second mass proportion of the explosive particles within each concentric circle region. Based on the virtual concentric circles corresponding to the concentric circle regions whose second mass proportions fall within a first preset range, determine the inner boundary of the first fog region, and define the radius of the virtual concentric circles within the first preset range as the inner radius of the first fog region. Based on the virtual concentric circles corresponding to the concentric circle regions whose second mass proportions fall within a second preset range, determine the outer boundary of the first fog region, and define the radius of the virtual concentric circles within the second preset range as the outer radius of the first fog region.

[0042] The second preset range is larger than the first preset range. For example, if the first preset range is 5%~6% and the second preset range is 95%~96%, the virtual concentric circle corresponding to the concentric circle region containing particles accounting for 5%~6% of the total mass is used as the inner boundary of the first cloud and fog region, and the virtual concentric circle corresponding to the concentric circle region containing particles accounting for 95%~96% of the total mass is used as the outer boundary of the first cloud and fog region.

[0043] In another embodiment, if there are at least two virtual concentric circles corresponding to concentric circle regions with a second mass proportion within a first preset range, the virtual concentric circle corresponding to the concentric circle region with the smallest second mass proportion within the first preset range is determined as the inner boundary of the first fog region, and the radius of the virtual concentric circle corresponding to the smallest second mass proportion within the first preset range is determined as the inner radius of the first fog region. If there are at least two virtual concentric circles corresponding to concentric circle regions with a second mass proportion within a second preset range, the virtual concentric circle corresponding to the concentric circle region with the largest second mass proportion within the second preset range is determined as the outer boundary of the first fog region, and the radius of the virtual concentric circle corresponding to the largest second mass proportion within the second preset range is determined as the inner radius of the first fog region. This embodiment ensures that the first fog region includes the vast majority of explosive particles.

[0044] In some embodiments of this application, the method for determining the first jet region includes: determining the packing density gradient corresponding to each concentric circle within the first cloud region (which can be...). Figure 4 The radial packing density gradient (in the first cloud region) represents the rate of change of the packing density of the explosive particles along the radial direction within the first cloud region. The packing density gradient can be directly obtained using numerical simulation software. A virtual concentric circle with the minimum packing density gradient is defined as the inner boundary of the first jet region, and the radius of this virtual concentric circle is defined as the inner radius of the first jet region. The outer boundary of the first cloud region is defined as the outer boundary of the first jet region. The area enclosed by the inner and outer boundaries of the first jet region is defined as the first jet region. In one example, such as Figure 4 As shown, Figure 4 This is a schematic diagram of the first cloud and fog zone and the first jet zone provided in the embodiments of this application. Figure 4 Figure (a) shows a local (if the first cloud region is a ring, only the 0-90° range is shown) Voronoi scatter plot of the packing density of the explosive particles within the first cloud region. Figure 4 As shown in Figure (a), in the figure The local packing density of the explosive particles can be represented. Figure (a) shows the variation of the local packing density with the radius of the virtual concentric circles within the first cloud region from 0° to 90°. As can be seen from the figure, the gray dashed line (YW10) represents the boundary within the first cloud region, and the gray dashed line corresponds to... The inner radius of the first cloud and fog region is represented by the blue dashed line (YW11), which indicates the inner boundary of the first jet region. The red dashed line (YW12) represents the inner radius of the first jet region's inner boundary, while the red dashed line represents the outer boundary of the first cloud region or the outer boundary of the first jet region. The outer radius represents the outer boundary of the first cloud zone or the outer boundary of the first jet zone, and different colors represent the packing density of the explosive particles. Figure 4 Figure (b) shows the radial packing density gradient of explosive particles in the first cloud area (only the 0~40° range is shown in the figure). This can represent the packing density gradient. Figure (b) shows how the packing density gradient varies with the radius of the virtual concentric circles within the first cloud region from 0° to 40°. Figure 4 As shown in Figure (b), based on the location corresponding to the minimum packing density gradient ( (Corresponding position), determine the inner boundary of the first jet region.

[0045] In this embodiment, the first cloud zone is determined based on the ratio of the mass of the explosive particles in the concentric circles to the mass of all explosive particles, a first preset range, and a second preset range. This ensures that the first cloud zone includes the vast majority of explosive particles. The first jet zone is determined based on the packing density gradient and the first cloud zone, which defines a clear boundary for assessing the jetting behavior of the explosive particles, thereby improving the accuracy of the assessment.

[0046] S202, Based on the mass of the explosive particles in the first jet region and the mass of the explosive particles in the first cloud region, calculate the first mass ratio of the explosive particles in the first jet region. In at least one embodiment of this application, the first mass proportion can represent the distribution of explosive particles dispersed into the first jet region. A larger first mass proportion indicates more explosive particles becoming jet particles, and a greater jet intensity. The electronic device can determine the first mass proportion as the ratio of the mass of explosive particles in the first jet region to the mass of explosive particles in the first cloud region. For example, if the mass of explosive particles in the first jet region is *a*, and the total mass of all explosive particles in the first cloud region is *b*, then the first mass proportion of explosive particles in the first jet region is *x* = *a* / *b*. The mass of explosive particles in the first jet region and the total mass of explosive particles can be directly obtained using numerical simulation software. For example, to obtain the mass of explosive particles in the first cloud region, the post-processing module of the numerical simulation software can be used. Filtering conditions can be set in the post-processing module, such as *r* being the radius of the location of the explosive particles, and the inner radius of the inner boundary of the first cloud region being... The outer radius of the first cloud and fog zone is , obtain <r< The explosive particles inside the ring were extracted using numerical simulation software. <r< The total mass of the explosive particles within the first cloud zone (within the circular ring) is calculated.

[0047] This embodiment calculates the first mass ratio of the explosive particles in the first jet region by comparing the mass of the explosive particles in the first jet region with the mass of the explosive particles in the first cloud region. The first mass ratio quantifies the proportion of explosive particles that become part of the jet, revealing the distribution of the jet and providing data reference for calculating the jet intensity.

[0048] S203, determine the local packing density of explosive particles at each location of the first jet within the first jet region and the minimum packing density of explosive particles within the first jet region.

[0049] In at least one embodiment of this application, when the first mass weight is the same, the higher the degree of aggregation of explosive particles within a certain azimuth angle range of the first jet region, the higher the jet intensity within that azimuth angle range. Explosive particles strongly aggregate at the azimuth angle of the jet, resulting in a relatively high local packing density. The first jet region is divided into multiple linear angles of equal angle according to a preset angle, such as 0.5°, 1°, etc. Based on the packing density corresponding to each linear angle, the local packing density and minimum packing density of explosive particles at each jet location within the first jet region are determined. The packing density value within each linear angle can also be directly obtained using numerical simulation software.

[0050] S204, based on the first mass specific gravity, local packing density and minimum packing density, the first jet intensity is obtained.

[0051] In at least one embodiment of this application, the first jet intensity represents the jet intensity obtained under numerical simulation, and the formula for determining the first jet intensity is: in, Indicates the first jet intensity. This represents the local packing density of the i-th explosive particle at the location of the jet within the first jet region. This represents the minimum packing density within the first jet region. This indicates the total number of explosive particles in the first jet region. This indicates the first mass ratio of the explosive particles within the first jet zone. This embodiment can quickly obtain a quantitative value characterizing the jet intensity by combining the distribution and packing density of the jet and based on the formula for determining the intensity of the first jet.

[0052] In several embodiments of this application, by determining the first cloud zone and the first jet zone based on the distribution range of the explosive particles when the explosive particles are driven by a high-pressure gas cloud, a clear boundary is defined for evaluating the jet intensity of the explosive particles. Based on the mass of the explosive particles in the first jet zone and the mass of the explosive particles in the first cloud zone, the first mass proportion of the explosive particles in the first jet zone is calculated. The first mass proportion quantifies the proportion of explosive particles entering the first jet zone, revealing the density of explosive particles in the jet morphology, determining the local packing density of explosive particles at each jet location in the first jet zone, and the minimum packing density of explosive particles in the first jet zone. Then, based on the first mass proportion, local packing density, and minimum packing density, the first jet intensity is obtained. By integrating the distribution morphology and packing density of the explosive particles, the jet intensity value is directly quantified, which can accurately and efficiently quantitatively characterize the jet intensity of the explosive particles, providing a reliable basis for optimizing the design of explosive devices.

[0053] like Figure 5 The diagram shown is a flowchart of a method for calculating the intensity of an explosive particle jet according to another embodiment of this application. This method can be applied to electronic devices; the order of the steps in the flowchart can be changed, and some steps can be omitted.

[0054] S501, acquire a grayscale stained image of the explosive particle dispersion, and determine the second cloud and fog region and the second jet region based on the grayscale values ​​in the grayscale stained image.

[0055] In at least one embodiment of this application, since it is difficult to directly obtain the packing density and first mass ratio in the experimental test explosion system, a grayscale stained image can be obtained by image processing of the high-speed photographic image of the explosion dispersion. The grayscale stained image can characterize the spatial distribution of the explosion particles.

[0056] Electronic equipment acquires high-speed photographs of explosions and processes them into grayscale images. Using the grayscale image of the device before the explosion and subsequent dispersion as a reference background image, the grayscale image is denoised by subtracting the pixel values ​​of each grayscale image from the reference background image. The denoised image is then divided into multiple concentric circular regions, with equal spacing between adjacent concentric circles. The grayscale values ​​within the concentric rings formed by adjacent concentric circles are summed and averaged to obtain the total grayscale radial distribution and the average grayscale distribution within the concentric rings, resulting in a grayscale distribution map. Based on pseudo-color mapping rules, the grayscale values ​​of the grayscale distribution map are mapped to a color space using piecewise linear functions or mathematical formulas (such as logarithmic mapping), resulting in a grayscale colored image. Taking the RGB color space as an example, in the RGB model, each color is composed of red (R), green (G), and blue (B) components combined in different proportions. Through pseudo-color mapping rules, the grayscale values ​​of the grayscale distribution map are mapped to different R, G, and B component values, thus determining a unique color representation in the color space. For example, low gray values ​​are mapped to a bluish tone, while high gray values ​​are mapped to a reddish tone.

[0057] Determine the first total grayscale value within each concentric circle region and the second total grayscale value of the entire grayscale-colored image. Based on the ratio of the first total grayscale value to the second total grayscale value, determine the grayscale ratio. The concentric circles corresponding to the concentric circle regions with grayscale ratios within a third preset range are defined as the inner boundary of the second fog region, and the radii of the concentric circles corresponding to the concentric circle regions within the third preset range are defined as the inner radius of the second fog region. The outer boundary of the second fog region is determined based on the concentric circles corresponding to the concentric circle regions with grayscale ratios within a fourth preset range, and the radii of the concentric circles corresponding to the concentric circle regions within the fourth preset range are defined as the outer radius of the second fog region. The second fog region is defined based on the area enclosed by the inner and outer boundaries of the second fog region.

[0058] The third preset range is smaller than the fourth preset range. For example, if the third preset range is 5%~6% and the fourth preset range is 95%~96%, the concentric circles corresponding to the concentric areas with a grayscale ratio of 5%~6% are used as the inner boundary of the second cloud area, and the concentric circles corresponding to the concentric areas with a grayscale ratio of 95%~96% are used as the outer boundary of the second cloud area.

[0059] In another embodiment, if there are at least two concentric circles corresponding to concentric circle regions with grayscale ratios within a third preset range, the concentric circle corresponding to the concentric circle region with the smallest grayscale ratio within the third preset range is determined as the inner boundary of the second fog region, and the radius of the concentric circle corresponding to the smallest grayscale ratio within the third preset range is determined as the inner radius of the second fog region. If there are at least two concentric circles corresponding to concentric circle regions with grayscale ratios within a fourth preset range, the concentric circle corresponding to the concentric circle region with the largest grayscale ratio within the fourth preset range is determined as the outer boundary of the second fog region, and the radius of the concentric circle corresponding to the largest grayscale ratio within the fourth preset range is determined as the inner radius of the second fog region. The grayscale gradient on each concentric circle boundary is determined. The grayscale gradient represents the rate of change of grayscale values ​​in space and can be obtained by calculating the difference in grayscale values ​​between pixels on each concentric circle boundary and pixels on adjacent concentric circle boundaries. The concentric circle boundary corresponding to the minimum grayscale gradient is defined as the inner boundary of the second jet region, and the outer boundary of the second cloud region is defined as the outer boundary of the second jet region. For example... Figure 6 As shown, Figure 6 This is a schematic diagram of a grayscale staining image of explosive particles provided in an embodiment of this application. Figure 6 The figure shows the range from 0° to 40°, illustrating how the grayscale value changes with the radius of the concentric circles within this range. This can be represented by grayscale values. As shown in the image, the larger the grayscale value, the stronger the aggregation of explosive particles. The white dashed line corresponds to... This represents the outer radius of the second cloud / fog region or the outer radius of the second jet region. The inner boundary of the second jet region is determined by the black dashed line corresponding to the minimum grayscale gradient. This represents the outer radius of the inner boundary of the second jet region.

[0060] This embodiment determines the second cloud and fog region by the grayscale ratio of the first total grayscale to the second total grayscale, the third preset range, and the fourth preset range, ensuring that the second cloud and fog region includes most of the grayscale data. Based on the minimum grayscale gradient, it ensures that the explosive particles aggregate most strongly at the boundary of the second jet region, thereby ensuring that the jet in the second jet region is the most abundant. This provides a clear boundary for assessing the jet situation of explosive particles, thereby improving the accuracy of the assessment.

[0061] S502, the grayscale stained image is divided into multiple grid regions, and the first grayscale value of each grid region in the second jet region and the second grayscale value of each grid region in the second cloud region are determined.

[0062] In this embodiment, the grayscale image is divided into multiple grid regions according to a preset grid side length. The preset grid side length can be set and adjusted according to actual conditions. For example, the preset grid can be a square grid, and the preset grid side length can be 0.1mm or 0.2mm. In another example, the preset grid can also be a grid of other shapes, such as a rectangle, etc., and this application does not impose specific limitations on this.

[0063] S503, determine the third gray value based on the minimum gray value in the second jet region.

[0064] S504, calculate the second jet intensity based on the first gray value, the second gray value, and the third gray value.

[0065] In at least one embodiment of this application, the second jet intensity represents the jet intensity obtained under experimental testing. The formula for determining the second jet intensity is: in, Indicates the intensity of the second jet. This indicates the total number of grid regions within the second jet region. This represents the first grayscale value of the i-th grid region within the second jet region. This represents the third grayscale value within the second jet region. The radius representing the location of the explosive particle. This indicates the inner radius of the second jet region. This indicates the total number of grid regions within the second cloud / fog area. This represents the second grayscale value corresponding to the j-th grid region within the second cloud / fog area. This represents the inner radius of the second cloud / fog region. To better understand the evaluation of the above-mentioned explosive particle jet intensity indicators, the following section combines... Figure 7 Further explanation. For example... Figure 7 Figures (a)-(f) shown below illustrate scatter plots of explosive particles from six explosion dispersion experiments. The diagram shows the jet intensity of the explosive particles. The area enclosed by the red line and the white dashed line is the second cloud region, and the area enclosed by the black line and the white dashed line is the second jet region. As can be seen from the diagram, from... Figure 7 In Figure (a), the finger-shaped intense jet in the second jet region gradually... Figure 7 In Figure (f), the filamentous jet in the second jet region gradually decreases in intensity. Figure 7 The (g)-(l) plots shown depict scatter plots of explosive particles from six numerical simulations of explosion dispersion. From... Figure 7 The finger-shaped intense jet in diagram (g) gradually reaches... Figure 7The fine filamentary jet in Figure (l) shows a gradual decrease in jet intensity, demonstrating the accuracy and reliability of the jet intensity index. This embodiment uses grayscale values ​​instead of bulk density to characterize the spatial distribution of explosive particles. Combined with the calculation formula for the second jet intensity, it efficiently calculates the value representing the jet intensity. The numerical simulation and experimental testing demonstrated the accuracy and reliability of the index. Figure 8 This is a structural diagram of the device for calculating the intensity of an explosive particle jet provided in this application embodiment. The device 11 for calculating the intensity of an explosive particle jet may include multiple functional modules composed of computer program segments. The computer programs of each program segment in the device 11 for calculating the intensity of an explosive particle jet may be stored in the memory of an electronic device and executed by at least one processor to perform (see details). Figure 2 (Description) Function for calculating the intensity of explosive particle jets.

[0066] In this embodiment, the explosive particle jet intensity calculation device 11 can be divided into multiple functional modules according to its function. The functional modules may include a determination module 110 and a calculation module 111. The determination module 110 is used to determine a first cloud zone based on the mass distribution of the explosive particles when the explosive particles are driven by a high-pressure gas cloud, and to determine a first jet zone based on the packing density gradient of the explosive particles and the first cloud zone. The calculation module 111 is used to calculate a first mass ratio of the explosive particles in the first jet zone based on the mass of the explosive particles in the first jet zone and the mass of the explosive particles in the first cloud zone. The determination module 110 is also used to determine the local packing density of the explosive particles at each jet location in the first jet zone and the minimum packing density of the explosive particles in the first jet zone. The calculation module 111 is also used to obtain a first jet intensity based on the first mass ratio, local packing density, and minimum packing density.

[0067] In several embodiments of this application, by determining the first cloud zone and the first jet zone based on the distribution range of the explosive particles when the explosive particles are driven by a high-pressure gas cloud, a clear boundary is defined for evaluating the jet intensity of the explosive particles. Based on the mass of the explosive particles in the first jet zone and the mass of the explosive particles in the first cloud zone, the first mass ratio of the explosive particles in the first jet zone is calculated. The first mass ratio quantifies the proportion of explosive particles entering the first jet zone, revealing the density of explosive particles in the jet morphology. The local packing density of explosive particles at each jet location in the first jet zone and the minimum packing density of explosive particles in the first jet zone are determined. Then, based on the first mass ratio, local packing density, and minimum packing density, the first jet intensity is obtained. By integrating the distribution morphology and packing density of the explosive particles, the jet intensity value is directly quantified, which can accurately and efficiently quantitatively characterize the jet intensity of the explosive particles, providing a reliable basis for optimizing the design of explosive devices.

[0068] The module referred to in this application is a series of computer program segments that can be executed by at least one processor and perform a fixed function, and which are stored in memory. In this embodiment, the limitation of the calculation device 11 for the intensity of the explosive particle jet can be found in the limitation of the calculation method for the intensity of the explosive particle jet above, and will not be repeated in detail here.

[0069] The display screen 100 can be a touch screen, specifically a touch-sensitive liquid crystal display device. Alternatively, the display screen 100 can be a non-touch screen. The display screen 100 can visualize information, for example, to display parameters such as the trajectory of explosive particles, the changing trend of flow field pressure, and the dynamic changing trend of explosive particle packing density gradient.

[0070] Communication module 101 may include a wired communication module and / or a wireless communication module. The wired communication module may provide one or more wired communication solutions such as Universal Serial Bus (USB) and Controller Area Network (CAN). The wireless communication module may provide one or more wireless communication solutions such as Wireless Fidelity (Wi-Fi), Bluetooth (BT), mobile communication networks, Frequency Modulation (FM), Near Field Communication (NFC), and Infrared (IR). The memory 102 may include one or more random access memory (RAM) and one or more non-volatile memory (NVM). The RAM can be directly read and written by the processor 103, and can be used to store executable programs (such as machine instructions) of the operating system or other running programs, as well as user and application data.

[0071] Random access memory can include static random-access memory (SRAM), dynamic random-access memory (DRAM), synchronous dynamic random-access memory (SDRAM), double data rate synchronous dynamic random-access memory (DDR SDRAM), etc.

[0072] Non-volatile memory can also store executable programs and user and application data, and can be pre-loaded into random access memory for direct reading and writing by the processor 103. Non-volatile memory can include disk storage devices and flash memory.

[0073] Memory 102 is used to store one or more computer programs. The one or more computer programs are configured to be executed by processor 103. The one or more computer programs include multiple instructions that, when executed by processor 103, enable a method for calculating the intensity of an explosive particle jet to be executed on electronic device 10.

[0074] In other embodiments, the electronic device 10 also includes an external memory interface for connecting to an external memory to expand the storage capacity of the electronic device 10.

[0075] Processor 103 may include one or more processing units, such as an application processor (AP), a modem processor, a graphics processing unit (GPU), an image signal processor (ISP), a controller, a video codec, a digital signal processor (DSP), a baseband processor, and / or a neural network processing unit (NPU). These different processing units may be independent devices or integrated into one or more processors.

[0076] The processor 103 provides computing and control capabilities. For example, the processor 103 is used to execute computer programs stored in the memory 102 to implement the above-described method for calculating the intensity of the explosive particle jet.

[0077] I / O interface 104 is used to provide a channel for user input or output. For example, I / O interface 104 can be used to connect various input and output devices, such as mouse, keyboard, touch device, display screen, etc., so that users can enter information or visualize information.

[0078] Bus 105 is used at least to provide a channel for communication between communication modules 101, memory 102, processor 103, and I / O interface 104 in electronic device 10.

[0079] This application also provides a computer-readable storage medium storing a computer program, which includes program instructions. When the program instructions are executed, the method implemented can refer to the methods in the above embodiments of this application.

[0080] The computer-readable storage medium can be the internal memory of the electronic device described in the above embodiments, such as the hard disk or memory of the electronic device. Alternatively, the computer-readable storage medium can be an external storage device of the electronic device, such as a plug-in hard disk, smart media card (SMC), secure digital card (SD), flash card, etc., provided on the electronic device.

[0081] In some embodiments, a computer-readable storage medium may include a stored program area and a stored data area, wherein the stored program area may store an operating system, an application program required for at least one function, etc.; and the stored data area may store data created based on the use of the electronic device, etc.

[0082] In the above embodiments, the descriptions of each embodiment have different focuses. For parts that are not described in detail or recorded in a certain embodiment, please refer to the relevant descriptions of other embodiments.

[0083] Those skilled in the art will recognize that the units and algorithm steps of the various examples described in conjunction with the embodiments disclosed herein can be implemented in electronic hardware, or a combination of computer software and electronic hardware. Whether these functions are implemented in hardware or software depends on the specific application and design constraints of the technical solution. Those skilled in the art can use different methods to implement the described functions for each specific application, but such implementation should not be considered beyond the scope of this application.

[0084] In the embodiments provided in this application, it should be understood that the disclosed apparatus / terminal devices and methods can be implemented in other ways. For example, the apparatus / terminal device embodiments described above are merely illustrative. For instance, the division of modules or units is only a logical functional division, and in actual implementation, there may be other division methods. For example, multiple units or components may be combined or integrated into another system, or some features may be ignored or not executed. Furthermore, the coupling or direct coupling or communication connection shown or discussed may be through some interfaces; the indirect coupling or communication connection between devices or units may be electrical, mechanical, or other forms.

[0085] The units described as separate components may or may not be physically separate. The components shown as units may or may not be physical units; that is, they may be located in one place or distributed across multiple network units. Some or all of the units can be selected to achieve the purpose of this embodiment according to actual needs.

[0086] The above embodiments are only used to illustrate the technical solutions of this application, and are not intended to limit them. Although this application has been described in detail with reference to the foregoing embodiments, those skilled in the art should understand that modifications can still be made to the technical solutions described in the foregoing embodiments, or equivalent substitutions can be made to some of the technical features. Such modifications or substitutions do not cause the essence of the corresponding technical solutions to deviate from the spirit and scope of the technical solutions of the embodiments of this application, and should all be included within the protection scope of this application.

Claims

1. A method for calculating the intensity of an explosive particle jet, characterized in that, The method includes: In the case where the explosive particles are driven by a high-pressure gas mass, a first cloud region is determined based on the mass distribution of the explosive particles, and a first jet region is determined based on the packing density gradient of the explosive particles and the first cloud region. Based on the mass of the explosive particles in the first jet region and the mass of the explosive particles in the first cloud region, the first mass ratio of the explosive particles in the first jet region is calculated. Determine the local packing density of explosive particles at each location of the first jet within the first jet region and the minimum packing density of explosive particles within the first jet region. Based on the first mass specific gravity, the local bulk density, and the minimum bulk density, the first jet intensity is obtained, wherein the first jet intensity is calculated using the following formula: in, Indicates the intensity of the first jet. This represents the local packing density of the i-th explosive particle at the location of the jet within the first jet region. This represents the minimum packing density of explosive particles within the first jet region. This indicates the total number of explosive particles in the first jet region. This indicates the first mass ratio of the explosive particles within the first jet region.

2. The method for calculating the intensity of an explosive particle jet according to claim 1, characterized in that, Determining the first cloud region based on the mass distribution of the explosive particles includes: The distribution range of the explosive particles is divided into multiple concentric circle regions based on multiple virtual concentric circles, wherein the center of the multiple virtual concentric circles is the center of the high-pressure gas mass. The second mass ratio of the explosive particles in each concentric circle region is determined based on the ratio of the mass of the explosive particles in each concentric circle region to the total mass of all explosive particles. The inner boundary of the first cloud and fog area is determined based on the virtual concentric circles corresponding to the concentric circle areas where the second mass specific gravity is within the first preset range, and the outer boundary of the first cloud and fog area is determined based on the virtual concentric circles corresponding to the concentric circle areas where the second mass specific gravity is within the second preset range, wherein the second preset range is greater than the first preset range. The first cloud and fog region is determined based on the area enclosed by the inner and outer boundaries of the first cloud and fog region.

3. The method for calculating the intensity of an explosive particle jet according to claim 2, characterized in that, The step of determining the first jet region based on the packing density gradient of the explosive particles and the first cloud region includes: The inner boundary of the first jet region is determined based on the virtual concentric circles where the packing density gradient of the explosive particles is minimized. The outer boundary of the first jet region is determined based on the outer boundary of the first cloud region; The first jet region is determined based on the area enclosed by the inner and outer boundaries of the first jet region.

4. The method for calculating the intensity of an explosive particle jet according to claim 1, characterized in that, The calculation of the first mass ratio of the explosive particles in the first jet region based on the mass of the explosive particles in the first jet region and the mass of the explosive particles in the first cloud region includes: The first mass ratio is determined based on the ratio of the mass of the explosive particles in the first jet zone to the mass of the explosive particles in the first cloud zone.

5. The method for calculating the intensity of an explosive particle jet according to claim 1, characterized in that, The method further includes: Obtain a grayscale stained image of the dispersion of the explosive particles, and determine the second cloud and fog region and the second jet region based on the grayscale values ​​in the grayscale stained image; The grayscale stained image is divided into multiple grid regions, and the first grayscale value of each grid region in the second jet region and the second grayscale value of each grid region in the second cloud region are determined. The third gray value is determined based on the minimum gray value within the second jet region; The second jet intensity is calculated based on the first gray value, the second gray value, and the third gray value.

6. The method for calculating the intensity of an explosive particle jet according to claim 5, characterized in that, The second jet intensity is determined by the following formula: in, Indicates the intensity of the second jet. This indicates the total number of grid regions within the second jet region. This represents the first grayscale value of the i-th grid region within the second jet region. This represents the third grayscale value within the second jet region. The radius representing the location of the explosive particle. This indicates the inner radius of the second jet region. This indicates the total number of grid regions within the second cloud / fog area. This represents the second grayscale value corresponding to the j-th grid region within the second cloud / fog area. This represents the inner radius of the second cloud / fog region.

7. A device for calculating the intensity of an explosive particle jet, characterized in that, The device includes: The determining unit is used to determine a first cloud region based on the mass distribution of the explosive particles when the explosive particles are driven to move by a high-pressure gas mass, and to determine a first jet region based on the packing density gradient of the explosive particles and the first cloud region. The calculation unit is used to calculate the first mass ratio of the explosive particles in the first jet zone based on the mass of the explosive particles in the first jet zone and the mass of the explosive particles in the first cloud zone. The determining unit is also used to determine the local packing density of the explosive particles at the location of each jet in the first jet region, and the minimum packing density of the explosive particles in the first jet region. The calculation unit is further configured to obtain a first jet intensity based on the first mass specific gravity, the local bulk density, and the minimum bulk density, wherein the first jet intensity is calculated using the following formula: in, Indicates the intensity of the first jet. This represents the local packing density of the i-th explosive particle at the location of the jet within the first jet region. This represents the minimum packing density of explosive particles within the first jet region. This indicates the total number of explosive particles in the first jet region. This indicates the first mass ratio of the explosive particles within the first jet region.

8. An electronic device, characterized in that, include: A memory, a processor, and a computer program stored in the memory and executable on the processor, wherein the processor, when executing the computer program, implements a method for calculating the intensity of an explosive particle jet as described in any one of claims 1 to 6.

9. A computer-readable storage medium, characterized in that, The computer-readable storage medium stores a computer program that, when executed by a processor, implements a method for calculating the intensity of an explosive particle jet as described in any one of claims 1 to 6.

Images