Quantitative prediction method for state of three-dimensional PETN charging instantaneous detonation
Through the quantitative forecast method of instantaneous detonation state of three-dimensional PETN charge, the problem of high computational complexity of existing methods is solved by using normalized processing and CJ detonation parameters, and a fast and quantitative detonation state forecast is achieved, which is suitable for high-precision forecasting in far and near fields.
Patent Information
- Application Number
- CN202510349335.2
- Authority / Receiving Office
- CN · China
- Patent Type
- Applications(China)
- Current Assignee / Owner
- Filing Date
- 2025-03-24
- Publication Date
- 2025-07-11
AI Technical Summary
The existing PETN charge detonation status forecasting methods have problems such as high computational complexity, difficult to determine parameters and time-consuming, and cannot meet the quantitative high-precision forecasting requirements of near-field effects.
The state quantitative forecast method of instantaneous detonation of three-dimensional PETN charge is adopted. By normalizing the distance between the detonation wave front and the starting point, combined with the CJ detonation parameters of PETN, the density, radial velocity, pressure and internal energy per unit mass at any position and at any time are quickly calculated.
实现了快速、定量的爆轰状态预报,计算复杂度和运算量小,预报速度快,精度与复杂爆轰计算模型相当,适用于远近场的高精度预报。
Smart Images
Figure CN120297040A_ABST
Abstract
Description
Technical Field
[0001] The present application relates to the technical field of detonation, in particular to a method for quantitatively predicting the instantaneous detonation state of a three-dimensional PETN charge. Background Art
[0002] Explosion phenomena have wide applications in both military and civilian fields. Whether the explosive charge is detonated 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. When the explosive is subjected to an external action, it loses its stability and explodes. This process from a stable state to a state of instability is called the initiation process of the explosive. The propagation of the detonation wave is caused by the shock wave impacting and compressing 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] PETN (pentaerythritol tetranitrate) is an extremely powerful explosive with strong detonation ability and high mechanical sensitivity. It is the most stable one among nitrate ester explosives and also one of the most important high-energy explosives. To evaluate and predict the explosion damage effect of PETN, it is necessary to predict the detonation state of the explosive charge filled with PETN. Currently, there are mainly two methods: (1) Introduce complex detonation calculation models, which can detail the prediction of the detonation chemical reaction process starting from the initiation 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. Using the detonation calculation model can achieve a quantitative and refined prediction method for the detonation state with high calculation accuracy. Since it can describe the generation, development, propagation process of the detonation wave, and the coupling effect between the detonation wave and the surrounding flow field, it can be applied to scenarios considering the near-field effect of the explosion. However, various existing detonation calculation models are often very complex, require a large number of input parameters, and it is also difficult to determine the specific values of these parameters. The calculation complexity and computational workload are large, and it takes a long time. (2) Another method ignores the specific detonation process of the explosive, that is, does not consider the details of the explosive detonation, and simplifies the impact of the explosive on the surrounding environment as the effect generated by the expansion of a mass of high-temperature gas. This method believes that the explosive detonation is completed instantaneously, and the detonation product gas becomes a mass of high-temperature, high-density, and high-pressure gas, thus constructing a simplified instantaneous detonation model. This method can simplify the calculation process compared with the first method 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 considering the far-field effect of the explosion. When the near-field effect of the explosion needs to be considered or a quantitative and high-precision prediction is required, this method is not applicable. It can be seen from this that the two commonly used existing methods for predicting the detonation state have their own defects and cannot meet the requirements in actual academic research and engineering application scenarios. Summary of the Invention
[0004] In view of the above problems and technical requirements, the present application proposes a method for quantitatively predicting the state of instantaneous detonation of a three-dimensional PETN charge. The technical solution of the present application is as follows:
[0005] A method for quantitatively predicting the state of instantaneous detonation of a three-dimensional PETN charge, the state quantitative prediction method comprising:
[0006] Normalize the distance r between any position Q in the calculation domain and the initiation point by the detonation wave front at any time t after the central initiation of the three-dimensional spherical PETN charge to obtain the normalized spacing of position Q at time t
[0007] According to Calculate the normalized density of position Q at time t Combined with the CJ detonation parameters of PETN, the density prediction result ρ at position Q at time t is obtained. t (r);
[0008] According to The normalized radial velocity at position Q at time t is calculated. Combined with the CJ detonation parameters of PETN, the radial velocity prediction result u at position Q at time t is obtained. t (r);
[0009] According to The normalized pressure at position Q at time t is calculated. Combined with the CJ detonation parameters of PETN, the pressure prediction result p at position Q at time t is obtained. t (r);
[0010] According to The normalized internal energy per unit mass at position Q at time t is calculated. Combined with the CJ detonation parameters of PETN, the internal energy per unit mass prediction result e at position Q at time t is obtained. t (r).
[0011] A further technical solution is to obtain the normalized spacing at position Q at time t. It includes:
[0012] Based on the detonation velocity D = 8050 m / s of PETN, the distance R0(t) = D × t between the detonation wave front and the initiation point at any time t after the initiation of the three-dimensional spherical PETN charge is determined, and after normalization, the normalized spacing at position Q at time t is obtained.
[0013] A further technical solution is that the CJ detonation parameters of PETN include the detonation product density ρ of PETN. CJ = 2388 kg / m 3 、the detonation product radial velocity u CJ = 2100 m / s、the detonation pressure p CJ = 30.42 GPa and the internal energy per unit mass e of the detonation product. CJ = 10.38 J / kg;
[0014] Obtaining the density prediction result ρ at position Q at time t. t (r), the radial velocity prediction result u t (r), the pressure prediction result p t (r) and the internal energy per unit mass prediction result e t (r) includes calculation according to Calculation.
[0015] The beneficial technical effects of the present application are as follows:
[0016] The present application discloses a method for quantitatively predicting the instantaneous detonation state of a three-dimensional PETN charge. This method constructs a dimensionless normalized distance between any position in the computational domain and the initiation point at any time after initiation using the distance between the detonation wave front and the initiation point, and provides a quantitative prediction function based on the normalized distance. By using the quantitative prediction function provided by the present application and combining with the CJ detonation parameters of PETN, the prediction results of density, radial velocity, pressure, and internal energy per unit mass at any position in the computational domain at any time can be quickly and quantitatively calculated, without the need to introduce a complex detonation calculation model. The calculation complexity and computational amount are relatively small, and the prediction speed is relatively fast. Moreover, the accuracy of the quantitative prediction can reach the same level as that of the detonation calculation model. Therefore, the method of the present application can be used for quantitative and high-precision prediction of the far and near fields in the central initiation scenario of a three-dimensional PETN charge. Brief Description of the Drawings
[0017] Figure 1 is a schematic flowchart of the state quantitative prediction method according to an embodiment of the present application.
[0018] Figure 2 is a comparison diagram of the time history curves of the prediction results obtained by using the method of the present application and the calculation results of the AUTODYN software in a calculation example.
[0019] Figure 3 is a distribution cloud map of the density prediction results at different times in the target prediction area after central initiation quantitatively predicted by using the state quantitative prediction method of the present application in a simulation example.
[0020] Figure 4 is Figure 3 a distribution cloud map of the prediction results of the X component of the radial velocity at different times in the target prediction area after central initiation quantitatively predicted by using the state quantitative prediction method of the present application in a simulation example.
[0021] Figure 5 is Figure 3 a distribution cloud map of the prediction results of the Y component of the radial velocity at different times in the target prediction area after central initiation quantitatively predicted by using the state quantitative prediction method of the present application in a simulation example.
[0022] Figure 6 is Figure 3 a distribution cloud map of the prediction results of the Z component of the radial velocity at different times in the target prediction area after central initiation quantitatively predicted by using the state quantitative prediction method of the present application in a simulation example.
[0023] Figure 7 is Figure 3The distribution cloud map of the pressure prediction results at different times within the target prediction area after central initiation, quantitatively predicted using the state quantitative prediction method of the present application in the simulation example.
[0024] Figure 8 Yes Figure 3 The distribution cloud map of the internal energy per unit mass prediction results at different times within the target prediction area after central initiation, quantitatively predicted using the state quantitative prediction method of the present application in the simulation example. Detailed implementation manners
[0025] The following further describes the detailed implementation manners of the present application with reference to the accompanying drawings.
[0026] The present application discloses a state quantitative prediction method for the instantaneous detonation of a three-dimensional PETN charge. This method is applicable to the application scenario of central initiation of a three-dimensional spherical PETN charge, where central initiation means the initiation point is located at the center of the three-dimensional spherical PETN charge. This method can be used to quickly and quantitatively predict the detonation state of a three-dimensional spherical charge structure filled with PETN after central initiation, so as to quickly and quantitatively calculate the predicted results of the state parameters at any position Q in the calculation domain at any time t after central initiation. The quantitative prediction method includes:
[0027] Please refer to Figure 1 the schematic diagram of the method flow shown, and first normalize the distance r between any position Q in the calculation domain and the initiation point using the detonation wave front at any time t after the central initiation of the three-dimensional spherical PETN charge to obtain the normalized spacing of position Q at time t after central initiation Specifically, according to the detonation velocity D = 8050 m / s of PETN, determine the distance R0(t) = D × t between the detonation wave front and the initiation point at any time t after the central initiation of the three-dimensional spherical PETN charge, and then normalize to obtain the normalized spacing of position Q at time t
[0028] Based on the normalized spacing of any position Q at time t The state parameters of position Q predicted by the state quantitative prediction method provided by the present application at time t include the density prediction result ρ t (r), the radial velocity prediction result u t (r), the pressure prediction result p t (r), and the internal energy per unit mass prediction result e t (r), which are introduced as follows:
[0029] 1. Density prediction result ρ t (r)
[0030] First, calculate according to to obtain the normalized density of position Q at time t
[0031] Then, combined with the CJ detonation parameters of PETN, the density prediction result ρ at position Q at time t is obtained: t (r), the CJ detonation parameters of PETN used in this step include the detonation product density ρ of PETN CJ =2388kg / m 3 , then we get the density prediction result at position Q at time t
[0032] 2. Radial velocity prediction result u t (r)
[0033] First follow The normalized radial velocity at position Q at time t is calculated as
[0034] Then, combined with the CJ detonation parameters of PETN, the radial velocity prediction result u at position Q at time t is obtained: t (r), the CJ detonation parameters of PETN used in this step include the radial velocity u of the detonation products of PETN CJ =2100m / s, then the radial velocity prediction result at position Q at time t is obtained
[0035] 3. Pressure prediction results p t (r)
[0036] First follow Calculate the normalized pressure at position Q at time t
[0037] Then, combined with the CJ detonation parameters of PETN, the pressure prediction result p at position Q at time t is obtained. t (r), the CJ detonation parameters of PETN used in this step include the detonation pressure p of PETN CJ =30.42GPa, then the pressure prediction result of position Q at time t is obtained
[0038] 4. Prediction result of unit mass internal energy e t (r)
[0039] First follow Calculate the normalized unit mass internal energy at position Q at time t
[0040] Then, combined with the CJ detonation parameters of PETN, the prediction result of unit mass internal energy at position Q at time t is obtained: t(r), the CJ detonation parameters of PETN used in this step include the internal energy per unit mass e of the detonation products of PETN CJ = 10.38 J / kg, then the prediction result of the internal energy per unit mass at position Q at time t is obtained
[0041]
[0042] It should be noted that in the scenario of quantitative prediction of detonation state, the quantitative prediction of the state of the internal area of the detonation wave front is mainly concerned. Therefore, the calculation domain referred to in this application is mainly the area between the initiation point and the detonation wave front. So the distance r between any position Q in the calculation domain and the initiation point does not exceed the distance R0(t) between the detonation wave front and the initiation point. The dimensionless normalized spacing of any position in the corresponding calculation domain at any time However, in fact, this method can also be used to quantitatively predict the state parameters in the area outside the detonation wave front. Then for any position in the area of
[0043] This application provides a quantitative prediction function for the state parameters at any position Q in the calculation domain at time t after the central initiation of a three-dimensional spherical PETN charge. Using this quantitative prediction function and combining with the CJ detonation parameters of PETN, the prediction results of the state parameters at any position in the calculation domain at any time can be quickly and quantitatively calculated without introducing a complex detonation calculation model. Therefore, the calculation complexity and the amount of computation are both small, and the prediction speed is fast.
[0044] To verify the accuracy of the method of this application for calculating the instantaneous detonation wave of PETN charge, the prediction results of this application are respectively compared with the calculation results of AUTODYN software. This calculation example is designed as follows: Calculate the propagation of the detonation wave and the process of multiple reflections inside the solid wall after the central initiation of a PETN charge in a three-dimensional closed space of a cube. The side length of the cube is 20 cm. At a distance of 10 cm from the center of the cube, the time history curves of the flow field density obtained by the prediction results of this application method and the calculation results of AUTODYN software are compared as Figure 2 shown Figure 2 In , the black solid line is the prediction result obtained by using the method of this application, and the red dashed line is the calculation result of AUTODYN software.
[0045] It can be seen from the comparison results that the prediction results of the state parameters quantitatively calculated by the quantitative prediction function provided by this application are basically consistent within the error range with the prediction results of the state parameters obtained by using the detonation calculation model, indicating that the accuracy of quantitative prediction using the method of this application can reach the same level as that of the detonation calculation model. Therefore, the method of this application can be used for quantitative high-precision prediction in the far and near fields, and from Figure 2It can also be seen that the ability of the method of the present application to capture the peak pressure of the detonation wave is actually better than that of commercial software.
[0046] In a simulation example, a three-dimensional coordinate system XYZ is established with the initiation point of a three-dimensional spherical PETN charge as the coordinate origin. Due to symmetry, only the state quantitative prediction of the target prediction area composed of X ∈ [0m 3 , 0.16m 3 , Y ∈ [0m 3 , 0.16m 3 , Z ∈ [0m 3 , 0.16m 3 is analyzed. In this simulation example, a uniform grid is adopted, and dX = dY = 1mm. After a central detonation occurs at the coordinate origin of the three-dimensional spherical PETN charge, the distance R0(t) between the detonation wave front and the initiation point at t = 5 μs after central detonation is 40.25 mm, the distance R0(t) between the detonation wave front and the initiation point at t = 10 μs after central detonation is 80.5 mm, the distance R0(t) between the detonation wave front and the initiation point at t = 15 μs after central detonation is 120.75 mm, and the distance R0(t) between the detonation wave front and the initiation point at t = 20 μs after central detonation is 161 mm.
[0047] Using the method provided by the present application, the density prediction results ρ t (r) of the target prediction area at t = 5 μs, t = 10 μs, t = 15 μs, and t = 20 μs after central detonation can be quickly and quantitatively predicted as shown in Figure 3 (a), (b), (c), and (d) therein.
[0048] Using the method provided by the present application, the X-direction component U of the radial velocity prediction results u t (r) of the target prediction area at t = 5 μs, t = 10 μs, t = 15 μs, and t = 20 μs after central detonation can be quickly and quantitatively predicted as shown in Figure 4 (a), (b), (c), and (d) therein. The Y-direction component V of the radial velocity prediction results u t (r) at t = 5 μs, t = 10 μs, t = 15 μs, and t = 20 μs is as shown in Figure 5 (a), (b), (c), and (d) therein. The Z-direction component W of the radial velocity prediction results u t (r) at t = 5 μs, t = 10 μs, t = 15 μs, and t = 20 μs is as shown in Figure 6 (a), (b), (c), and (d) therein.
[0049] The method provided in this application can be used to quickly and quantitatively predict the pressure prediction results p of the target prediction area at t = 5μs, t = 10μs, t = 15μs, and t = 20μs after the central detonation. t (r) Figure 7 As shown in (a), (b), (c) and (d) in FIG.
[0050] The method provided in this application can be used to quickly and quantitatively predict the unit mass internal energy prediction results e of the target forecast area at t = 5μs, t = 10μs, t = 15μs, and t = 20μs after the central detonation. t (r) Figure 8 As shown in (a), (b), (c) and (d) in FIG.
[0051] 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 quantitatively predicting the instantaneous detonation state of a three-dimensional PETN charge, characterized in that, The state quantitative prediction method includes: Normalize the distance r between any position Q in the computational domain and the initiation point by the detonation wave front at any time t after the center initiation of a three-dimensional spherical PETN charge, to obtain the normalized spacing of position Q at time t According to calculate the normalized density at position Q at time t and combine with the CJ detonation parameters of PETN to obtain the density prediction result ρ t (r) at position Q at time t; According to calculate the normalized radial velocity at position Q at time t and combine with the CJ detonation parameters of PETN to obtain the radial velocity prediction result u t (r); According to Calculate the normalized pressure at position Q at time t And combine with the CJ detonation parameters of PETN to obtain the pressure prediction result p(r) at position Q at time t t (r); According to calculate the normalized internal energy per unit mass at position Q at time t and combine with the CJ detonation parameters of PETN to obtain the predicted result e of the internal energy per unit mass at position Q at time t t (r).
2. The state quantitative prediction method according to claim 1, characterized in that, Obtain the normalized spacing at position Q at time t Comprising: Based on the detonation velocity D = 8050 m / s of PETN, the distance R0(t) between the detonation wave front and the initiation point at any time t after the center initiation of the three-dimensional spherical PETN charge is determined as R0(t) = D × t, and the normalized spacing at position Q at time t is obtained through normalization processing 3. The state quantitative prediction method according to claim 1, characterized in that The CJ detonation parameters of PETN include the detonation product density ρ CJ = 2388 kg / m 3 , the radial velocity u of the detonation products CJ = 2100 m / s, the detonation pressure p CJ = 30.42 GPa, and the specific internal energy e of the detonation products CJ = 10.38 J / kg; Obtain the density prediction result ρ at position Q at time t t (r), the radial velocity prediction result u t (r), the pressure prediction result p t (r), and the internal energy prediction result e per unit mass t (r) includes calculating according to Calculate.