A method for predicting the detonation state of centrosymmetric PBX9404 charge

Abstract

The present application discloses a method for predicting the detonation state of a centrosymmetric PBX9404 charge, and relates to the field of detonation technology. The method uses the distance between the detonation wavefront and the detonation point to normalize the distance between different positions and the detonation point of the charge into a normalized spacing, and provides a fitting function between the normalized state parameters and the normalized spacing under different PBX9404 detonation conditions. The normalized state parameters can be obtained by bringing the normalized spacing at different positions into the fitting function, and then the detonation state parameters at different positions can be obtained by combining the normalized state parameters and the CJ state quantity of the PBX9404. The method can realize the rapid and accurate prediction of the detonation state in the centrosymmetric PBX9404 charge scenario, has excellent calculation accuracy, stability and engineering applicability, and can be used for both near-field explosion prediction and far-field explosion prediction.

Description

Technical Field

[0001] The present application relates to the field of detonation technology, and in particular to a method for predicting the detonation state of a centrosymmetric PBX9404 charge. Background Art

[0002] The explosion phenomenon has a wide range of applications in military and civilian fields. Whether the explosive bag explodes in water or in the air, its chemical reaction process includes complex phenomena such as flame, deflagration, and detonation. The detonation process mainly includes two basic stages: initiation and propagation. The explosive loses stability and explodes after being affected by the outside world. This process from a stable state to a lost stable state is called the initiation process of the explosive. The propagation of the detonation wave is caused by the impact compression of the shock wave on the explosive. When the intensity of the shock wave exceeds a certain value, it will cause a chemical reaction of the explosive material behind the shock wave front, and the energy released by the chemical reaction is used to support the continued propagation of the shock wave, thus forming the detonation process of the explosive.

[0003] In order to evaluate and predict the detonation damage effect, it is often necessary to predict the detonation state. There are currently two main methods: (1) Introducing a complex detonation calculation model that can predict in detail the detonation chemical reaction process starting from the detonation point, mainly involving the generation and transmission process of the detonation wave. Common detonation calculation models include the "ignition-growth" reaction model proposed by Lee and Tarver. The use of detonation calculation models can realize a quantitative and refined prediction method for the detonation state with high calculation accuracy. Since it can describe the generation, development, and propagation process of the detonation wave, as well as the coupling effect between the detonation wave and the surrounding flow field, it can be applied in scenarios considering the near-field effect of the explosion. However, the various existing detonation calculation models are often very complex, requiring a large number of input parameters, and it is also difficult to determine the specific values ​​of these parameters. The calculation complexity and amount of calculation are large, and it takes a long time. (2) Another method ignores the specific detonation process of the explosive, that is, it does not consider the details of the detonation of the explosive, and simplifies the impact of the explosive on the surrounding environment into the effect of a mass of high-temperature gas expansion. This method assumes that the detonation of the explosive is completed instantaneously, and the detonation product gas becomes a mass of high-temperature, high-density, and high-pressure gas, thereby constructing a simplified instantaneous detonation model. Compared with the first method, this method can simplify the calculation process and is easy to implement in engineering. However, since this method does not consider the details of the detonation process, it can only be used for qualitative engineering estimation in scenarios where the far-field effects of the explosion are considered. When the near-field effects of the explosion need to be considered or quantitative high-precision predictions are required, this method cannot be applied. It can be seen from this that the two existing commonly used detonation state prediction methods have their own defects and cannot meet the needs of actual academic research and engineering application scenarios. Summary of the invention

[0004] In view of the above problems and technical requirements, this application proposes a method for predicting the detonation state of a centrally symmetrical PBX9404 charge. The technical solution of this application is as follows:

[0005] A method for predicting the detonation state of a centrosymmetric PBX9404 charge, the method comprising:

[0006] The CJ state quantity of PBX9404 is used to construct the normalized state parameters at any position Q in the calculation domain where the drug package is located. The drug package is filled with PBX9404, the charge structure of the drug package is a centrally symmetrical structure, and the detonation point is located at the center of the drug package; the normalized state parameters at position Q include the normalized density Normalized speed Normalized pressure and normalized internal energy The CJ state variables of PBX9404 include ρ CJ =2488kg / m 3 、u CJ =2285m / s, p CJ =37GPa, e CJ =8.15×10 6 J / kg; the detonation state parameters at position Q include density ρ(r), velocity u(r), pressure p(r) and internal energy e(r);

[0007] According to the detonation velocity D of PBX9404, the normalized spacing at any position Q in the calculation domain where the charge is located at any time t is constructed. Where r is the distance r between position Q and the detonation point of the charge, and R0 = D × t is the distance between the detonation wave front and the detonation point at time t;

[0008] According to the detonation characteristics under the PBX9404 detonation condition type, each normalized state parameter at position Q at any time t is established With normalized spacing The fitting function relationship between Among them, the current PBX9404 detonation condition type corresponds to the charge structure of the explosive package, and the boundary distance r0 of the detonation static zone and the values ​​of the parameters a0, b0, b1, b2, b3, and b4 in the fitting function relationship of each normalized state parameter are all related to the detonation characteristics under the current PBX9404 detonation condition type;

[0009] Using each normalized state parameter With normalized spacing The fitting function relationship between them can be used to predict the detonation state parameters at different positions in the calculation domain where the charge is located at different times.

[0010] Its further technical solution is that when the charge structure of the explosive package is a one-dimensional plane charge structure, the current PBX9404 detonation condition type is a PBX9404 one-dimensional plane detonation condition;

[0011] The boundary distance of the detonation static zone corresponding to the one-dimensional plane detonation condition of PBX9404 is r0=0.5;

[0012] Normalized density corresponding to the one-dimensional plane detonation condition of PBX9404 The parameter a0 in the fitting function relationship is 0.663, the normalized speed The parameters a0 in the fitting function relationship are 0 and the normalized pressure The parameters a0 in the fitting function relationship are 0.277 and the normalized internal energy The parameter a0 in the fitting function relationship is 1.491.

[0013] Its further technical solution is that when the charge structure of the explosive package is a one-dimensional plane charge structure, the current PBX9404 detonation condition type is the PBX9404 one-dimensional plane detonation condition; each normalized state parameter corresponding to the PBX9404 one-dimensional plane detonation condition is The fitting function relationship within the range is:

[0014]

[0015] Its further technical solution is that when the charge structure of the explosive package is a two-dimensional circular charge structure, the current PBX9404 detonation condition type is a PBX9404 one-dimensional cylindrical detonation condition;

[0016] The boundary distance of the detonation static zone corresponding to the one-dimensional cylindrical detonation condition of PBX9404 is r0=0.47;

[0017] Normalized density corresponding to the one-dimensional cylindrical detonation condition of PBX9404 The parameter a0 in the fitting function relationship is 0.627, the normalized speed The parameters a0 in the fitting function relationship are 0 and the normalized pressure The parameters a0 in the fitting function relationship are 0.232 and the normalized internal energy The parameter a0 in the fitting function relationship is 0.452.

[0018] Its further technical solution is that when the charge structure of the explosive package is a two-dimensional circular charge structure, the current PBX9404 detonation condition type is a PBX9404 one-dimensional cylindrical detonation condition; each normalized state parameter corresponding to the PBX9404 one-dimensional cylindrical detonation condition is The fitting function relationship within the range is:

[0019]

[0020] Its further technical solution is that when the charge structure of the explosive package is a three-dimensional spherical charge structure, the current PBX9404 detonation condition type is a PBX9404 one-dimensional spherical detonation condition;

[0021] The boundary distance of the detonation static zone corresponding to the one-dimensional spherical detonation condition of PBX9404 is r0=0.46;

[0022] Normalized density corresponding to the one-dimensional spherical detonation condition of PBX9404 The parameter a0 in the fitting function relationship is 0.606, the normalized speed The parameters a0 in the fitting function relationship are 0 and the normalized pressure The parameters a0 in the fitting function relationship are 0.208 and the normalized internal energy The parameter a0 in the fitting function relationship is 0.429.

[0023] Its further technical solution is that when the charge structure of the explosive package is a three-dimensional spherical charge structure, the current PBX9404 detonation condition type is a PBX9404 one-dimensional spherical detonation condition; each normalized state parameter corresponding to the PBX9404 one-dimensional spherical detonation condition is The fitting function relationship within the range is:

[0024]

[0025] A further technical solution is that the explosion velocity D of PBX9404 is 8800 m / s.

[0026] The beneficial technical effects of this application are:

[0027] The present application discloses a method for predicting the detonation state of a centrosymmetric PBX9404 charge. The method targets a classic explosion scenario of a common explosive type such as PBX9404 in a common charge structure of a centrosymmetric charge. By analyzing the propagation characteristics of the detonation wave after the detonation of the centrosymmetric PBX9404 charge, it is found that the detonation wave has good similarity inside the explosive package, thereby providing a high-precision instantaneous detonation model for the explosion scenario through the constructed normalized variables. The fitting formula provided by the present application can directly and quickly calculate the normalized state parameters at different positions, and then the detonation state parameters at different positions can be obtained by combining the CJ state quantity of PBX9404. The fitting function provided in the present application is established by taking into account the details of the detonation process, and therefore has high prediction accuracy. The fitting function is simple in form, so the amount of calculation is small and the calculation time is short, which makes the prediction method of the present application have good performance in terms of prediction accuracy and prediction speed. It can be used for near-field explosion prediction as well as far-field explosion prediction, and has better calculation accuracy, stability and engineering applicability. It has important significance and value both in academic research and engineering application.

[0028] This application also specifically gives the parameter values ​​in the fitting function of each normalized state parameter under different PBX9404 detonation conditions, which improves the calculation accuracy and efficiency of charge detonation, and can provide key technical support for shock wave and bubble movement and load propagation under near-field and far-field explosions of PBX9404 charges in air and underwater environments. BRIEF DESCRIPTION OF THE DRAWINGS

[0029] Figure 1 It is a method flow chart of a detonation state prediction method according to an embodiment of the present application.

[0030] Figure 2 These are the different normalized spacings under three different PBX9404 detonation conditions. The normalized density under Graph of the curve.

[0031] Figure 3 These are the different normalized spacings under three different PBX9404 detonation conditions. Normalized speed under Graph of the curve.

[0032] Figure 4 These are the different normalized spacings under three different PBX9404 detonation conditions. Normalized pressure under Graph of the curve.

[0033] Figure 5 These are the different normalized spacings under three different PBX9404 detonation conditions. Normalized internal energy Graph of the curve.

[0034] Figure 6 It is a distribution curve diagram of detonation state parameters at different positions of the detonation point at different times in a detonation example of a PBX9404 one-dimensional plane detonation condition.

[0035] Figure 7 This is a density distribution cloud diagram at different times and different positions from the detonation point in another detonation example of the PBX9404 one-dimensional cylindrical detonation condition.

[0036] Figure 8 yes Figure 7 In the detonation example, the distribution cloud diagram of the velocity component along the X direction at different positions and at different times relative to the detonation point.

[0037] Fig. 9 yes Figure 7 In the detonation example, the distribution cloud diagram of the component of the velocity along the Y direction at different positions and at different times relative to the detonation point.

[0038] Fig.10 yes Figure 7 In the detonation example, the pressure distribution cloud diagram at different positions and at different times relative to the detonation point.

[0039] Fig.11 yes Figure 7 In the detonation example, the internal energy distribution cloud diagram at different positions and at different times relative to the detonation point. DETAILED DESCRIPTION

[0040] The specific implementation of the present application is further described below in conjunction with the accompanying drawings.

[0041] This application discloses a method for predicting the detonation state of a centrally symmetrical PBX9404 charge. Figure 1 The method flow chart is shown. The method is aimed at the application scenario of the centrally symmetrical PBX9404 charge, and can be used to calculate and predict the detonation state of the centrally symmetrical structure charge bag filled with PBX9404 after detonation. The charge structure of the charge bag targeted by this application is a centrally symmetrical structure, the charge bag is filled with PBX9404, and the detonation point of the charge bag is located at the center of the charge bag.

[0042] There are three main types of charge structures for the centrally symmetrical charge, corresponding to different PBX9404 detonation conditions:

[0043] (1) When the charge structure of the charge is a one-dimensional plane charge structure, the current PBX9404 detonation condition type is the PBX9404 one-dimensional plane detonation condition.

[0044] (2) When the charge structure of the explosive package is a two-dimensional circular charge structure, the current PBX9404 detonation condition type is the PBX9404 one-dimensional cylindrical detonation condition.

[0045] (3) When the charge structure of the charge is a three-dimensional spherical charge structure, the current PBX9404 detonation condition type is the PBX9404 one-dimensional spherical detonation condition.

[0046] When predicting the detonation state of the centrally symmetrical PBX9404 charge, it is necessary to obtain the detonation state parameters of any position Q at any time t after detonation. The detonation state parameters include density, velocity, pressure and internal energy. Regardless of the centrally symmetrical charge structure mentioned above, the internal detonation wave structure of the charge after central detonation has obvious structural similarities at different times, so it can be normalized. Therefore, this application introduces the normalized spacing and normalized state parameters, the method comprises the following steps:

[0047] Using the CJ state quantity of PBX9404, the detonation state parameters at any position Q in the calculation domain where the charge is located are normalized to construct the normalized state parameters at position Q. The normalized state parameters at position Q include the normalized density Normalized speed Normalized pressure and normalized internal energy and: Among them, the detonation state parameters at position Q include density ρ(r), velocity u(r), pressure p(r), and internal energy e(r). The CJ state quantities of PBX9404 include ρ CJ =2488kg / m 3 、u CJ =2285m / s, p CJ =37GPa, e CJ =8.15×10 6 J / kg.

[0048] Using the distance R0=D×t between the detonation wave front and the detonation point at the current time t, the distance r between any position Q in the calculation domain where the explosive package is located and the detonation point of the explosive package at any time t is normalized to obtain the normalized distance between the position Q and the detonation point at any time t Normalized spacing is a dimensionless parameter and D is the detonation velocity of PBX9404 and D=8800m / s.

[0049] This application studies the detonation characteristics under different PBX9404 detonation conditions, divides the detonation zone into a detonation static zone close to the detonation point and a detonation sparse zone far from the detonation point, and directly performs function fitting on the normalized state parameters and normalized spacing of the positions of the static zone and the sparse zone, thereby establishing each normalized state parameter of the position Q at any time t With normalized spacing The fitting function relationship between them, where the normalized state parameter is the normalized density Normalized speed Normalized pressure and normalized internal energy Regardless of the PBX9404 detonation condition type under the above-mentioned charge structure, each normalized state parameter With normalized spacing The functional form of the fitting function relationship between is:

[0050]

[0051] Among them, the current PBX9404 detonation condition type corresponds to the charge structure of the explosive package, the boundary distance r0 of the detonation static zone, and the values ​​of parameters a0, b0, b1, b2, b3, and b4 in the fitting function relationship of each normalized state parameter are all related to the detonation characteristics under the current PBX9404 detonation condition type.

[0052] The following are the three types of charge structures:

[0053] (1) When the charge structure of the explosive package is a one-dimensional plane charge structure, the current PBX9404 detonation condition type is the PBX9404 one-dimensional plane detonation condition. At this time:

[0054] The boundary distance r0 of the detonation static zone corresponding to the one-dimensional plane detonation condition of PBX9404 is 0.5.

[0055] Normalized density corresponding to the one-dimensional plane detonation condition of PBX9404 The parameter a0 in the fitting function relationship is 0.663, the normalized speed The parameters a0 in the fitting function relationship are 0 and the normalized pressure The parameters a0 in the fitting function relationship are 0.277 and the normalized internal energy The parameter a0 in the fitting function relationship is 1.491.

[0056] The normalized state parameters corresponding to the one-dimensional plane detonation condition of PBX9404 are The fitting function relationship within the range is:

[0057]

[0058] It can be concluded that under the one-dimensional plane detonation condition of PBX9404, the four different normalized state parameters and the normalized spacing at different positions are shown in Table 1:

[0059] Table 1. Fitting function relationship of PBX9404 under one-dimensional plane detonation condition

[0060]

[0061] (2) When the charge structure of the explosive package is a two-dimensional circular charge structure, the current PBX9404 detonation condition type is the PBX9404 one-dimensional cylindrical detonation condition. At this time:

[0062] The boundary distance of the detonation static zone corresponding to the one-dimensional cylindrical detonation condition of PBX9404 is r0=0.47;

[0063] Normalized density corresponding to the one-dimensional cylindrical detonation condition of PBX9404 The parameter a0 in the fitting function relationship is 0.627, the normalized speed The parameters a0 in the fitting function relationship are 0 and the normalized pressure The parameters a0 in the fitting function relationship are 0.232 and the normalized internal energy The parameter a0 in the fitting function relationship is 0.452.

[0064] The normalized state parameters corresponding to the one-dimensional cylindrical detonation condition of PBX9404 are The fitting function relationship within the range is:

[0065]

[0066] It can be concluded that under the one-dimensional cylindrical detonation condition of PBX9404, the four different normalized state parameters and the normalized spacing at different positions are shown in Table 2 below:

[0067] Table 2. Fitting function relationship of PBX9404 under one-dimensional cylindrical detonation condition

[0068]

[0069] (3) When the charge structure of the explosive package is a three-dimensional spherical charge structure, the current PBX9404 detonation condition type is the PBX9404 one-dimensional spherical detonation condition. At this time:

[0070] The boundary distance r0 of the detonation static zone corresponding to the one-dimensional spherical detonation condition of PBX9404 is 0.46.

[0071] Normalized density corresponding to the one-dimensional spherical detonation condition of PBX9404 The parameter a0 in the fitting function relationship is 0.606, the normalized speed The parameters a0 in the fitting function relationship are 0 and the normalized pressure The parameters a0 in the fitting function relationship are 0.208 and the normalized internal energy The parameter a0 in the fitting function relationship is 0.429.

[0072] The normalized state parameters corresponding to the one-dimensional spherical detonation condition of PBX9404 are The fitting function relationship within the range is:

[0073]

[0074] It can be concluded that under the one-dimensional spherical detonation condition of PBX9404, the four different normalized state parameters and the normalized spacing at different positions are shown in Table 2 below:

[0075] Table 3. Fitting function relationship of PBX9404 under one-dimensional spherical detonation condition

[0076]

[0077]

[0078] Normalized density under the above three different PBX9404 detonation conditions Normalized spacing over position The curve diagram is as follows Figure 2 The normalized velocity under the above three different PBX9404 detonation conditions Normalized spacing over position The curve diagram is as follows Figure 3 Normalized pressure under the above three different PBX9404 detonation conditions Normalized spacing over position The curve diagram is as follows Figure 4 The normalized internal energy of the above three different PBX9404 detonation conditions is shown in Figure 2. Normalized spacing over position The curve diagram is as follows Figure 5 shown.

[0079] Then, each normalized state parameter is used With normalized spacing The fitting function relationship between them predicts the detonation state parameters at different positions in the calculation domain where the charge is located at different times. Specifically, after determining the detonation condition type of PBX9404, the normalized state parameters can be determined. With normalized spacing The parameters in the fitting function relationship between the two are then used as the normalized spacing of any position Q in the calculation domain where the drug package is located at any time t Substituting into the fitting function relationship, we can obtain the detonation state parameters at position Q at time t.

[0080] For example, in one example, the charge structure of the explosive package is a one-dimensional planar charge structure. The detonation velocity of PBX9404 is D = 8800 m / s. Therefore, it can be known that the distance R0 between the detonation wave front and the detonation point at t = 1 μs after detonation is 8.8 mm, the distance R0 between the detonation wave front and the detonation point at t = 3 μs after detonation is 26.4 mm, the distance R0 between the detonation wave front and the detonation point at t = 5 μs after detonation is 44 mm, the distance R0 between the detonation wave front and the detonation point at t = 7 μs after detonation is 61.6 mm, and the distance R0 between the detonation wave front and the detonation point at t = 9 μs after detonation is 79.2 mm. Combined with the CJ state quantity ρ of PBX9404 CJ =2488kg / m 3 、u CJ =2285m / s, p CJ =37GPa, e CJ =8.15×10 6 J / kg. Based on the fitting function in Table 1, the detonation state parameters at different times at different positions within the range of 0.1m from the detonation point can be obtained. Figure 6 As shown in (a), (b), (c), and (d), the results are basically consistent with the prediction results obtained by the existing detonation calculation model, indicating that the prediction method of the present application has high prediction accuracy, and after verification by different examples, the prediction method of the present application has good prediction accuracy under various PBX9404 detonation conditions.

[0081] For example, in another example, the charge structure of the explosive package is a two-dimensional circular charge structure, the detonation point is located at the origin of the coordinate system, and the calculation domain is [-0.20, 0.20] × [-0.20, 0.20] m 2 , using a uniform grid. Combined with the detonation velocity of PBX9404 being D=8800m / s, it can be seen that the distance R0 between the detonation wave front and the detonation point at t=5μs after detonation is 44mm, the distance R0 between the detonation wave front and the detonation point at t=10μs after detonation is 88mm, the distance R0 between the detonation wave front and the detonation point at t=15μs after detonation is 132mm, and the distance R0 between the detonation wave front and the detonation point at t=20μs after detonation is 176mm. Based on the fitting function in Table 1, the density distribution cloud diagram at different positions in the calculation domain at different times can be obtained as follows Figure 7The distribution cloud diagram of the velocity component along the X direction at different positions in the computational domain at different times is shown in Figure 8 As shown, the distribution cloud diagram of the velocity component along the Y direction at different positions in the calculation domain at different times is as follows Fig. 9 The pressure distribution cloud diagram of different locations in the computational domain at different times is shown in Fig.10 As shown, the distribution cloud diagram of the internal energy at different locations in the calculation domain at different times is as follows Fig.11 shown.

[0082] This application directly provides a quantitative calculation formula for detonation state parameters under different PBX9404 detonation conditions based on the detonation characteristics of the centrally symmetrical PBX9404 charge, so that the detonation process can be quickly predicted by bringing it into the solution. And the quantitative calculation formula for detonation state parameters obtained by fitting has a high prediction accuracy.

[0083] The above is only a preferred embodiment of the present application, and the present application is not limited to the above embodiments. It is understood that other improvements and changes directly derived or associated by those skilled in the art without departing from the spirit and concept of the present application should be considered to be included in the protection scope of the present application.

Claims

1. A method for predicting the detonation state of a centrosymmetric PBX9404 charge, characterized in that: The detonation state prediction method comprises: The detonation state parameters at any position Q in the calculation domain where the charge is located are normalized by using the CJ state quantity of PBX9404 to obtain normalized state parameters. The charge is filled with PBX9404, the charge structure of the charge is centrally symmetrical and the detonation point is located at the center of the charge. The normalized state parameters at position Q include normalized density Normalized speed Normalized pressure and normalized internal energy The CJ state variables of PBX9404 include ρ CJ =2488kg / m 3 、u CJ =2285m / s, p CJ =37GPa, e CJ =8.15×10 6 J / kg; the detonation state parameters at position Q include density ρ(r), velocity u(r), pressure p(r) and internal energy e(r); According to the detonation velocity D of PBX9404, the normalized spacing at any position Q in the calculation domain where the charge is located at any time t is constructed. Where r is the distance r between position Q and the detonation point of the charge, and R0 = D × t is the distance between the detonation wave front and the detonation point at time t; According to the detonation characteristics under the PBX9404 detonation condition type, each normalized state parameter at position Q at any time t is established With normalized spacing The fitting function relationship between Among them, the current PBX9404 detonation condition type corresponds to the charge structure of the explosive package, and the boundary distance r0 of the detonation static zone and the values ​​of the parameters a0, b0, b1, b2, b3, and b4 in the fitting function relationship of each normalized state parameter are all related to the detonation characteristics under the current PBX9404 detonation condition type; Using each normalized state parameter With normalized spacing The fitting function relationship between them can be used to predict the detonation state parameters at different positions in the calculation domain where the charge is located at different times.

2. The detonation state prediction method according to claim 1, characterized in that: When the charge structure of the explosive package is a one-dimensional plane charge structure, the current PBX9404 detonation condition type is a PBX9404 one-dimensional plane detonation condition; The boundary distance of the detonation static zone corresponding to the one-dimensional plane detonation condition of PBX9404 is r0=0.5; Normalized density corresponding to the one-dimensional plane detonation condition of PBX9404 The parameter a0 in the fitting function relationship is 0.663, the normalized speed The parameters a0 in the fitting function relationship are 0 and the normalized pressure The parameters a0 in the fitting function relationship are 0.277 and the normalized internal energy The parameter a0 in the fitting function relationship is 1.

491.

3. The detonation state prediction method according to claim 1, characterized in that: When the charge structure of the explosive package is a one-dimensional plane charge structure, the current PBX9404 detonation condition type is the PBX9404 one-dimensional plane detonation condition; the normalized state parameters corresponding to the PBX9404 one-dimensional plane detonation condition are The fitting function relationship within the range is:

4. The detonation state prediction method according to claim 1, characterized in that: When the charge structure of the explosive package is a two-dimensional circular charge structure, the current PBX9404 detonation condition type is a PBX9404 one-dimensional cylindrical detonation condition; The boundary distance of the detonation static zone corresponding to the one-dimensional cylindrical detonation condition of PBX9404 is r0=0.47; Normalized density corresponding to the one-dimensional cylindrical detonation condition of PBX9404 The parameter a0 in the fitting function relationship is 0.627, the normalized speed The parameters a0 in the fitting function relationship are 0 and the normalized pressure The parameters a0 in the fitting function relationship are 0.232 and the normalized internal energy The parameter a0 in the fitting function relationship is 0.

452.

5. The detonation state prediction method according to claim 1, characterized in that: When the charge structure of the explosive package is a two-dimensional circular charge structure, the current PBX9404 detonation condition type is a PBX9404 one-dimensional cylindrical detonation condition; the normalized state parameters corresponding to the PBX9404 one-dimensional cylindrical detonation condition are The fitting function relationship within the range is:

6. The detonation state prediction method according to claim 1, characterized in that: When the charge structure of the explosive package is a three-dimensional spherical charge structure, the current PBX9404 detonation condition type is a PBX9404 one-dimensional spherical detonation condition; The boundary distance of the detonation static zone corresponding to the one-dimensional spherical detonation condition of PBX9404 is r0=0.46; Normalized density corresponding to the one-dimensional spherical detonation condition of PBX9404 The parameter a0 in the fitting function relationship is 0.606, the normalized speed The parameters a0 in the fitting function relationship are 0 and the normalized pressure The parameters a0 in the fitting function relationship are 0.208 and the normalized internal energy The parameter a0 in the fitting function relationship is 0.

429.

7. The detonation state prediction method according to claim 1, characterized in that: When the charge structure of the explosive package is a three-dimensional spherical charge structure, the current PBX9404 detonation condition type is a PBX9404 one-dimensional spherical detonation condition; The normalized state parameters corresponding to the one-dimensional spherical detonation condition of PBX9404 are The fitting function relationship within the range is:

8. The detonation state prediction method according to claim 1, characterized in that: The explosion speed of PBX9404 is D=8800m / s.

Images