A rapid simulation and assessment method and system for large-scale shock wave damage

Abstract

This invention discloses a rapid simulation and evaluation method and system for large-scale shock wave damage. The method includes: S1, calculation of shock wave power parameters; S2, rapid acquisition of explosion shock wave power field parameters; S3, establishment of target data structure, and establishment of a target data description structure suitable for explosion effect analysis based on mirrored virtual sources; S4, damage analysis of the target by the explosion shock wave; S5, visualization of the target damage scene; S6, using a dynamic physics engine, integrating the engine with the explosion shock damage algorithm, combining graphical interface technology and graphical visualization technology to practically integrate the explosion shock damage algorithm, forming a rapid simulation analysis of explosion wave damage. This invention approximates the linearization of the propagation law of explosion waves, introduces virtual sources to handle the propagation, reflection, and diffraction problems of explosion waves in complex scenes, and forms a source effect superposition algorithm based on the basic conservation laws of fluid mechanics to calculate the nonlinear superposition effect of explosion waves from multiple virtual sources.

Description

Technical Field

[0001] This invention relates to the field of damage assessment technology, and more specifically, to a rapid simulation assessment method and system for large-scale shock wave damage. Background Technology

[0002] The shock wave generated by an explosion is a significant destructive element. Its propagation and effect on targets are complex hydrodynamic phenomena, and it has long been an important research topic in the field of damage performance assessment. Currently, many practical application problems related to explosion shock waves involve large-scale and complex scenarios, requiring rapid acquisition of explosion wave analysis and calculation results.

[0003] (1) Simulation and calculation of explosion wave propagation in complex environments and large-scale scenes. Complex environments include urban blocks /

[0004] In environments such as buildings and tunnels, blast waves exhibit complex reflection and diffraction phenomena, and the spatial scale involved is also large. When calculating and analyzing the effects of explosions, empirical formulas cannot be directly applied, and CFD analysis methods also suffer from problems such as excessive computation time and lack of flexibility.

[0005] (2) Simulation of Explosion Effects under Numerous Working Conditions. In engineering design, it is often necessary to perform explosion effect simulations under numerous different working conditions for the same scenario to identify vulnerable parts of the engineering structure and optimize it. Traditionally, CFD analysis methods are used to calculate these problems. However, due to the excessively long calculation time for a single working condition, a compromise can only be made by reducing the number of working conditions. Consequently, the number of simulated working conditions is very limited, and the lack of sample calculation results leads to limitations in the analytical conclusions.

[0006] (3) Rapid prediction of explosion effects. Since explosives can detonate at any time and time is of the essence, it is necessary to obtain calculation and analysis results within tens of seconds or minutes, quickly providing the damage range and effects of the blast wave, and giving an accurate assessment. Traditional CFD methods clearly cannot meet these requirements.

[0007] To address the challenge of rapid simulation and evaluation of large-scale explosion shockwave damage, it is necessary to develop a rapid simulation technology for shockwave damage. This technology should consider both fluid dynamics principles and employ approximate modeling methods to significantly improve the speed and efficiency of explosion shockwave damage simulation. Furthermore, combining graphical visualization and interactive technologies, a corresponding software system should be developed to support experimental simulation and evaluation. Therefore, it is essential to provide a rapid simulation and evaluation method and system for large-scale shockwave damage to solve the aforementioned technical problems. Summary of the Invention

[0008] The purpose of this invention is to provide a rapid simulation and evaluation method and system for large-scale shock wave damage, so as to overcome the defects of the existing technology.

[0009] To achieve the above objectives, the technical solution adopted by the present invention is as follows:

[0010] A rapid simulation evaluation method for large-scale shock wave damage includes the following steps:

[0011] S1. The simplified Kingery model is used to calculate the shock wave power parameters characterizing the damage effect, and a mathematical description model of the propagation process is established.

[0012] S2. Based on the mirror source method, the reflection and diffraction effects of the explosion shock wave are processed by multiple virtual mirror explosion sources. The explosion shock wave power parameters suitable for the virtual mirror source method are established, and the power field parameters are obtained by superposition algorithm.

[0013] S3. Establish the target data structure and, based on the mirrored virtual source, establish a target data description structure suitable for explosion effect analysis;

[0014] S4. Calculate the relative relationship between the surface element and the explosion center by establishing a virtual ray method, and calculate the explosion shock wave reflection pressure and sweep pressure;

[0015] S5. A spherical surface is used to represent the peak pressure front of the explosion shock wave. Color fusion technology is used to make the color of the spherical surface correspond to a fusion factor. As the explosion shock wave propagates, the color of the spherical surface becomes darker in the early stage to represent the increase of the shock wave pressure.

[0016] S6. A dynamic physics engine is adopted and integrated with the explosion impact damage algorithm. Graphical interface technology and graphical visualization technology are combined to practically integrate the explosion impact damage algorithm and form a rapid simulation analysis of explosion wave damage.

[0017] Further, step S1 specifically includes:

[0018] S11. Obtain the shock wave overpressure, positive pressure duration, and specific impulse of the explosive air shock wave;

[0019] S12. The shock wave power parameters characterizing the damage effect are calculated based on the shock wave overpressure, barometric pressure duration and specific impulse, using a simplified Kingery-Bulmash model.

[0020] S13. Integrate steps S11 and S12 into a calculation to form a dedicated program for calculating the power parameters of the explosion wave.

[0021] Furthermore, in step S11:

[0022] The formula for calculating shock wave overpressure is: Δp m =p m -p0;

[0023] The formula for calculating specific impulse is:

[0024] In the formula, τ + This represents the duration of positive pressure.

[0025] Further, step S12 specifically includes:

[0026] Assuming a known equivalent TNT quantity W, and setting the distance from the charge to R, first calculate the proportional distance Z:

[0027]

[0028] Then calculate the temporary variable T:

[0029] T = lnZ

[0030] In the simplified Kingery-Bulmash model, the power parameter F of a certain shock wave is expressed as a polynomial of T with the exponent e:

[0031] F = exp(A + BT + CT) 2 +DT 3 +ET 4 +FT 5 +GT 6 ).

[0032] Furthermore, step S13 specifically includes:

[0033] S131. Calculate the peak overpressure p at a certain distance using the pressure time history algorithm. op duration of positive pressure t d The specific impulse II in the barotropic region is described by the Friedlander equation, which describes the pressure-time history curve of the explosion shock wave:

[0034]

[0035] In the formula, p op It is the peak overpressure, t a It is the arrival time of the shock wave, t d It is the duration of positive pressure, and α is the waveform coefficient;

[0036] Determine the waveform coefficient α, perform time integration over pressure p(t), and from t a Start, integrate up to t a +t d The integral result is guaranteed to be equal to the specific impulse II in the positive pressure zone, that is:

[0037]

[0038] The nonlinear equation for the waveform coefficient α is obtained and solved by iterative method;

[0039] S132. Determine the positive reflection overpressure or free field overpressure using the formula of the reflection pressure algorithm:

[0040] p a =p op (1+cosα-2cos 2 α)+p r cos 2 α

[0041] In the formula, pop is the free-field peak overpressure, pr is the positive reflection overpressure, and α is the incident angle. If the incident angle is 0 degrees, it is perpendicular incidence, and the positive reflection overpressure is obtained; if the incident angle is 90 degrees, it is glancing incidence, and the free-field overpressure is obtained.

[0042] S133. Display and output the explosion wave pressure waveform, time history curve, and spatiotemporal cloud map distribution curve under a certain equivalent setting;

[0043] S134. Calculation of reflected pressure and impulse based on the ConWep model.

[0044] Further, step S2 specifically includes:

[0045] S21. Perform explosion wave simulation calculations to obtain the effective TNT charge, shock wave overpressure field distribution, specific impulse distribution, and power field parameters of the target component under overpressure and impulse conditions.

[0046] S22. Optimize the explosion wave reflection algorithm by equating the reflection of the explosion wave generated by the real explosion source to the direct action of the corresponding virtual explosion source. Then, use the superposition rule of multiple virtual source effects to comprehensively calculate the interaction of explosion waves from each virtual source and calculate the parameter changes of overpressure and impulse at the monitoring point under the mutual influence of different explosion waves.

[0047] Furthermore, step S3 specifically includes:

[0048] S31. Establish the terrain geometry data structure;

[0049] S32. Use model files and parameter files to describe the target data description structure.

[0050] Further, step S4 specifically includes:

[0051] S41. Calculate the relative relationship between the surface element and the explosion center using the virtual ray method;

[0052] S42. Calculate the shock wave overpressure value Δp at the surface element. R :

[0053] Overpressure value Δp R This is the reflected pressure after considering the incident angle θ of the shock wave, and the overpressure value Δp. I For the shock wave to pass over the pressure, Δp R and Δp I Calculated by the following formula

[0054] Δp R =Δp2·cos 2 θ+Δp1·(1+cos 2 θ-2cosθ)

[0055]

[0056] In the formula, Δp2 and Δp1 are the shock wave positive reflection overpressure and free field overpressure at the point, respectively, p0 is the atmospheric pressure, and Δp1 is obtained according to the shock wave peak attenuation formula.

[0057] Furthermore, step S5 also includes:

[0058] The overpressure magnitude on the target is calculated using an overpressure field algorithm to determine the overpressure distribution, which is then displayed as a contour map. The contour map uses color changes to represent the pressure magnitude.

[0059] The 1D texture mapping technology provided by the graphics engine uses color bars as 1D textures to interpolate the color representing the pressure magnitude into the color display of the target surface element.

[0060] The present invention also provides a system for implementing the above-mentioned rapid simulation and evaluation method for large-scale shock wave damage, comprising:

[0061] The mathematical description unit for the power parameters of the explosion shock wave is used to calculate the power parameters of the shock wave that characterize the damage effect using a simplified Kingery model, and to establish a mathematical description model of the propagation process.

[0062] The explosive shock wave power field parameter rapid acquisition unit is used to process the reflection and diffraction effects of explosive shock waves using multiple virtual mirror explosion sources based on the mirror source method, establish explosive shock wave power parameters suitable for the virtual mirror source method, and obtain the power field parameters using a superposition algorithm;

[0063] The target data structure establishment unit is used to establish the target data structure and establish a target data description structure suitable for explosion effect analysis based on the mirror virtual source.

[0064] The damage analysis unit for the target caused by the explosion shock wave is used to calculate the relative relationship between the surface element and the explosion center by establishing a virtual ray, and to calculate the explosion shock wave reflection pressure and sweep pressure.

[0065] The target damage scene visualization unit is used to represent the peak pressure front of the explosion shock wave using a spherical surface. Color fusion technology is used to make the color of the spherical surface correspond to a fusion factor. As the explosion shock wave propagates, the color of the spherical surface initially becomes darker, representing that the shock wave pressure is increasing.

[0066] The target response simulation technology and system integration unit is used to integrate the dynamic physics engine with the explosion impact damage algorithm, and combine graphical interface technology and graphical visualization technology to practically integrate the explosion impact damage algorithm to form a rapid simulation analysis of explosion wave damage.

[0067] Compared with the prior art, the advantages of the present invention are as follows: The present invention approximates the propagation law of the explosion wave as linear, introduces virtual sources to handle the propagation, reflection and diffraction problems of the explosion wave in complex scenes, and forms a source effect superposition algorithm based on the basic conservation law of fluid mechanics to calculate the nonlinear superposition effect of multiple virtual source explosion waves, and establishes a rapid simulation evaluation method for shock wave damage in large scenes. Attached Figure Description

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

[0069] Figure 1 This is a diagram showing the relationship between the large-scale shock wave damage rapid simulation evaluation method of this invention and existing analysis methods.

[0070] Figure 2 This is a flowchart of the method for mathematical description and parameter calculation model of the explosive shock wave power parameters in the rapid damage simulation assessment method provided in this embodiment of the invention.

[0071] Figure 3 In the diagrams a and b, respectively, are the timing and point waveforms of the shock wave in the mathematical description and parameter calculation model of the explosive shock wave power parameters during the rapid damage simulation assessment method provided in this embodiment of the invention.

[0072] Figure 4 This is a diagram showing the calculation of reflection pressure and impulse based on the ConWep model during the rapid damage simulation assessment method provided in this embodiment of the invention.

[0073] Figure 5 This is a structural diagram showing the rapid acquisition of explosive shock wave force field parameters during the rapid damage simulation and evaluation method provided in this embodiment of the invention.

[0074] Figure 6 This is a flowchart illustrating the damage analysis of a target by an explosion shock wave during the rapid damage simulation and evaluation method provided in this embodiment of the invention.

[0075] Figure 7 This is a visualization example of the target damage scene in the rapid damage simulation and evaluation method provided in this embodiment of the invention.

[0076] Figure 8 These are stress cloud color bars used to visualize the target damage scene during the rapid damage simulation and evaluation method provided in this embodiment of the invention.

[0077] Figure 9 This is a diagram illustrating the practical application of the target response simulation and system integration in the rapid damage simulation evaluation method provided in this embodiment of the invention.

[0078] Figure 10 This is a structural diagram of the large-scale shock wave damage rapid simulation and evaluation system provided in this embodiment of the invention. Detailed Implementation

[0079] The preferred embodiments of the present invention will now be described in detail with reference to the accompanying drawings, so that the advantages and features of the present invention can be more easily understood by those skilled in the art, thereby providing a clearer and more explicit definition of the scope of protection of the present invention.

[0080] See Figure 1 As shown in the figure, this embodiment discloses a rapid simulation evaluation method for large-scale shock wave damage, including the following steps:

[0081] Step S1: Mathematical description and calculation model of the explosive shock wave power parameters. A simplified Kingery model is used to calculate the shock wave power parameters characterizing the damage effect, and a mathematical description model of the propagation process is established.

[0082] Step S2: Rapid acquisition of explosive shock wave power field parameters. Based on the mirror source method, the reflection and diffraction effects of the explosive shock wave are processed using multiple virtual mirror sources. Explosive shock wave power parameters suitable for the virtual mirror source method are established. Based on the fundamental laws of mass conservation, momentum conservation, and energy conservation, a superposition algorithm is used to obtain the power field parameters, providing a basis for judging the degree of damage for the damage assessment model.

[0083] Step S3: Establish a large-scale spatial target data structure, and based on the consideration of mirror virtual source, establish a target data description structure suitable for explosion effect analysis.

[0084] Step S4: Damage analysis of the target caused by the explosion shock wave. The relative relationship between the surface element and the explosion center is calculated by establishing a virtual ray, and the reflected pressure and passing pressure of the explosion shock wave are calculated by empirical formulas (using existing algorithms).

[0085] Step S5: Visualization of target damage scene. A sphere is used to represent the peak pressure front of the blast shock wave. Color fusion technology is used to make the color of the sphere correspond to a fusion factor. As the blast shock wave propagates, the color of the sphere initially becomes darker, representing the increase in shock wave pressure. That is, as the blast shock wave propagates, the fusion factor changes, and the degree of target damage is presented through the color of the sphere.

[0086] Step S6: Target Response Simulation and System Integration. A dynamic physics engine is employed and integrated with an explosion impact damage algorithm. Graphical interface and visualization technologies are combined to practically integrate the explosion impact damage algorithm, forming a rapid simulation analysis of explosion wave damage.

[0087] See Figure 2 As shown, step S1 in this embodiment specifically includes:

[0088] Step S11: Based on the generation mechanism and calculation formula of the explosion shock wave, obtain the shock wave overpressure, positive pressure duration and specific impulse of the explosion air shock wave.

[0089] Step S12: Calculate the shock wave power parameters characterizing the damage effect based on the shock wave overpressure, barometric pressure duration, and specific impulse using a simplified Kingery-Bulmash model.

[0090] Step S13: Integrate steps S11 and S12 into a calculation to form a dedicated program for calculating the explosive wave power parameters.

[0091] Preferably, in step S11, the explosive detonates in the air, transforming into high-temperature and high-pressure detonation products in a very short time. Because the initial pressure and density of the air are low, the detonation products expand rapidly, leading to a decrease in pressure and density, forming rarefaction waves within the detonation products. Simultaneously, the expansion of the detonation products strongly compresses the surrounding air, forming air shock waves.

[0092] The change of shock wave pressure over time is as follows Figure 3 As shown in a and b, the formula for calculating shock wave overpressure (the magnitude of pressure exceeding atmospheric pressure) is: Δp m =p m -p0;

[0093] Specific impulse I + The calculation formula is: In the formula, τ + This represents the duration of positive pressure.

[0094] Preferably, step S12 specifically includes:

[0095] Assuming a known equivalent TNT weight W (in kg), and setting the distance from the charge to the charge to R (in m), first calculate the proportional distance Z:

[0096]

[0097] Then calculate the temporary variable T (by taking the natural logarithm of Z):

[0098] T = lnZ

[0099] In the simplified Kingery-Bulmash model, the power parameter F of a certain shock wave is expressed as a polynomial of T with the exponent e:

[0100] F = exp(A + BT + CT) 2 +DT 3 +ET 4 +FT 5 +GT 6 ).

[0101] Except for the peak overpressure and shock wave velocity, other parameters need to be proportionalized to the equivalent. The coefficients of the T polynomial also differ for different power parameters, as shown in Tables 1 to 5.

[0102] Calculate the arrival time t of the shock wave a The parameters (in milliseconds) are shown in Table 1 (5).

[0103] Table 1. Polynomial coefficients for calculating the arrival time of the shock wave.

[0104]

[0105] at this time:

[0106]

[0107] Calculate the peak overpressure p of the shock wave so The parameters (in kPa) are shown in Table 2.

[0108] Table 2 shows the polynomial coefficients for calculating the peak overpressure.

[0109]

[0110] at this time:

[0111] p so =exp(A+BT+CT) 2 +DT 3 +ET 4+FT 5 +GT 6 )

[0112] Calculate the duration τ of the shock wave barypressure + The parameters (in milliseconds) are shown in Table 3 (7).

[0113] Table 3 shows the polynomial coefficients for calculating the duration of positive pressure.

[0114]

[0115] at this time:

[0116]

[0117] Calculate the positive pressure impulse I + The parameters (in kPa·ms) are shown in Table 4 (8).

[0118] Table 4 Polynomial coefficients for calculating the positive pressure impulse

[0119]

[0120] at this time:

[0121]

[0122] The parameters for calculating the shock wave velocity D (in km / s) are shown in Table 5.

[0123] Table 5 Polynomial coefficients for calculating shock wave velocity

[0124]

[0125] Preferably, step S13 specifically includes:

[0126] Step S131: Calculate the peak overpressure p at a certain distance using the pressure time history algorithm. op duration of positive pressure t d The specific impulse II in the barotropic region is described by the Friedlander equation, which describes the pressure-time history curve of the explosion shock wave:

[0127]

[0128] In the formula, p op It is the peak overpressure, t a It is the arrival time of the shock wave, t d α is the duration of positive pressure, and α is the waveform coefficient.

[0129] Determine the waveform coefficient α, perform time integration over pressure p(t), and from t a Start, integrate up to t a +td The integral result is guaranteed to be equal to the specific impulse II in the positive pressure zone, that is:

[0130]

[0131] The nonlinear equation for the waveform coefficient α is obtained and solved by iterative method.

[0132] Step S132: Determine the positive reflection overpressure or free field overpressure using the formula of the reflection pressure algorithm (i.e., the ConWep reflection model):

[0133] p a =p op (1+cosα-2cos 2 α)+p r cos 2 α

[0134] In the formula, pop is the peak overpressure of the free field, pr is the positive reflection overpressure, and α is the incident angle. If the incident angle is 0 degrees, it is a perpendicular incident, and the positive reflection overpressure is obtained; if the incident angle is 90 degrees, it is a glancing incident, and the free field overpressure is obtained.

[0135] Step S133: Pressure waveform and time history curve. Considering the influence of altitude, the output shows the explosion wave pressure waveform, time history curve, and spatiotemporal cloud map distribution curve under a certain equivalent setting.

[0136] Step S134: Calculation of reflected pressure and impulse based on the ConWep model, as follows... Figure 4 As shown.

[0137] See Figure 5 As shown, step S2 in this embodiment specifically includes:

[0138] Step S21: Perform relevant calculations for explosion wave simulation to obtain power field parameters such as effective TNT charge, shock wave overpressure field distribution, specific impulse distribution, and the overpressure and impulse state of the target component.

[0139] Step S22: Optimize the explosion wave reflection algorithm by equating the reflection of the explosion wave generated by the real explosion source to the direct action of the corresponding virtual explosion source. Then, use the superposition rule of multiple virtual source effects to comprehensively calculate the interaction of explosion waves from each virtual source and calculate the changes in parameters such as overpressure and impulse at the monitoring point under the mutual influence of different explosion waves.

[0140] Step S3 in this embodiment specifically includes:

[0141] Step S31: Establish the terrain geometry (elevation) data structure.

[0142] Step S32: The target model mainly consists of buildings, and the target data description structure is described using model files and parameter files.

[0143] See Figure 6 As shown, step S4 in this embodiment specifically includes:

[0144] Step S41: Calculate the relative relationship between the surface element and the explosion center using the virtual ray method.

[0145] Step S42: Calculate the shock wave overpressure value Δp at the surface element. R .

[0146] Overpressure value Δp R This is the reflected pressure after considering the incident angle θ of the shock wave, and the overpressure value Δp. I For the shock wave to pass over the pressure, Δp R and Δp I It is calculated using the following formula:

[0147] Δp R =Δp2·cos 2 θ+Δp1·(1+cos 2 θ-2cosθ)

[0148]

[0149] In the formula, Δp2 and Δp1 are the shock wave positive reflection overpressure and free field overpressure at the point, respectively, p0 is the atmospheric pressure, and Δp1 is obtained according to the shock wave peak attenuation formula. Note that this algorithm implicitly considers the target surface as a rigid wall, that is, the shock wave is reflected on a rigid wall, which is acceptable for calculating the initial load on most metal surfaces and reinforced concrete surfaces.

[0150] See Figure 6 As shown, step S5 of this embodiment further includes:

[0151] Step S51: Color fusion technology. Considering that the peak pressure on the array surface gradually decreases, the spherical color corresponds to a fusion factor. As the explosion shock wave propagates, the fusion factor changes, making the initial color of the spherical surface more intense, representing higher shock wave pressure. As time progresses, the spherical surface color gradually fades and becomes translucent, representing the gradual decrease in pressure.

[0152] Step S52: Pressure Distribution Cloud Map. The overpressure magnitude on the target is calculated using the overpressure field algorithm described above, and the calculated overpressure distribution is displayed as a cloud map. The cloud map uses color changes to represent pressure magnitude, so it is crucial to define a color bar that represents pressure magnitude. The defined color bar is as follows: Figure 8 As shown.

[0153] Step S53: The 1D texture mapping technology provided by the graphics engine uses color bars as 1D textures to interpolate the colors representing pressure levels into the color display of the target surface elements. This not only allows different pressure values ​​to be displayed in different colors, but also ensures that each displayed color corresponds to a pressure value. Figure 7 The image shows the distribution of the overpressure field on the target at different times. It can be seen that the overpressure on the nose of the aircraft is greater, while the overpressure on the tail is smaller.

[0154] See Figure 9 As shown, this embodiment demonstrates a rapid simulation analysis of the damage caused by the explosion wave generated in step S6. Figure 8 This provides a calculation example of the dynamic response of ground vehicles, including the initial scenario and simulation results.

[0155] See Figure 10 As shown, the present invention also provides a system for implementing the above-mentioned rapid simulation and evaluation method for large-scale shock wave damage, comprising:

[0156] The mathematical description unit 1 for the explosive shock wave power parameters is used to calculate the shock wave power parameters characterizing the damage effect using a simplified Kingery model, and to establish a mathematical description model of the propagation process.

[0157] Unit 2 for rapid acquisition of explosive shock wave power field parameters is used to process the reflection and diffraction effects of explosive shock waves using multiple virtual mirror explosion sources based on the mirror source method, establish explosive shock wave power parameters suitable for the virtual mirror source method, and obtain the power field parameters using a superposition algorithm.

[0158] Target data structure establishment unit 3 is used to establish the target data structure and establish a target data description structure suitable for explosion effect analysis based on the mirror virtual source.

[0159] The damage analysis unit 4 of the explosion shock wave to the target is used to calculate the relative relationship between the surface element and the explosion center by establishing a virtual ray method, and to calculate the explosion shock wave reflection pressure and sweep pressure.

[0160] The target damage scene visualization unit 5 is used to represent the peak pressure front of the explosion shock wave using a spherical surface. Color fusion technology is used to make the color of the spherical surface correspond to a fusion factor. As the explosion shock wave propagates, the color of the spherical surface initially becomes darker, representing that the shock wave pressure is increasing.

[0161] The target response simulation technology and system integration unit 6 is used to integrate the dynamic physics engine with the explosion impact damage algorithm, and combine graphical interface technology and graphical visualization technology to practically integrate the explosion impact damage algorithm to form a rapid simulation analysis of explosion wave damage.

[0162] This invention approximates the linearization of the propagation law of explosion waves, introduces virtual sources to handle the propagation, reflection and diffraction problems of explosion waves in complex scenes, and forms a source effect superposition algorithm based on the basic conservation law of fluid mechanics to calculate the nonlinear superposition effect of explosion waves from multiple virtual sources, thus establishing a rapid simulation and evaluation method for shock wave damage in large scenes.

[0163] Although embodiments of the present invention have been described in conjunction with the accompanying drawings, the patent owner may make various modifications or alterations within the scope of the appended claims, as long as they do not exceed the protection scope described in the claims of the present invention, they shall be within the protection scope of the present invention.

Claims

1. A rapid simulation and evaluation method for large-scale shock wave damage, characterized in that: Includes the following steps: S1. The simplified Kingery model is used to calculate the shock wave power parameters characterizing the damage effect, and a mathematical description model of the propagation process is established. S2. Based on the mirror source method, the reflection and diffraction effects of the explosion shock wave are processed by multiple virtual mirror explosion sources. The explosion shock wave power parameters suitable for the virtual mirror source method are established, and the power field parameters are obtained by superposition algorithm. S3. Establish the target data structure and, based on the mirrored virtual source, establish a target data description structure suitable for explosion effect analysis; S4. Calculate the relative relationship between the surface element and the explosion center by establishing a virtual ray method, and calculate the explosion shock wave reflection pressure and sweep pressure; S5. A spherical surface is used to represent the peak pressure front of the explosion shock wave. Color fusion technology is used to make the color of the spherical surface correspond to a fusion factor. As the explosion shock wave propagates, the color of the spherical surface becomes darker in the early stage to represent the increase of the shock wave pressure. S6. A dynamic physics engine is adopted and integrated with the explosion impact damage algorithm. Graphical interface technology and graphical visualization technology are combined to practically integrate the explosion impact damage algorithm and form a rapid simulation analysis of explosion wave damage.

2. The rapid simulation and evaluation method for large-scale shock wave damage according to claim 1, characterized in that, Step S1 specifically includes: S11. Obtain the shock wave overpressure, positive pressure duration, and specific impulse of the explosive air shock wave; S12. The shock wave power parameters characterizing the damage effect are calculated based on the shock wave overpressure, barometric pressure duration and specific impulse, using a simplified Kingery-Bulmash model. S13. Integrate steps S11 and S12 into a calculation to form a dedicated program for calculating the power parameters of the explosion wave.

3. The rapid simulation and evaluation method for large-scale shock wave damage according to claim 2, characterized in that, In step S11: The formula for calculating shock wave overpressure is: Δp m =p m -p0; The formula for calculating specific impulse is: In the formula, τ + This represents the duration of positive pressure.

4. The rapid simulation and evaluation method for large-scale shock wave damage according to claim 2, characterized in that, Step S12 specifically includes: Assuming a known equivalent TNT quantity W, and setting the distance from the charge to R, first calculate the proportional distance Z: Then calculate the temporary variable T: T = lnZ In the simplified Kingery-Bulmash model, the power parameter F of a certain shock wave is expressed as a polynomial of T with the exponent e: F=exp(A+BT+CT 2 +DT 3 +ET 4 +FT 5 +GT 6 )。 5. The rapid simulation and evaluation method for large-scale shock wave damage according to claim 2, characterized in that, Step S13 specifically includes: S131. Calculate the peak overpressure p at a certain distance using the pressure time history algorithm. op duration of positive pressure t d The specific impulse II in the barotropic region is described by the Friedlander equation, which describes the pressure-time history curve of the explosion shock wave: In the formula, p op It is the peak overpressure, t a It is the arrival time of the shock wave, t d It is the duration of positive pressure, and α is the waveform coefficient; Determine the waveform coefficient α, perform time integration over pressure p(t), and from t a Start, integrate up to t a +t d The integral result is guaranteed to be equal to the specific impulse II in the positive pressure zone, that is: The nonlinear equation for the waveform coefficient α is obtained and solved by iterative method; S132. Determine the positive reflection overpressure or free field overpressure using the formula of the reflection pressure algorithm: p a =p op (1+cosα-2cos 2 a)+p r cos 2 a In the formula, p op α is the peak overpressure, pr is the positive reflection overpressure, and α is the incident angle. If the incident angle is 0 degrees, it is perpendicular incidence, and the result is the positive reflection overpressure; if the incident angle is 90 degrees, it is glancing incidence, and the result is the free field overpressure. S133. Display and output the explosion wave pressure waveform, time history curve, and spatiotemporal cloud map distribution curve under a certain equivalent setting; S134. Calculation of reflection pressure and impulse based on the ConWep model.

6. The rapid simulation and evaluation method for large-scale shock wave damage according to claim 1, characterized in that, Step S2 specifically includes: S21. Perform explosion wave simulation calculations to obtain the effective TNT charge, shock wave overpressure field distribution, specific impulse distribution, and power field parameters of the target component under overpressure and impulse conditions. S22. Optimize the explosion wave reflection algorithm by equating the reflection of the explosion wave generated by the real explosion source to the direct action of the corresponding virtual explosion source. Then, use the superposition rule of multiple virtual source effects to comprehensively calculate the interaction of explosion waves from each virtual source and calculate the parameter changes of overpressure and impulse at the monitoring point under the mutual influence of different explosion waves.

7. The rapid simulation and evaluation method for large-scale shock wave damage according to claim 1, characterized in that, Step S3 specifically includes: S31. Establish the terrain geometry data structure; S32. Use model files and parameter files to describe the target data description structure.

8. The rapid simulation and evaluation method for large-scale shock wave damage according to claim 1, characterized in that, Step S4 specifically includes: S41. Calculate the relative relationship between the surface element and the explosion center using the virtual ray method; S42. Calculate the shock wave overpressure value Δp at the surface element. R : Overpressure value Δp R This is the reflected pressure after considering the incident angle θ of the shock wave, and the overpressure value Δp. I For the shock wave to pass over the pressure, Δp R and Δp I Calculated by the following formula Δp R =Δp2·cos 2 θ+Δp1·(1+cos 2 θ-2cosθ) In the formula, Δp2 and Δp1 are the shock wave positive reflection overpressure and free field overpressure at the surface element, respectively, p0 is the atmospheric pressure, and Δp1 is obtained according to the shock wave peak attenuation formula.

9. The rapid simulation and evaluation method for large-scale shock wave damage according to claim 1, characterized in that, Step S5 also includes: The overpressure magnitude on the target is calculated using an overpressure field algorithm to determine the overpressure distribution, which is then displayed as a contour map. The contour map uses color changes to represent the pressure magnitude. The 1D texture mapping technology provided by the graphics engine uses color bars as 1D textures to interpolate the color representing the pressure magnitude into the color display of the target surface element.

10. A system for implementing the rapid simulation and evaluation method for large-scale shock wave damage as described in any one of claims 1-9, characterized in that, include: The mathematical description unit for the power parameters of the explosion shock wave is used to calculate the power parameters of the shock wave that characterize the damage effect using a simplified Kingery model, and to establish a mathematical description model of the propagation process. The explosive shock wave power field parameter rapid acquisition unit is used to process the reflection and diffraction effects of explosive shock waves using multiple virtual mirror explosion sources based on the mirror source method, establish explosive shock wave power parameters suitable for the virtual mirror source method, and obtain the power field parameters using a superposition algorithm; The target data structure establishment unit is used to establish the target data structure and establish a target data description structure suitable for explosion effect analysis based on the mirror virtual source. The damage analysis unit for the target caused by the explosion shock wave is used to calculate the relative relationship between the surface element and the explosion center by establishing a virtual ray, and to calculate the explosion shock wave reflection pressure and sweep pressure. The target damage scene visualization unit is used to represent the peak pressure front of the explosion shock wave using a spherical surface. Color fusion technology is used to make the color of the spherical surface correspond to a fusion factor. As the explosion shock wave propagates, the color of the spherical surface initially becomes darker, representing that the shock wave pressure is increasing. The target response simulation technology and system integration unit is used to integrate the dynamic physics engine with the explosion impact damage algorithm, and combine graphical interface technology and graphical visualization technology to practically integrate the explosion impact damage algorithm to form a rapid simulation analysis of explosion wave damage.

Images