Method and processor for calculating shock wave pressure for internal blast in a closed structure

Abstract

The application discloses a kind of for closed structure inside explosion Shockwave pressure calculation method and processor, the Shockwave pressure calculation method includes: with each wall of closed structure as mirror surface, mirror image is carried out to actual explosion source and according to the quantity of actual explosion source determines mirror image explosion source set;Determine the coordinate of observation point of closed structure inside explosion;According to the coordinate of observation point and mirror image explosion source set determines virtual influence explosion source set;Based on actual explosion source and virtual influence explosion source set, calculate the Shockwave pressure of observation point when closed structure inside explosion occurs.The Shockwave pressure calculation method for closed structure inside explosion has simple method, easy to realize, is advantageous to promote the Shockwave pressure calculation efficiency and accuracy of closed structure inside explosion.

Description

TECHNICAL FIELD

[0001] The present application relates to the field of explosion technology, in particular to a shock wave pressure calculation method and processor for internal explosion of a closed structure. BACKGROUND

[0002] The explosion of explosives in a structure is different from the explosion in an open space. The explosion shock wave is reflected on different walls to enhance and interact with each other. The time history curve of the shock wave pressure presents the characteristics of multiple peaks, and the rules are complex. The internal walls of the structure will also be subjected to multiple pressure impacts. There are mainly two ways to calculate the internal explosion pressure of the structure in the prior art: one is numerical simulation, and the other is empirical model estimation. Among them, the numerical simulation method has large calculation amount and long calculation time, and the empirical estimation model is difficult to handle the complex reflection, diffraction and other situations of the explosion shock wave, which leads to difficulty in ensuring the accuracy of the calculation results. SUMMARY

[0003] The purpose of the present application is to provide a shock wave pressure calculation method and processor for internal explosion of a closed structure, which has the advantages of simple method, easy implementation, and is beneficial to improve the calculation efficiency and accuracy of the shock wave pressure of the internal explosion of the closed structure.

[0004] In order to achieve the above-mentioned purpose, one aspect of the present application provides a shock wave pressure calculation method for internal explosion of a closed structure, which comprises:

[0005] mirroring the actual explosive source with each wall surface of the closed structure as a mirror surface and determining a mirror explosive source set according to the number of actual explosive sources;

[0006] determining the coordinates of the observation point of the internal explosion of the closed structure;

[0007] determining a virtual influence explosive source set according to the coordinates of the observation point and the mirror explosive source set;

[0008] calculating the shock wave pressure of the observation point during the internal explosion of the closed structure based on the actual explosive source and the virtual influence explosive source set.

[0009] In the embodiments of the present application, mirroring the actual explosive source with each wall surface of the closed structure as a mirror surface and determining a mirror explosive source set according to the number of actual explosive sources comprises:

[0010] determining the number of mirror images of the actual explosive source;

[0011] determining the mirror explosive source set according to the number of mirror images and the number of actual explosive sources.

[0012] In the embodiments of the present application, the number of mirror images is 2.

[0013] In embodiments of the present invention, determining the virtual impact source set based on the coordinates of the observation point and the mirrored explosion source set includes:

[0014] The influence domain segmentation plane corresponding to the secondary mirror explosion source is determined based on the coordinates of the observation point and the coordinates of the secondary mirror explosion source in the mirror explosion source set.

[0015] The coordinates of the intersection point are determined based on the influence domain segmentation plane and the mirror plane corresponding to the secondary mirror explosion source.

[0016] The set of virtual impact sources is determined based on the coordinates of the intersection point and the coordinates of the observation point.

[0017] In embodiments of the present invention, determining the set of virtual impact explosion sources based on the coordinates of the intersection point and the coordinates of the observation point includes:

[0018] The first and second coordinate values ​​are determined based on the coordinates of the observation point and the coordinates of the secondary mirror explosion source.

[0019] The first and second preset ranges are determined based on the set of dimensional parameters of the sealed structure.

[0020] Determine whether the first coordinate value belongs to the first preset range;

[0021] Determine whether the second coordinate value belongs to the second preset range;

[0022] If the first coordinate value belongs to the first preset range and the second coordinate value belongs to the second preset range, the secondary mirror explosion source is determined to be a virtual influence explosion source.

[0023] The set of virtual impact sources is determined based on the virtual impact sources.

[0024] In embodiments of the present invention, determining the first coordinate value and the second coordinate value of the intersection point includes:

[0025] The first and second coordinate values ​​are determined based on the coordinates of the observation point and the coordinates of the secondary mirror explosion source.

[0026] In embodiments of the present invention, determining the first preset range and the second preset range based on the set of size parameters of the sealed structure includes:

[0027] The first preset range is determined based on the first parameter in the set of size parameters;

[0028] The second preset range is determined based on the second parameter in the set of size parameters, wherein the first parameter is a parameter belonging to the same coordinate direction as the first coordinate value, and the second parameter is a parameter belonging to the same coordinate direction as the second coordinate value.

[0029] In an embodiment of the present invention, calculating the shock wave pressure at the observation point during an explosion inside a confined structure based on the actual explosion source and a set of virtual impact explosion sources includes:

[0030] Calculate the first shock wave pressure of the actual explosion source;

[0031] Calculate the second shock wave pressure of the virtual impact explosion source set;

[0032] The shock wave pressure at the observation point is determined by the sum of the first shock wave pressure and the second shock wave pressure.

[0033] A second aspect of the present invention provides a processor configured to execute the above-described method for calculating shock wave pressure in the event of an explosion inside a sealed structure.

[0034] A third aspect of the present invention provides a machine-readable storage medium storing instructions for causing a machine to execute the above-described method for calculating the shock wave pressure of an explosion inside a sealed structure.

[0035] The above technical solution involves mirroring the actual explosion source to obtain a set of mirrored explosion sources. Then, based on the set of mirrored explosion sources, a set of virtually influenced explosion sources is further obtained. Finally, the shock wave pressure at the observation point during the explosion inside the sealed structure is calculated based on the actual explosion source and the set of virtually influenced explosion sources. This method is simple, easy to implement, and helps to improve the efficiency and accuracy of calculating the shock wave pressure during an explosion inside a sealed structure. Attached Figure Description

[0036] The accompanying drawings are provided to further illustrate embodiments of the present invention and form part of the specification. They are used together with the following detailed description to explain the embodiments of the present invention, but do not constitute a limitation thereof. In the drawings:

[0037] Figure 1 This is a flowchart illustrating the shock wave pressure calculation method in an embodiment of the present invention;

[0038] Figure 2 This is a schematic diagram of the mirrored explosion source in an embodiment of the present invention;

[0039] Figure 3 This is a schematic diagram of influence domain segmentation in an embodiment of the present invention. Specific Implementation

[0040] The specific embodiments of the present invention will be described in detail below with reference to the accompanying drawings. It should be understood that the specific embodiments described herein are for illustration and explanation only and are not intended to limit the scope of the present invention.

[0041] An embodiment of the present invention provides a method for calculating the shock wave pressure of an explosion inside a sealed structure, such as...Figure 1 As shown, the shock wave pressure calculation method includes the following steps:

[0042] Step S101: Mirror the actual explosion source and determine the set of mirrored explosion sources formed by the actual explosion source based on the number of actual explosion sources.

[0043] In one embodiment of the present invention, step S101: determining the set of mirrored explosion sources formed by the actual explosion source based on the actual explosion source further includes steps S201-S202, wherein:

[0044] Step S201: Determine the number of times the actual explosion source is mirrored.

[0045] When an explosive shock wave collides with the rigid wall of a sealed structure, it reflects a shock wave. Therefore, the pressure at a certain point in the sealed structure is the superposition of the shock wave pressures from the actual explosion source and the mirrored explosion source obtained by mirroring. In this embodiment, the number of actual explosion sources B0 is selected as one. By mirroring the actual explosion source B0 with each wall of the pipe gallery structure as a mirror, multiple primary mirrored explosion sources can be obtained. By mirroring the primary mirrored explosion source again with each wall of the pipe gallery structure as a mirror, multiple secondary mirrored explosion sources can be obtained. The more times the mirroring is performed, the more mirrored explosion sources are obtained, and the more accurate the subsequent calculation results will be. However, the amount of calculation will increase exponentially, which will lead to a significant decrease in calculation efficiency. Therefore, in order to ensure both the accuracy of the calculation results and high calculation efficiency, the number of mirroring is selected as two in the internal explosion pressure calculation method of this example.

[0046] Step S202: Determine the set of mirrored sources based on the number of mirroring attempts and the actual number of sources.

[0047] In this embodiment, before executing step S202, a set of dimensional parameters for the sealed structure can be obtained first. Further, the sealed structure in this embodiment can be a pipe gallery structure with a rectangular cross-section. The internal explosion pressure calculation method is implemented by a processor in a computer, and the set of dimensional parameters includes the length a of the pipe gallery structure. Figure 2 (middle y direction), width b ( Figure 2 (middle x direction) and high c ( Figure 2 Before implementation, operators can measure the length (a), width (b), and height (c) of the utility tunnel structure using a length measuring tool (such as a tape measure), and then input these three parameters into the computer processor.

[0048] In this embodiment, the actual explosion source B0 is mirrored once, and the mirrored explosion source B is obtained after one mirroring. iThere are 6, namely B1, B2, B3, B4, B5, and B6. The coordinates of the actual explosion source B0 are (Ex, Ey, Ez) (in this coordinate system, the lower left corner of the rear end of the pipe gallery structure is the origin); further, according to the principle of explosion source mirroring, the coordinates of B1 are (Ex, Ey, -Ez), the coordinates of B2 are (Ex, -Ey, Ez), the coordinates of B3 are (2*a-Ex, Ey, Ez), the coordinates of B4 are (Ex, 2*b-Ey, Ez), the coordinates of B5 are (-Ex, Ey, Ez), and the coordinates of B6 are (Ex, Ey, 2*c-Ez); a single mirrored explosion source B i The secondary image obtained after mirroring reveals source B ij (Where i = 1, 2, 3, ..., 6; j = 1, 2, 3, ..., 6, and i ≠ j) There are 30 such sources, among which the secondary mirror source B is obtained by mirroring the z = 0 plane. 12 B 13 B 14 B 15 and B 16 The coordinates of all points are (Ex, Ey, -Ez); the secondary mirror image source B is obtained by mirroring the plane with respect to the y=0 plane. 21 B 23 B 24 B 25 and B 26 The coordinates are all (Ex, -Ey, Ez); the secondary mirror image source B is obtained by mirroring the plane x = a. 31 B 32 B 34 B 35 and B 36 The coordinates are all (2*a-Ex, Ey, Ez); the secondary mirror image source B is obtained by mirroring the plane y=b. 41 B 42 B 43 B 45 and B 46 The coordinates are all (Ex, 2*b-Ey, Ez); the secondary mirror image source B is obtained by mirroring the plane at x=0. 51 B 52 B 53 B 54 and B 56 The coordinates are all (-Ex, Ey, Ez); the secondary mirror image source B is obtained by mirroring the plane z=c. 61 B 62 B 63 B 64 and B 65 The coordinates are all (Ex, Ey, 2*c-Ez). The mirrored explosion source set B in this embodiment includes all the aforementioned primary mirrored explosion sources B. iand secondary mirror source B ij .

[0049] Step S102: Determine the coordinates of the observation point for the explosion inside the sealed structure.

[0050] Specifically, in this embodiment, there is one observation point with coordinates (x, y, z).

[0051] Step S103: Determine the set of virtual impact explosion sources based on the coordinates of the observation point and the set of mirror explosion sources.

[0052] In one embodiment of the present invention, step S103: determining the virtual impact explosion source set based on the observation point and the mirror explosion source set further includes steps S301-S303, wherein:

[0053] Step S301: Determine the influence domain segmentation plane corresponding to the secondary mirrored explosion source based on the coordinates of the observation point and the coordinates of the secondary mirrored explosion source in the mirrored explosion source set.

[0054] In this embodiment, the coordinates of the observation point are (x, y, z). Secondary mirrored explosion sources with different coordinate positions in the mirrored explosion source set have different influence domains. Therefore, secondary mirrored explosion sources with different coordinate positions correspond to different influence domain dividing planes. The influence domain dividing plane is determined based on the influence domain dividing line. Specifically, the coordinates of the observation point and the coordinates of a certain secondary mirrored explosion source in the mirrored explosion source set are obtained first. Then, the straight line where the two coordinate points are located is determined based on the coordinates of the observation point and the coordinates of the certain secondary mirrored explosion source. This straight line is the influence domain dividing line. The expression of the influence domain dividing line can be solved based on the coordinates of the observation point and the coordinates of the secondary mirrored explosion source. The plane that passes through the influence domain dividing line and is perpendicular to the vertical section of the pipe gallery structure is the influence domain dividing plane. The influence domain dividing plane divides the pipe gallery structure into two parts in space. The space part that is close to the primary mirrored explosion source corresponding to the secondary mirrored explosion source is the influence domain corresponding to the secondary mirrored explosion source. The influence domain segmentation plane corresponding to each secondary mirror explosion source can be determined using the above method. Secondary mirror explosion sources with the same coordinate position have the same influence domain segmentation plane, but two secondary mirror explosion sources with the same coordinate position may correspond to different influence domains.

[0055] Step S302: Determine the coordinates of the intersection point based on the influence domain segmentation plane and the mirror plane corresponding to the secondary mirror explosion source.

[0056] Specifically, such as Figure 3 As shown, there is an intersection point S between the influence domain segmentation line corresponding to the secondary mirror explosion source and the mirror plane corresponding to the secondary mirror explosion source. In this embodiment, the secondary mirror explosion source B... 12 B 13 B 14 B15 and B 16 The coordinates of the corresponding intersection point S1 are (x1, y1, z1); the secondary mirror image source B 21 B 23 B 24 B 25 and B 26 The coordinates of the corresponding intersection point S2 are (x2, y2, z2); the secondary mirror image source B 31 B 32 B 34 B 35 and B 36 The coordinates of the corresponding intersection point S3 are (x3, y3, z3); ​​the secondary mirror image source B 41 B 42 B 43 B 45 and B 46 The coordinates of the corresponding intersection point S4 are (x4, y4, z4); the secondary mirror image source B 51 B 52 B 53 B 54 and B 56 The coordinates of the corresponding intersection point S5 are (x5, y5, z5); the secondary mirror image source B 61 B 62 B 63 B 64 and B 65 The coordinates of the corresponding intersection point S6 are (x6, y6, z6).

[0057] Step S303: Determine the set of virtual impact explosion sources based on the coordinates of the intersection point and the coordinates of the observation point.

[0058] In one embodiment of the present invention, step S303: determining the set of virtual impact explosion sources based on the coordinates of the intersection point and the coordinates of the observation point further includes steps S401-S406, wherein:

[0059] Step S401: Determine the first and second coordinate values ​​of the intersection point.

[0060] In one embodiment of the present invention, step S401: determining the first coordinate value and the second coordinate value of the intersection point includes:

[0061] The first and second coordinate values ​​are determined based on the coordinates of the observation point and the coordinates of the secondary mirror explosion source.

[0062] Specifically, the secondary mirror image reveals source B. 12 B 13 B 14 B 15 and B 16The coordinates of all points are (Ex, Ey, -Ez), the mirror plane is z = 0, the coordinates of the intersection point S1 are (x1, y1, z1), where x1 = x + (-z)*(Ex - x) / (Ez - z), y1 = y + (-z)*(Ey - y) / (Ez - z), z1 = 0, the first coordinate value of the intersection point S1 is x1, and the second coordinate value of the intersection point S1 is y1;

[0063] Secondary mirror image source B 21 B 23 B 24 B 25 and B 26 The coordinates of all points are (Ex, -Ey, Ez), the mirror plane is y = 0, and the coordinates of the intersection point S2 are (x2, y2, z2), where x2 = x + (-y)*(Ex - x) / (Ey - y), y2 = 0, z2 = z + (-y)*(Ez - z) / (Ey - y), the first coordinate value of the intersection point S2 is x2, and the second coordinate value of the intersection point S2 is z2.

[0064] Secondary mirror image source B 31 B 32 B 34 B 35 and B 36 The coordinates of all points are (2*a-Ex, Ey, Ez), the mirror plane is x=a, and the coordinates of the intersection point S3 are (x3, y3, z3), where x3=a, y3=y+(ax)*(Ey-y) / (Ex-x), z3=z+(ax)*(Ez-z) / (Ex-x), the first coordinate value of the intersection point S3 is y3, and the second coordinate value of the intersection point S3 is z3;

[0065] Secondary mirror image source B 41 B 42 B 43 B 45 and B 46 The coordinates of all points are (Ex, 2*b-Ey, Ez), the mirror plane is y=b, and the coordinates of the intersection point S4 are (x4, y4, z4), where x4=x+(by)*(Ex-x) / (Ey-y), y4=b, z4=z+(by)*(Ez-z) / (Ey-y), the first coordinate value of the intersection point S4 is x4, and the second coordinate value of the intersection point S4 is z4.

[0066] Secondary mirror image source B 51 B 52 B 53 B 54 and B 56The coordinates of all points are (-Ex, Ey, Ez), the mirror plane is x5 = 0, and the coordinates of the intersection point S5 are (x5, y5, z5), where x5 = 0, y5 = y + (-x)*(Ey - y) / (Ex - x), z5 = z + (-x)*(Ez - z) / (Ex - x), the first coordinate value of the intersection point S5 is y5, and the second coordinate value of the intersection point S5 is z5.

[0067] Secondary mirror image source B 61 B 62 B 63 B 64 and B 65 The coordinates of all points are (Ex, Ey, 2*c-Ez), the mirror plane is z6=c, the coordinates of the intersection point S6 are (x5, y5, z5), where x6=x+(cz)*(Ex-x) / (Ez-z), y6=y+(cz)*(Ey-y) / (Ez-z), z6=c, the first coordinate value of the intersection point S6 is x6, and the second coordinate value of the intersection point S6 is y6.

[0068] Step S402: Determine the first preset range and the second preset range based on the set of size parameters.

[0069] In one embodiment of the present invention, step S402: determining the first preset range and the second preset range based on the set of size parameters of the sealed structure includes steps S501-S502, wherein:

[0070] Step S501: Determine the first preset range based on the first parameter in the size parameter set;

[0071] Step S502: Determine the second preset range based on the second parameter in the size parameter set, wherein the first parameter is a parameter belonging to the same coordinate direction as the first coordinate value, and the second parameter is a parameter belonging to the same coordinate direction as the second coordinate value.

[0072] Specifically, in this embodiment, the first coordinate value corresponding to the intersection point S1 is x1, the second coordinate value corresponding to the intersection point S1 is y1, the first preset range corresponding to the intersection point S1 is 0-a, and the second preset range corresponding to the intersection point S1 is 0-b.

[0073] If the first coordinate value of the intersection point S1 is x1 and the second coordinate value of the intersection point S1 is y1, then the first preset range of the intersection point S1 is 0-a and the second preset range of the intersection point S1 is 0-b.

[0074] The first coordinate value corresponding to the intersection point S2 is x2, and the second coordinate value corresponding to the intersection point S2 is z2. Then the first preset range corresponding to the intersection point S2 is 0-a, and the second preset range corresponding to the intersection point S2 is 0-c.

[0075] If the first coordinate value corresponding to the intersection point S3 is y3 and the second coordinate value corresponding to the intersection point S3 is z3, then the first preset range corresponding to the intersection point S3 is 0 - b, and the second preset range corresponding to the intersection point S3 is 0 - c;

[0076] If the first coordinate value corresponding to the intersection point S4 is x4 and the second coordinate value corresponding to the intersection point S4 is z4, then the first preset range corresponding to the intersection point S4 is 0 - a, and the second preset range corresponding to the intersection point S4 is 0 - c;

[0077] If the first coordinate value corresponding to the intersection point S5 is y5 and the second coordinate value corresponding to the intersection point S5 is z5, then the first preset range corresponding to the intersection point S5 is 0 - b, and the second preset range corresponding to the intersection point S5 is 0 - c;

[0078] If the first coordinate value corresponding to the intersection point S6 is x6 and the second coordinate value corresponding to the intersection point S6 is y6, then the first preset range corresponding to the intersection point S6 is 0 - a, and the second preset range corresponding to the intersection point S6 is 0 - b.

[0079] Step S403: Determine whether the first coordinate value belongs to the first preset range;

[0080] Step S404: Determine whether the second coordinate value belongs to the second preset range;

[0081] Step S405: When the first coordinate value belongs to the first preset range and the second coordinate value belongs to the second preset range, determine that the secondary mirror explosion source is a virtual influence explosion source.

[0082] Specifically, when the first coordinate value x1 corresponding to the intersection point S1 satisfies the condition 0 < x1 < a, and the second coordinate value y1 corresponding to the intersection point S1 satisfies the condition 0 < y1 < b, then determine that the secondary mirror explosion sources B 12 、B 13 、B 14 、B 15 and B 16 are virtual influence explosion sources; otherwise, the secondary mirror explosion sources B 12 、B 13 、B 14 、B 15 and B 16 are not virtual influence explosion sources;

[0083] When the first coordinate value x2 corresponding to the intersection point S2 satisfies the condition 0 < x2 < a, and the second coordinate value z2 corresponding to the intersection point S2 satisfies the condition 0 < z2 < c, then determine that the secondary mirror explosion sources B 21 、B 23 、B 24 、B 25 and B 26 are virtual influence explosion sources; otherwise, the secondary mirror explosion sources B21 and B 23 and B 24 and B 25 and B 26 is not a virtual impact explosion source;

[0084] When the first coordinate value y3 corresponding to the intersection point S3 satisfies the condition 0 < y3 < b, and the second coordinate value z3 corresponding to the intersection point S3 satisfies the condition 0 < z3 < c, then the secondary mirror explosion source B is determined 31 and B 32 and B 34 and B 35 and B 36 is a virtual impact explosion source; otherwise, the secondary mirror explosion source B 31 and B 32 and B 34 and B 35 and B 36 is not a virtual impact explosion source;

[0085] When the first coordinate value x4 corresponding to the intersection point S4 satisfies the condition 0 < x4 < a, and the second coordinate value z4 corresponding to the intersection point S4 satisfies the condition 0 < z4 < c, then the secondary mirror explosion source B is determined 41 and B 42 and B 43 and B 45 and B 46 is a virtual impact explosion source; otherwise, the secondary mirror explosion source B 41 and B 42 and B 43 and B 45 and B 46 is not a virtual impact explosion source;

[0086] When the first coordinate value y5 corresponding to the intersection point S5 satisfies the condition 0 < y5 < b, and the second coordinate value z5 corresponding to the intersection point S5 satisfies the condition 0 < z5 < c, then the secondary mirror explosion source B is determined 51 and B 52 and B 53 and B 54 and B 56 is a virtual impact explosion source; otherwise, the secondary mirror explosion source B 51 and B 52 and B 53 and B 54 and B 56 is not a virtual impact explosion source;

[0087] When the first coordinate value x6 corresponding to the intersection point S6 satisfies the condition 0 < x6 < a, and the second coordinate value y6 corresponding to the intersection point S6 satisfies the condition 0 < y6 < b, then the secondary mirror explosion source B is determined 61 and B 62 and B63 B 64 and B 65 The virtual source is the source of the explosion; otherwise, the secondary mirror image is the source B. 61 B 62 B 63 B 64 and B 65 It is not a virtual source of the outbreak.

[0088] Step S406: Determine the set of virtual impact sources based on the virtual impact sources.

[0089] Specifically, the virtual impact source set includes all primary mirrored sources and secondary mirrored sources that meet the requirements (i.e., the first coordinate value of the intersection associated with the secondary mirrored source belongs to the first preset range and the second coordinate value belongs to the second preset range).

[0090] Step S104: Calculate the shock wave pressure at the observation point when an explosion occurs inside a closed structure, based on the actual explosion source and the set of virtual impact explosion sources.

[0091] In one embodiment of the present invention, step S104: calculating the shock wave pressure at the observation point based on the actual explosion source and the set of virtual impact explosion sources further includes steps S601-S603, wherein:

[0092] Step S601: Calculate the first shock wave pressure of the actual explosion source;

[0093] Specifically, the pressure of the first shock wave is calculated according to the following formula:

[0094]

[0095] Wherein, P1 is the first shock wave pressure; P0 is the ambient pressure at the actual explosion source location; ΔP1 is the peak overpressure at the actual explosion source location, in MPa; and ρ1 is the shock wave density at the actual explosion source location. This represents the shock wave velocity at the actual location of the explosion source. This represents the increase in shock wave density at the actual location of the explosion source.

[0096] Calculate ΔP1 using the following formula:

[0097]

[0098] Where Z1 is the first proportional distance of the actual explosion source and can be calculated according to the following formula:

[0099]

[0100] Where R1 is the distance between the observation point and the actual blast source, in meters; W1 is the equivalent TNT at the actual blast source location, in kilograms.

[0101] W1 can be calculated using the following formula:

[0102]

[0103] Among them, t + The actual time of positive pressure action from the explosion source is expressed in seconds (s).

[0104]

[0105] Where i1 is the specific impulse of the actual explosion source, in N·s / m 2 .

[0106]

[0107] Further, calculate according to the following formula

[0108]

[0109] Where ρ0 is the density of the ambient medium at the actual location of the explosion source, and Δρ1 is the density increment of the ambient medium at the actual location of the explosion source.

[0110] Step S602: Calculate the second shock wave pressure of the virtual impact explosion source set;

[0111] Specifically, in this embodiment, the virtual impact sources are distributed at k different locations. It is possible that there are multiple virtual impact sources at a certain location. The shock wave pressure at a certain location can be calculated according to the following formula:

[0112]

[0113] Among them, P2 k ΔP represents the shock wave pressure superimposed by the virtual impact sources at the k-th location; N represents the total number of virtual impact sources at the k-th location; P0' represents the environmental pressure of the virtual impact sources at the k-th location; ΔP 2m ρ represents the peak overpressure generated by the m-th virtual impact source exploding alone at the k-th location, expressed in MPa. 2m Let m be the shock wave density of the m-th virtual impact source when it explodes alone at the k-th location; Let be the velocity of the shock wave generated by the m-th virtual impact blast source detonating alone at the k-th location; Let be the density increment of the shock wave generated by the m-th virtual impact source exploding alone at the k-th location.

[0114] ΔP is calculated using the following formula. 2m :

[0115]

[0116] Among them, Z 2m The first proportional distance to the m-th virtual impact source at the k-th location can be calculated using the following formula:

[0117]

[0118] Among them, R 2m W is the distance between the observation point and the virtual impact source at the k-th location, in meters. 2m The equivalent TNT yield of the m-th virtual impact blast source exploding alone at the k-th location is expressed in kg.

[0119] W1 can be calculated using the following formula:

[0120]

[0121] Among them, t + ' is the time of positive pressure effect when the m-th virtual influence explosion source explodes alone at the k-th location, in seconds.

[0122]

[0123] Where i2 is the specific impulse of the m-th virtual impact source detonating alone at the k-th position, in N·s / m. 2 .

[0124]

[0125] Further, calculate according to the following formula

[0126]

[0127] Where, ρ 0k ' represents the density of the environmental medium at the k-th location, Δρ 2k This represents the density increment of the environmental medium at the k-th location.

[0128] The second shock wave pressure of the virtual impact explosion source set is calculated using the following formula:

[0129]

[0130] Where M represents the number of virtual impact explosion source locations.

[0131] Step S603: Determine the shock wave pressure at the observation point based on the sum of the first shock wave pressure and the second shock wave pressure.

[0132] Specifically, in this embodiment, the shock wave pressure at a certain observation point after multiple reflections from the actual explosion source is obtained by adding the first shock wave pressure P1 and the second shock wave pressure P2. Calculating the shock wave pressure at each moment during the explosion process allows for the plotting of the shock wave pressure time history curve. Differentiating the shock wave pressure at each moment (i.e., the derivative of the shock wave pressure with respect to time) allows for the plotting of the impulse time history curve. Analyzing the pressure time history curve and the impulse time history curve reveals the evolution of shock wave pressure and impulse over time during the explosion process.

[0133] Another embodiment of the present invention provides a processor configured to execute the shock wave pressure calculation method for an explosion inside a sealed structure as described in the above embodiments.

[0134] Another embodiment of the present invention provides a machine-readable storage medium storing instructions for causing a machine to execute the shock wave pressure calculation method for an explosion inside a sealed structure as described in the above embodiments.

[0135] Those skilled in the art will understand that embodiments of this application can be provided as methods, systems, or computer program products. Therefore, this application can take the form of a completely hardware embodiment, a completely software embodiment, or an embodiment combining software and hardware aspects. Furthermore, this application can take the form of a computer program product embodied on one or more computer-usable storage media (including but not limited to disk storage, CD-ROM, optical storage, etc.) containing computer-usable program code.

[0136] This application is described with reference to flowchart illustrations and / or block diagrams of methods, apparatus (systems), and computer program products according to embodiments of this application. It will be understood that each block of the flowchart illustrations and / or block diagrams, and combinations of blocks in the flowchart illustrations and / or block diagrams, can be implemented by computer program instructions. These computer program instructions can be provided to a processor of a general-purpose computer, special-purpose computer, embedded processor, or other programmable data processing apparatus to produce a machine, such that the instructions, which execute via the processor of the computer or other programmable data processing apparatus, generate instructions for implementing the flowchart... Figure 1 One or more processes and / or boxes Figure 1 A device that provides the functions specified in one or more boxes.

[0137] These computer program instructions may also be stored in a computer-readable storage medium that can direct a computer or other programmable data processing device to function in a particular manner, such that the instructions stored in the computer-readable storage medium produce an article of manufacture including instruction means, which are implemented in a process Figure 1 One or more processes and / or boxesFigure 1 The function specified in one or more boxes.

[0138] These computer program instructions may also be loaded onto a computer or other programmable data processing equipment to cause a series of operational steps to be performed on the computer or other programmable equipment to produce a computer-implemented process, thereby providing instructions that execute on the computer or other programmable equipment for implementing the process. Figure 1 One or more processes and / or boxes Figure 1 The steps of the function specified in one or more boxes.

[0139] In a typical configuration, a computing device includes one or more processors (CPU), input / output interfaces, network interfaces, and memory.

[0140] Memory may include non-persistent memory in computer-readable media, such as random access memory (RAM) and / or non-volatile memory, such as read-only memory (ROM) or flash RAM. Memory is an example of computer-readable media.

[0141] Computer-readable media includes both permanent and non-permanent, removable and non-removable media that can store information using any method or technology. Information can be computer-readable instructions, data structures, modules of programs, or other data. Examples of computer storage media include, but are not limited to, phase-change memory (PRAM), static random access memory (SRAM), dynamic random access memory (DRAM), other types of random access memory (RAM), read-only memory (ROM), electrically erasable programmable read-only memory (EEPROM), flash memory or other memory technologies, CD-ROM, digital versatile optical disc (DVD) or other optical storage, magnetic tape, magnetic magnetic disk storage or other magnetic storage devices, or any other non-transferable medium that can be used to store information accessible by a computing device. As defined herein, computer-readable media does not include transient computer-readable media, such as modulated data signals and carrier waves.

[0142] It should also be noted that the terms "comprising," "including," or any other variations thereof are intended to cover non-exclusive inclusion, such that a process, method, article, or apparatus that comprises a list of elements includes not only those elements but also other elements not expressly listed, or elements inherent to such process, method, article, or apparatus. Unless otherwise specified, an element defined by the phrase "comprising one..." does not exclude the presence of other identical elements in the process, method, article, or apparatus that includes that element.

[0143] The above are merely embodiments of this application and are not intended to limit the scope of this application. Various modifications and variations can be made to this application by those skilled in the art. Any modifications, equivalent substitutions, improvements, etc., made within the spirit and principles of this application should be included within the scope of the claims of this application.

Claims

1. A method for calculating shock wave pressure in an explosion inside a sealed structure, characterized in that, The sealed structure is a tube gallery structure with a rectangular cross-section, and the method for calculating the shock wave pressure includes: The actual explosion source is mirrored using each wall surface of the sealed structure as a mirror, and the set of mirrored explosion sources is determined according to the number of actual explosion sources. Determine the coordinates of the observation point for the explosion inside the sealed structure; The influence domain segmentation plane corresponding to the secondary mirror explosion source is determined based on the coordinates of the observation point and the coordinates of the secondary mirror explosion sources in the mirror explosion source set. The coordinates of the intersection point are determined based on the influence domain segmentation plane and the mirror plane corresponding to the secondary mirror explosion source. The set of virtual impact explosion sources is determined based on the coordinates of the intersection point and the coordinates of the observation point; The shock wave pressure at the observation point during an explosion inside the sealed structure is calculated based on the actual explosion source and the set of virtual impact explosion sources; wherein... The step of determining the influence domain segmentation plane corresponding to the secondary mirrored explosion source based on the coordinates of the observation point and the coordinates of the secondary mirrored explosion sources in the mirrored explosion source set includes: Based on the coordinates of the observation point and the coordinates of one of the secondary mirrored explosion sources in the mirrored explosion source set, determine the coordinate point corresponding to the coordinates of the secondary mirrored explosion source and the straight line where the observation point is located. The straight line is the influence domain dividing line. The plane that passes through the dividing line of the influence domain and is perpendicular to the vertical section of the pipe gallery structure is defined as the dividing plane of the influence domain.

2. The method for calculating shock wave pressure in an explosion inside a sealed structure according to claim 1, characterized in that, The step of mirroring the actual explosion source using each wall surface of the sealed structure as a mirror and determining the mirrored explosion source set based on the number of actual explosion sources includes: Determine the number of times the actual explosion source is mirrored; The set of mirrored explosion sources is determined based on the number of mirrored explosions and the number of actual explosion sources.

3. The method for calculating shock wave pressure in an explosion inside a sealed structure according to claim 2, characterized in that, The number of mirroring operations is 2.

4. The method for calculating shock wave pressure in an explosion inside a sealed structure according to claim 3, characterized in that, Determining the set of virtual impact sources based on the coordinates of the intersection point and the coordinates of the observation point includes: The first coordinate value and the second coordinate value are determined based on the coordinates of the observation point and the coordinates of the secondary mirror explosion source; The first preset range and the second preset range are determined based on the set of size parameters of the sealed structure; Determine whether the first coordinate value belongs to the first preset range; Determine whether the second coordinate value belongs to the second preset range; If the first coordinate value belongs to the first preset range and the second coordinate value belongs to the second preset range, the secondary mirror explosion source is determined to be the virtual influence explosion source; The set of virtual impact sources is determined based on the virtual impact sources.

5. The method for calculating shock wave pressure in an explosion inside a sealed structure according to claim 4, characterized in that, Determining the first preset range and the second preset range based on the set of size parameters of the sealed structure includes: The first preset range is determined based on the first parameter in the set of size parameters; The second preset range is determined based on the second parameter in the set of size parameters, wherein the first parameter is a parameter belonging to the same coordinate direction as the first coordinate value, and the second parameter is a parameter belonging to the same coordinate direction as the second coordinate value.

6. The method for calculating shock wave pressure in an explosion inside a sealed structure according to claim 1, characterized in that, The calculation of the shock wave pressure at the observation point during an explosion inside the sealed structure based on the actual explosion source and the set of virtual impact explosion sources includes: Calculate the first shock wave pressure of the actual explosion source; Calculate the second shock wave pressure of the virtual impact explosion source set; The shock wave pressure at the observation point is determined based on the sum of the first shock wave pressure and the second shock wave pressure.

7. A processor, characterized in that, It is configured to perform the shock wave pressure calculation method for an explosion inside a sealed structure as described in any one of claims 1 to 6.

8. A machine-readable storage medium storing instructions thereon, characterized in that, The instructions are used to cause the machine to execute the shock wave pressure calculation method for an explosion inside a sealed structure according to any one of claims 1 to 6.

Images