Damage efficiency simulation and evaluation system based on model simulation

Abstract

The invention provides a damage efficiency simulation and evaluation system based on model simulation, and relates to the technical field of explosive action simulation analysis. Three-dimensional scanning data of a target simulation explosion object structure is obtained, and a physical field simulation model of a target object is reconstructed based on a finite element grid; and the stress behavior and the stress distribution state of the structure in the actual blasting environment can be precisely restored. Explosive particle size distribution parameters, terrain geometric factors, shock wave reflection and energy dissipation physical parameters, a first correction module and a second correction module are introduced, an explosive pressure field and equivalent stress are corrected, the limitation that a traditional simulation model ignores environmental complexity is broken through, and evaluation distortion or misjudgment is avoided. In the simulation process, a stress-concentrated fragile area in the structure is automatically identified, and a damage grade distribution diagram is dynamically calibrated in combination with a corrected pressure field and stress evaluation, so that a hierarchical early warning instruction is generated, and a scientific basis is provided for engineering protection design, emergency response and risk control.

Description

Technical Field

[0001] The present invention relates to the technical field of explosive action simulation analysis, and particularly to a damage effectiveness simulation and evaluation system based on model simulation. Background Art

[0002] In scenarios such as mine exploitation, tunnel excavation, foundation blasting, and military attack and defense, blasting operations are often carried out under complex environmental conditions, such as drastic terrain undulations, compact target structure layouts, the presence of reflective obstacles or non-uniform media, etc. These environmental factors will significantly affect the propagation mode of explosion energy, and thus change the stress state of the target structure and its damage mode.

[0003] Existing blasting simulation and evaluation technologies usually rely on uniform simplified models, ignoring the complexity of the environmental terrain and its effects on phenomena such as the propagation path of shock waves, the energy attenuation rate, the reflection and superposition effect, and the formation of local overpressure areas. For example, mountains, walls, or underground structures may all trigger secondary reflections and shock wave focusing, resulting in the structural stress in local areas being much higher than the theoretical estimate, leading to "unexpected damage"; In summary, the existing technologies have problems such as ignoring the reflection, diffraction, and focusing of shock waves after encountering obstacles (such as walls, rock masses, etc.), and being unable to predict the local high-pressure damage area formed due to energy aggregation, resulting in insufficient accuracy.

[0004] The above information disclosed in the background art section is only used to enhance the understanding of the background of the present disclosure, and thus it may include information that does not constitute the prior art known to those of ordinary skill in the art. Summary of the Invention

[0005] The purpose of the present invention is to provide a damage effectiveness simulation and evaluation system based on model simulation to solve the problems raised in the above background art.

[0006] To achieve the above purpose, the present invention provides the following technical solutions: A damage effectiveness simulation and evaluation system based on model simulation, comprising: A target structure analysis module, configured to obtain three-dimensional scan data of the physical structure of the target simulated explosion, and based on the finite element mesh, reconstruct the physical field simulation model of the target object, establish a simulation domain and divide it into several shell elements. During the process of simulating gunpowder explosion, construct the internal pressure field of the explosive in the i-th shell element and the equivalent stress , and preset a pressure threshold Y, and compare the internal pressure field of the explosive in the i-th shell element with the pressure threshold Y to obtain a first evaluation result. According to the first evaluation result, construct an "initial damage level distribution map"; The first correction module is used to collect the particle size distribution parameters of explosive particles, obtain the explosive energy data set, and construct the energy release rate of the i-th shell unit The internal pressure field of the explosive for the i-th shell element Make corrections to obtain the corrected internal pressure field of the explosive for the i-th shell element ; and correct the "initial distribution map of damage level" to obtain the "corrected distribution map of damage level"; The second correction module is used to collect the terrain geometry information of the i-th shell unit, construct the environmental response data set, and construct the environmental impact index of the i-th shell unit. , and the equivalent stress of the i-th shell element Make corrections to obtain the corrected equivalent stress of the i-th shell element And evaluate, further identify the brittle fracture point area in the physical field simulation model, and dynamically link the brittle fracture point response with the damage level assessment standard and the corresponding early warning instructions during the damage simulation process.

[0007] Furthermore, the target structure analysis module includes a structure data acquisition unit and a finite element network construction unit: The structural data acquisition unit is used to use 3D-LiDAR or industrial CT scanning equipment to collect three-dimensional scanning data of the target simulated explosion physical structure, and after denoising and unified format conversion, a finite element network construction unit is used to establish a physical field simulation model based on the three-dimensional scanning data using Gmsh or AnsysMesh tools. After importing the three-dimensional scanning data, an initial three-dimensional grid is generated, and a simulation domain is established in the initial three-dimensional network, and the simulation domain includes: a target structure domain, an explosion source domain, an environmental medium domain, and a reflection boundary domain; and the simulation domain is divided into a plurality of shell units, which are marked in the physical field simulation model as: , to Represents the first shell element to the The status mark of each shell element at the time point sd.

[0008] Furthermore, the target structure analysis module also includes an explosion triggering process unit, a blasting pressure field construction unit and a blasting intensity division unit: The explosion triggering process unit is used to preset initial conditions, set the explosive heat value, define the detonation point and detonation speed in the simulation domain, and the detonation speed is set to 6930m / s; The blasting pressure field construction unit is used to construct the internal pressure field of the explosive of the i-th shell unit by applying the JWL state equation. : ; In the formula, and They respectively represent the pressure empirical constants in the initial stage and the late stage of the explosion, which are obtained through the copper cylinder expansion test. is the first attenuation coefficient, which controls the sharp attenuation in the initial stage of the explosion. is the second attenuation coefficient, which controls the slow attenuation in the late stage of the explosion. For TNT nitroglycerin explosive, ; For Comp-B explosive, ; For RDX explosive, ; The meaning of the formula is: the internal pressure field of the explosive in the i-th shell element is applied as an external load to the surface of the shell element, and stress wave propagation is generated at each time step to form the internal force response of the structure; the internal pressure field of the explosive in the i-th shell element describes the instantaneous pressure change generated by the explosion of the explosive, reflecting the influence of gas expansion and shock wave. represents the initial volume of the i-th shell element, represents the current volume of the i-th shell element, represents the energy ratio of the explosive, represents the internal energy per unit mass of the i-th shell element, which is calculated by the following formula: ; Among them, Q is the explosion value of all the explosives in the i-th shell element, with the unit of J / kg, represents the mass of the explosive in the i-th shell element, represents the explosive density; The blasting strength division unit is used to preset the pressure threshold Y, and compare the internal pressure field of the explosive in the i-th shell element with the pressure threshold Y to obtain the first evaluation result, including: When the internal pressure field of the explosive in the i-th shell element ≥ pressure threshold Y × 120%, it is marked as the high-pressure area; When pressure threshold Y ≤ the internal pressure field of the explosive in the i-th shell element < pressure threshold Y × 120%, it is marked as the medium-pressure area; When the internal pressure field of the explosive in the i-th shell element < pressure threshold Y, it is marked as the low-pressure area; Traverse all shell elements to obtain their internal pressure fields of the explosive and construct the "initial damage level distribution map" according to the first evaluation result.

[0009] Furthermore, the target structure analysis module further includes an equivalent stress construction unit: The equivalent stress construction unit is used to construct the equivalent stress of the i-th shell element under the action of an explosion load. : ; Wherein, is the equivalent plastic strain value of the i-th shell element, represents the equivalent strain rate of the i-th shell element, represents the reference equivalent strain rate, which is set to , is the normalized temperature value of the i-th shell element, defined as , T represents the current temperature, T 0 represents the reference temperature, T m represents the melting point temperature of the material, represents the strain hardening index, represents the thermal softening index, represents the initial yield strength of the material, represents the strain hardening coefficient, represents the strain rate sensitivity coefficient; it is obtained by experimental fitting, specifically: When the material is AISI4340 high-strength steel, K 1 = 792 MPa; K 2 = 510 MPa; K 3 = 0.014; n 1 = 0.26; m 1 = 1.03; T m = 1793 K, T 0 = 300 K When the material is 6061-T6 aluminum alloy, K 1 = 324 MPa; K 2 = 114 MPa; K 3 = 0.002; n 1 = 0.42; m 1 = 1.3403; T m = 7751793 K, T 0 = 293 K When the material is Ti-6Al-4V titanium alloy, K 1 = 1098 MPa; K 2 = 1092 MPa; K 3 = 0.014; n 1 = 0.93; m 1 = 1.1; T m = 1878 K, T 0 = 298 K.

[0010] Furthermore, the first correction module includes a first acquisition unit, a first correction unit, and a pressure evaluation unit; The first acquisition unit is configured to acquire the particle size distribution parameters of explosive particles to obtain an explosive energy data set, where the explosive energy data set includes: the energy release characteristics under different particle size combinations. The value range of the particle size d is divided into Z intervals, and the z-th interval represents explosive particles with an average particle size of ; For the i-th shell element, the proportion of each particle size interval is sampled as , combined with the explosive mass of the i-th shell element, to construct the energy release rate of the i-th shell element through the following formula : ; In the formula, represents the mass proportion of the explosive of explosive particles with an average particle size of , satisfying ; represents the release rate of the explosive with a particle size of , in units of J / kg·s, and is calculated through the following formula: ; Among them, q represents the total energy release amount of unit mass of explosive. For TNT nitroglycerin explosive, q = 4.2×10 6 J / kg; for Comp-B explosive, q = 5.5×10 6 J / k; for RDX explosive, q = 5.7×10 6 J / k; represents the energy release time of explosive particles with a particle size of ; The first correction unit is configured to correct the energy release rate of the i-th shell element for the internal pressure field of the explosive in the i-th shell element to obtain the corrected internal pressure field of the explosive in the i-th shell element: ; In the formula, represents the reference energy release rate under standard particle size conditions for normalization processing.

[0011] Furthermore, the pressure evaluation unit is configured to re-compare and correct the corrected internal pressure field of the explosive in the i-th shell element with the pressure threshold Y to obtain a corrected first evaluation result, including: When the corrected internal pressure field ≥ the pressure threshold Y × 120%, marked as the high-pressure area; When the pressure threshold Y ≤ the corrected explosive internal pressure field of the i-th shell element < the pressure threshold Y × 120%, marked as the medium-pressure area; When the corrected explosive internal pressure field of the i-th shell element < the pressure threshold Y, marked as the low-pressure area; Traverse all shell elements to obtain the corrected explosive internal pressure field within the explosives , and construct a "damage level correction distribution map" by comparing with the corrected first evaluation result.

[0012] Furthermore, the second correction module includes a second acquisition unit and an associated unit; The second acquisition unit is used to collect the topographic geometric information of the i-th shell element and construct an environmental response data set, which includes the shock wave attenuation factor of the i-th shell element , reflectivity factor and energy dissipation factor ; The shock wave attenuation factor of the i-th shell element , reflectivity factor and energy dissipation factor are obtained specifically as follows: S21. During the propagation of the shock wave, it will attenuate due to the influence of the terrain. As the distance from the explosion center to the target position increases, the shock wave intensity gradually weakens. The shock wave attenuation of the i-th shell element is calculated through the following formula : ; Among them, is the distance from the i-th shell element to the explosion source, represents the attenuation exponent, and its value range is 1 - 2; S22. When the shock wave encounters different types of ground, a certain proportion of the wave energy will be reflected; the reflectivity is affected by the terrain surface material and the shock wave incident angle. The reflectivity factor of the i-th shell element is calculated through the following formula : ; In the formula, represents the incident angle of the shock wave with the ground surface of the i-th shell element; represents the reflection coefficient of different ground surface materials, specifically: when it is soil ; when it is concrete ; when it is vegetation ; S23. Characterize the proportion of energy lost by shock waves due to absorption and scattering in different media, and obtain the energy dissipation factor of the i-th shell element through the following formula : ; In the formula, represents the absorption coefficient of the base medium, including: when it is soil ; when it is concrete ; when it is vegetation ; represents the cover density correction coefficient, including: when it is bare land ; when it is grassland ; when it is vegetation .

[0013] Furthermore, the associated unit is used to extract the shock wave attenuation factor , reflectivity factor and energy dissipation factor of the i-th shell element in the environmental response dataset. After dimensionless processing, the environmental impact index of the i-th shell element is obtained through the following associated formula : ; In the formula, , and respectively represent the weight coefficients of the shock wave attenuation factor , reflectivity factor and energy dissipation factor of the i-th shell element, and the sum of the weights is 1.

[0014] Furthermore, the second correction module further includes a second correction unit, a simulation plastic deformation output unit and a warning unit; The second correction unit is used to correct the environmental impact index of the i-th shell element with respect to the equivalent stress of the i-th shell element to obtain the corrected equivalent stress of the i-th shell element: ; The simulation plastic deformation output unit is used to identify the brittle fracture point area through the following steps, including: S11. Preset the fracture strength threshold X. When the corrected equivalent stress of the i-th shell element ≥ the fracture strength threshold X, it is marked as a brittle fracture point candidate unit and the "high stress fracture" label is added; Among them, the fracture strength threshold X is obtained from experiments or material manuals, including setting the fracture strength of AISI 4340 to 1300–1500 MPa; S12. When ≥1, it indicates fracture caused by high strain rate explosion, and the unit is marked as a candidate unit for brittle fracture point, with the label "velocity instability fracture" added; S13. When 0.75 ≤ the normalized temperature value of the i-th shell element <1, it indicates that the material of the i-th shell element has a melting and softening phenomenon, and the unit is marked as a candidate unit for brittle fracture point, with the label "thermal softening fracture" added; S14. When any of the conditions in S11 - S13 is met, the i-th shell element is marked as the "brittle fracture point area" and the information is recorded in the database, including: unit number i, timestamp, spatial position coordinates, and the corresponding label; S15. Count the total number of brittle fracture point areas and generate the corresponding damage level, including: When the total number of points in the brittle fracture point area in the three-dimensional internal structure model of the target < 5, it indicates the appearance of minor brittle fracture points, the overall structure is intact, and the spatial distribution of brittle fracture points is discrete, generating a first-level damage level; When 5 ≤ the total number of points in the brittle fracture point area in the three-dimensional internal structure model of the target < 15, it indicates local brittle fracture, the local strength of the structure decreases, and there is a local aggregation area of brittle fracture points, generating a second-level damage level; When the total number of points in the brittle fracture point area in the three-dimensional internal structure model of the target ≥ 15, it indicates multi-point concentrated brittle fracture, and the main structure is damaged and concentrated on the main beam or support surface, generating a third-level damage level; And count the total number of brittle fracture point areas. If the total number of brittle fracture point areas forms a longitudinal or transverse continuous adjacent chain path, and the total number of points in the brittle fracture point area in the three-dimensional internal structure model of the target ≥ 15, it indicates that the brittle fracture point path is connected to form a chain path, there is a risk of fault or crack penetration, and the structure has failed, generating a fourth-level damage level.

[0015] Furthermore, the warning unit is used to generate corresponding warning instructions by combining the corrected first evaluation result and the corresponding damage level, including: Preset a damage quantity threshold L, and count the total number Nn of the i-th shell elements in the high-pressure area. When Nn < L, it indicates that it is within the tolerable range, continuous monitoring is carried out, no warning is triggered, and only the status is recorded; When Nn ≥ L and the damage level is the first-level damage level, it indicates that the overall damage of the target still belongs to "minor brittle fracture", and there is a risk of local strength decline, generating a first warning instruction; When Nn ≥ L and the damage level is the second-level damage level, it indicates that there is a risk of local strength decline, and the risk trend shows an increasing state compared with the first warning instruction, generating a second warning instruction; When Nn≥L, and the damage level is level 3, it indicates that there is an obvious damage trend of "local intensity reduction", and the third warning instruction is generated; When Nn≥L and the damage level is level four, it means that the structure has suffered "concentrated damage" or "connected cracking", and the fourth warning instruction is generated.

[0016] Compared with the prior art, the beneficial effects of the present invention are as follows: the system obtains high-precision geometric information of the target structure through three-dimensional scanning and finite element mesh technology, and constructs a fine physical simulation domain, which can accurately restore the force behavior and stress distribution state of the structure in the actual blasting environment. Physical parameters such as explosive particle size distribution parameters, terrain geometry factors, shock wave reflection and energy dissipation are introduced to construct a double-layer correction module to correct the explosive pressure field and equivalent stress, breaking through the limitation of traditional simulation models that ignore environmental complexity and avoiding evaluation distortion or misjudgment. The system can automatically identify the vulnerable areas of stress concentration in the structure during the simulation process, dynamically calibrate the damage level distribution map in combination with the corrected pressure field and stress assessment, and then generate hierarchical early warning instructions, providing a scientific basis for engineering protection design, emergency response and risk control. The system can automatically identify the vulnerable areas of stress concentration in the structure during the simulation process, dynamically calibrate the damage level distribution map in combination with the corrected pressure field and stress assessment, and then generate hierarchical early warning instructions, providing a scientific basis for engineering protection design, emergency response and risk control.

[0017] The present invention comprehensively considers the equivalent plastic strain value, normalized temperature value and strain rate change according to the equivalent stress, and describes the nonlinear yield behavior of the material under high temperature and high strain rate through multi-parameter control factors (such as thermal softening index, strain rate sensitivity coefficient, strain hardening index). Compared with the traditional ideal elastic-plastic model, the fitting accuracy of the mechanical response of the material under the extreme loading environment of explosion is significantly improved. The present invention fits the experimental parameters of typical engineering materials such as AISI4340 high-strength steel, 6061-T6 aluminum alloy and Ti-6Al-4V titanium alloy, clarifies the key constant values ​​such as initial yield strength, strain hardening coefficient, strain rate sensitivity coefficient, so as to realize the synchronous simulation of explosion response of multiple materials and multiple regions.

[0018] The first correction unit normalizes the energy release rate of the $i$-th shell element with the standard particle size reference rate and corrects the original pressure field of the $i$-th shell element, so that the finally obtained corrected internal pressure field of the explosive can more accurately reflect the detonation pressure fluctuations caused by particle size refinement or coarsening, thereby improving the ability to depict the shock response in the initial stage of explosion. By comparing the corrected pressure field with the pressure threshold $Y$ again through the pressure evaluation unit, a "first corrected evaluation result" is constructed. Compared with the initial distribution evaluation diagram, the risk of misjudgment caused by particle size distribution differences is eliminated, effectively avoiding problems such as high-pressure mislabeling and low-pressure missed detection, thereby improving the credibility and determination sensitivity of damage level evaluation.

[0019] The present invention also uses a second correction module to collect the topographic and geometric data of the location where the shell element is located in real time through a second acquisition unit, and constructs an environmental response data set including a shock wave attenuation factor, a reflectivity factor, and an energy dissipation factor, which can accurately simulate the complex behavior of shock waves during propagation in different landforms (such as mountains, hard ground, grasslands), and overcome the hypothetical errors caused by ignoring environmental factors in traditional explosion models. BRIEF DESCRIPTION OF THE DRAWINGS

[0020] Figure 1 It is a schematic diagram of the block diagram of the damage effectiveness simulation and evaluation system based on model simulation as a whole for the present invention. DETAILED DESCRIPTION OF THE INVENTION

[0021] To make the objectives, technical solutions, and advantages of the present invention clearer and more understandable, the present invention will be further described in detail below in conjunction with specific embodiments.

[0022] It should be noted that unless otherwise defined, the technical terms or scientific terms used in the present invention should have the ordinary meanings understood by those with ordinary skills in the field to which the present invention belongs. The "first", "second", and similar terms used in the present invention do not indicate any order, quantity, or importance, but are only used to distinguish different components. The terms such as "including" or "comprising" mean that the elements or objects appearing before this word cover the elements or objects listed after this word and their equivalents, without excluding other elements or objects. The terms such as "connected" or "coupled" are not limited to physical or mechanical connections, but may include electrical connections, whether direct or indirect. The terms such as "upper", "lower", "left", "right", etc. are only used to represent relative positional relationships, and when the absolute position of the object being described changes, the relative positional relationship may also change accordingly.

[0023] Embodiment 1: Please refer to Figure 1 , the present invention provides a technical solution: including: A target structure analysis module, which is used to obtain three-dimensional scan data of the physical structure of the target simulated explosion, and based on the finite element mesh, reconstruct the physical field simulation model of the target object, establish a simulation domain and divide it into several shell elements. During the process of simulating gunpowder explosion, construct the internal pressure field of the explosive in the i-th shell element and equivalent stress , and preset a pressure threshold Y, and compare the internal pressure field of the explosive in the i-th shell element with the pressure threshold Y to obtain the first evaluation result. According to the first evaluation result, construct the "initial damage level distribution map"; A first correction module, which is used to collect the particle size distribution parameters of explosive particles, obtain the explosive energy data set, and construct the energy release rate of the i-th shell element correct the internal pressure field of the explosive in the i-th shell element to obtain the corrected internal pressure field of the explosive in the i-th shell element ; and correct the "initial damage level distribution map" to obtain the "corrected damage level distribution map"; A second correction module, which is used to collect the terrain geometric information of the i-th shell element, construct the environmental response data set, and construct the environmental impact index of the i-th shell element , and correct the equivalent stress of the i-th shell element to obtain the corrected equivalent stress of the i-th shell element and evaluate, further identify the brittle fracture point area in the physical field simulation model, and dynamically link the brittle fracture point response with the damage level evaluation standard during the damage simulation process, as well as the corresponding warning instructions.

[0024] In this embodiment, the system obtains high-precision geometric information of the target structure through three-dimensional scanning and finite element mesh technology, and constructs a fine physical simulation domain, which can accurately restore the stress behavior and stress distribution state of the structure in the actual blasting environment. By introducing physical parameters such as explosive particle size distribution parameters, terrain geometric factors, shock wave reflection and energy dissipation, a double-layer correction module is constructed to correct the explosive pressure field and equivalent stress, breaking through the limitation of traditional simulation models that ignore environmental complexity and avoiding evaluation distortion or misjudgment. The system can automatically identify the vulnerable areas with stress concentration in the structure during the simulation process, combine the corrected pressure field and stress evaluation, dynamically calibrate the damage level distribution map, and then generate hierarchical warning instructions, providing a scientific basis for engineering protection design, emergency response and risk control. The system can automatically identify the vulnerable areas with stress concentration in the structure during the simulation process, combine the corrected pressure field and stress evaluation, dynamically calibrate the damage level distribution map, and then generate hierarchical warning instructions, providing a scientific basis for engineering protection design, emergency response and risk control.

[0025] Embodiment 2 Please refer toFigure 1 , the target structure analysis module includes a structure data acquisition unit and a finite element network construction unit: The structural data acquisition unit is used to use 3D-LiDAR or industrial CT scanning equipment to collect three-dimensional scanning data of the target simulated explosion physical structure, and after denoising and unified format conversion, a finite element network construction unit is used to establish a physical field simulation model based on the three-dimensional scanning data using Gmsh or AnsysMesh tools. After importing the three-dimensional scanning data, an initial three-dimensional grid is generated, and a simulation domain is established in the initial three-dimensional network, and the simulation domain includes: a target structure domain, an explosion source domain, an environmental medium domain, and a reflection boundary domain; and the simulation domain is divided into a plurality of shell units, which are marked in the physical field simulation model as: , to Represents the first shell element to the The status mark of each shell element at the time point sd.

[0026] The target structure analysis module also includes an explosion triggering process unit, a blasting pressure field construction unit and a blasting intensity division unit: The explosion triggering process unit is used to preset initial conditions, set the explosive heat value, define the detonation point and detonation speed in the simulation domain, and the detonation speed is set to 6930m / s; The blasting pressure field construction unit is used to construct the internal pressure field of the explosive of the i-th shell unit by applying the JWL state equation. : ; In the formula, and Respectively represent the pressure empirical constants at the initial and later stages of the explosion; obtained through the copper cylinder expansion test: is the first attenuation coefficient, which controls the rapid attenuation at the initial stage of the explosion. is the second attenuation coefficient, which controls the slow attenuation in the late stage of the explosion; For TNT nitroglycerin explosives,

[0027] For Comp-B explosives,

[0028] For RDX explosives,

[0029] The meaning of the formula is: The internal pressure field of the explosive of the i-th shell element As an external load applied to the surface of the shell element, stress wave propagation is generated at each time step to form the internal force response of the structure; the internal pressure field of the explosive of the i-th shell element Describes the instantaneous pressure change generated by the explosion of explosives, reflecting the effects of gas expansion and shock waves; Represents the initial volume of the i-th shell element, Represents the current volume of the i-th shell element, Represents the energy ratio of the explosive, Represents the internal energy per unit mass of the i-th shell element, calculated by the following formula: ; where Q is the explosion value of all explosives in the i-th shell element, in J / kg, Represents the explosive mass of the i-th shell element, Represents the explosive density; The blasting intensity division unit is used to preset the pressure threshold Y and compare the internal pressure field of the explosive in the i-th shell element with the pressure threshold Y to obtain the first evaluation result, including: When the internal pressure field of the explosive in the i-th shell element ≥ pressure threshold Y × 120%, it is marked as the high-pressure area; When pressure threshold Y ≤ the internal pressure field of the explosive in the i-th shell element < pressure threshold Y × 120%, it is marked as the medium-pressure area; When the internal pressure field of the explosive in the i-th shell element < pressure threshold Y, it is marked as the low-pressure area; Traverse all shell elements to obtain their internal pressure fields of explosives and construct an "initial damage level distribution map" according to the first evaluation result.

[0030] In this embodiment, the present invention uses 3D-LiDAR or industrial CT scanning equipment to obtain real structure three-dimensional point cloud data. After noise processing and unified format conversion, combined with Gmsh or AnsysMesh modeling tools, it realizes the construction of a complete simulation domain for the target structure domain, explosion source domain, environmental medium domain, and reflection boundary domain, and divides it into multiple shell elements, significantly improving the authenticity of the structure model and the simulation granularity, and avoiding the error accumulation problem caused by traditional idealized modeling. The system is based on the JWL (Jones-Wilkins-Lee) equation of state, and respectively establishes the evolution process of the explosive pressure field inside the shell elements for typical explosive materials such as TNT, Comp-B, and RDX, which can reflect the intense shock in the initial stage of the explosion and the pressure decay process in the later stage. Thus, the credibility of the simulation of the response of the explosion physical field to the target structure is improved. By setting the pressure threshold Y and establishing multi-segment pressure classification (high-pressure area, medium-pressure area, low-pressure area), the present invention can accurately calibrate each shell element, and then draw the "initial damage level distribution map", laying a foundation for subsequent energy correction, environmental response correction, and dynamic tracking of brittle breakpoints, and realizing the whole-process and hierarchical damage mapping from the occurrence of the explosion to the structural response. The system supports customized input of parameters for different types of explosives, including key variables such as energy density, detonation velocity, explosive density, and initiation position, and can be widely applied to multiple fields such as geological blasting, military blasting tests, tunnel boring simulations, and building demolition safety assessments, with good engineering versatility and algorithm compatibility.

[0031] Embodiment 3 Please refer to Figure 1 , the target structure analysis module further includes an equivalent stress construction unit: The equivalent stress construction unit is used to construct the equivalent stress of the ith shell element under the action of the explosion load : ; Among them, is the equivalent plastic strain value of the ith shell element, represents the equivalent strain rate of the ith shell element, represents the reference equivalent strain rate, which is set to , is the normalized temperature value of the ith shell element, defined as , T represents the current temperature, T 0 represents the reference temperature, T m represents the melting point temperature of the material, represents the strain hardening index, and the value range is 0–1; represents the thermal softening index, represents the initial yield strength of the material, represents the strain hardening coefficient, Denote the strain rate sensitivity coefficient; obtained by experimental fitting, specifically: when the material is AISI4340 high-strength steel, K 1 = 792 MPa; K 2 = 510 MPa; K 3 = 0.014; n 1 = 0.26; m 1 = 1.03; T m = 1793 K, T 0 = 300 K; When the material is 6061-T6 aluminum alloy, K 1 = 324 MPa; K 2 = 114 MPa; K 3 = 0.002; n 1 = 0.42; m 1 = 1.3403; T m = 7751793 K, T 0 = 293 K; When the material is Ti-6Al-4V titanium alloy, K 1 = 1098 MPa; K 2 = 1092 MPa; K 3 = 0.014; n 1 = 0.93; m 1 = 1.1; T m = 1878 K, T 0 = 298 K.

[0032] In this embodiment, the equivalent stress comprehensively considers the equivalent plastic strain value, the normalized temperature value, and the strain rate change, and characterizes the non-linear yield behavior of the material under high temperature and high strain rate through multi-parameter control factors (such as the thermal softening index, the strain rate sensitivity coefficient, and the strain hardening index). Compared with the traditional ideal elastic-plastic model, the fitting accuracy of the mechanical response of the material under the extreme explosion loading environment is significantly improved. The present invention fits the experimental parameters of typical engineering materials such as AISI4340 high-strength steel, 6061-T6 aluminum alloy, and Ti-6Al-4V titanium alloy, and determines the key constant values such as the initial yield strength, the strain hardening coefficient, and the strain rate sensitivity coefficient, so as to realize the synchronous simulation of the explosion responses of multiple materials and multiple regions.

[0033] Example 4 Please refer to Figure 1 , the first correction module includes a first acquisition unit, a first correction unit, and a pressure evaluation unit; The first acquisition unit is used to acquire the particle size distribution parameters of explosive particles and obtain an explosive energy data set, where the explosive energy data set includes: the energy release characteristics under different particle size combinations. The value range of the particle size d is divided into Z intervals, and the z-th interval represents explosive particles with an average particle size of ; For the i-th shell unit, the proportion of each particle size interval is sampled as , combined with the explosive mass of the i-th shell unit, and the energy release rate of the i-th shell unit is constructed through the following formula : ; In the formula, represents the mass proportion of the explosive of explosive particles with an average particle size of , satisfying ; represents the release rate of the explosive with a particle size of per unit mass, with the unit of J / kg·s, and is calculated through the following formula: ; where q represents the total energy release amount of the explosive per unit mass. For TNT nitroglycerin explosive, q = 4.2×10 6 J / kg; for Comp-B explosive, q = 5.5×10 6 J / k; for RDX explosive, q = 5.7×10 6 J / k; represents the energy release time of the explosive particles with a particle size of ; The first correction unit is used to correct the energy release rate of the i-th shell unit for the internal pressure field of the explosive in the i-th shell unit to obtain the corrected internal pressure field of the explosive in the i-th shell unit: ; In the formula, represents the reference energy release rate under standard particle size conditions for normalization processing.

[0034] Further, the pressure evaluation unit is used to re-compare and correct the corrected internal pressure field of the explosive in the i-th shell unit with the pressure threshold Y to obtain the corrected first evaluation result, including: When the corrected internal pressure field of the explosive in the i-th shell unit ≥ the pressure threshold Y×120%, it is marked as the high-pressure area; When the pressure threshold Y ≤ the corrected internal pressure field of the explosive in the i-th shell element < the pressure threshold Y × 120%, it is marked as the medium-pressure area; When the corrected internal pressure field of the explosive in the i-th shell element < the pressure threshold Y, it is marked as the low-pressure area; Traverse all shell elements to obtain the corrected internal pressure field of the explosive inside them , and construct a "damage level correction distribution map" by comparing with the corrected first evaluation result.

[0035] In this embodiment, the first correction unit normalizes the energy release rate of the i-th shell element with the standard particle size reference rate and corrects the original pressure field of the i-th shell element, so that the finally obtained corrected internal pressure field of the explosive can more accurately reflect the detonation pressure fluctuation caused by particle size refinement or coarsening, thereby improving the ability to depict the shock response in the initial stage of explosion. By comparing the corrected pressure field with the pressure threshold Y again through the pressure evaluation unit, a "corrected first evaluation result" is constructed. Compared with the initial distribution evaluation map, the risk of misjudgment caused by particle size distribution differences is eliminated, and problems such as high-pressure mislabeling and low-pressure missed detection are effectively avoided, thereby improving the credibility and determination sensitivity of damage level evaluation.

[0036] Embodiment 5 Please refer to Figure 1 , the second correction module includes a second acquisition unit and an associated unit; The second acquisition unit is used to acquire the topographic geometric information of the i-th shell element and construct an environmental response data set, and the environmental response data set includes the shock wave attenuation factor of the i-th shell element , reflectivity factor and energy dissipation factor ; The acquisition methods of the shock wave attenuation factor , reflectivity factor and energy dissipation factor of the i-th shell element are specifically as follows: S21. During the propagation of the shock wave, it will attenuate due to the influence of the terrain. As the distance from the explosion center to the target position increases, the shock wave intensity gradually weakens. The shock wave attenuation of the i-th shell element is calculated through the following formula : ; Among them, is the distance from the i-th shell element to the explosion source, represents the attenuation index, and its value range is 1-2; When the shock wave encounters different types of ground, a certain proportion of the wave energy will be reflected; the reflectivity is affected by the terrain surface material and the incident angle of the shock wave. The reflectivity factor of the i-th shell element is obtained by the following formula : ; In the formula, represents the incident angle of the shock wave with the ground surface of the i-th shell element; represents the reflection coefficient of different ground surface materials. Specifically, for soil ; for concrete ; for vegetation ; S23 represents the proportion of energy lost by the shock wave due to absorption and scattering in different media. The energy dissipation factor of the i-th shell element is obtained by the following formula : ; In the formula, represents the absorption coefficient of the base medium, including: for soil ; for concrete ; for vegetation ; represents the coverage density correction coefficient, including: for bare land ; for grassland ; for vegetation .

[0037] The associated unit is used to extract the shock wave attenuation factor , reflectivity factor and energy dissipation factor of the i-th shell element in the environmental response dataset. After dimensionless processing, the environmental impact index of the i-th shell element is obtained through the following associated formula : ; In the formula, , and respectively represent the weight coefficients of the shock wave attenuation factor , reflectivity factor and energy dissipation factor of the i-th shell element, and the sum of the weights is 1.

[0038] The second correction module further includes a second correction unit, a simulation plastic deformation output unit, and a warning unit; The second correction unit is used to use the environmental impact index of the i-th shell element for the equivalent stress Make corrections to obtain the corrected equivalent stress of the i-th shell element : ; The simulation plastic deformation output unit is used to identify the brittle fracture point area through the following steps, including: S11. Preset the fracture strength threshold X. When the corrected equivalent stress of the i-th shell element ≥ the fracture strength threshold X, it is marked as a brittle fracture point candidate unit and the "high stress fracture" label is added; Among them, the fracture strength threshold X is obtained from experiments or material manuals, including setting the fracture strength of AISI4340 to 1300–1500 MPa; S12. When ≥ 1, it means fracture caused by high explosion strain rate, and it is marked as a brittle fracture point candidate unit and the "velocity instability fracture" label is added; S13. When 0.75 ≤ the normalized temperature value of the i-th shell element <1, it means that there is a melting and softening phenomenon in the material of the i-th shell element, and it is marked as a brittle fracture point candidate unit and the "thermal softening fracture" label is added; when the normalized temperature value of the i-th shell element = 1, the stress is zeroed, but the softening has already started. Therefore, most materials enter the thermal softening sensitive area at = 0.75, which is the key window before dynamic instability; S14. When any of the conditions S11 - S13 is met, the i-th shell element is marked as the "brittle fracture point area" and the information is recorded in the database, including: element number i, timestamp, spatial position coordinates, and the corresponding label; S15. Count the total number of brittle fracture point areas and generate the corresponding damage level, including: When the total number of points in the brittle fracture point area in the three-dimensional internal structure model of the target object < 5, it means that small brittle fracture points appear, the overall structure is intact, and the spatial distribution of brittle fracture points is discrete, generating a first-level damage level; When 5 ≤ the total number of points in the brittle fracture point area in the three-dimensional internal structure model of the target object < 15, it means local brittle fracture, the local strength of the structure decreases, and there is a local aggregation area of brittle fracture points, generating a second-level damage level; When the total number of points in the brittle fracture point area in the three-dimensional internal structure model of the target object ≥ 15, it means multi-point concentrated brittle fracture, the main structure is damaged and concentrated on the main beam or support surface, generating a third-level damage level; And count the total number of brittle fracture point areas. If the total number of brittle fracture point areas forms a longitudinal or transverse continuous adjacent chain path, and the total number of points in the brittle fracture point area in the three-dimensional internal structure model of the target object ≥ 15, it means that the brittle fracture point path is connected to form a chain path, there is a risk of fault or crack penetration, and the structure has failed, generating a fourth-level damage level.

[0039] In this embodiment, the second correction module collects in real time the topographic and geometric data of the location where the shell unit is located through the second acquisition unit, constructs an environmental response data set including a shock wave attenuation factor, a reflectivity factor, and an energy dissipation factor, and can accurately simulate the complex behavior of shock waves during propagation in different landforms (such as mountains, hard ground, grasslands), overcoming the hypothetical errors caused by ignoring environmental factors in traditional explosion models.

[0040] The associated unit performs dimensionless normalization processing on the terrain influence data and calculates the "environmental influence index" based on preset weights (such as attenuation-dominated terrain, reflection-dominated surface, dissipation-dominated medium), realizing weighted modeling of the differential effects of different types of surface conditions on explosion effects, so that the damage analysis has response sensitivity to actual landform conditions.

[0041] Compared with the wave propagation assumption in an ideal homogeneous medium, the present invention dynamically calculates the distance attenuation, reflection loss, and energy dissipation occurring in the propagation path from the explosion source to the i-th shell unit, making the actually obtained impact value closer to the on-site measurement results, and is particularly suitable for explosion effect assessment in urban environments, forest areas, or military protection areas with complex landforms.

[0042] A combined correction method of a basic medium absorption coefficient and a coverage density correction coefficient is introduced into the module, which can distinguish natural conditions such as bare land, grassland, and vegetation, and can assign independent parameters according to different material types (such as soil, concrete, vegetation), supporting the joint prediction modeling of explosion effects in mixed areas of natural surfaces and artificial facilities. By applying the environmental influence index to the explosion effect result of the i-th shell unit through the second correction unit, the damage level can be automatically corrected, thereby improving the smoothness and spatial resolution of the damage level distribution map at the terrain boundary and avoiding damage level errors or prediction deviations in areas with drastic terrain changes.

[0043] Embodiment 6 Please refer to Figure 1 , the warning unit is used to generate corresponding warning instructions by combining the corrected first evaluation result and the corresponding damage level, including: Preset a damage quantity threshold L, and count the total number Nn of the i-th shell units in the high-pressure area. When Nn < L, it means that it is within the acceptable range, continue to monitor, do not trigger a warning, and only record the status; When Nn ≥ L and the damage level is the first-level damage level, it means that the overall damage of the target still belongs to "minor brittle fracture", and there is a risk of local strength decline. Generate a first warning instruction, including: reducing the explosive filling rate by 5% - 10%, extending the detonation delay by 5% - 10%, and keeping the detonation point unchanged or slightly adjusting it 1 - 2 meters away from the edge of the high-pressure area; When Nn ≥ L and the damage level is the second-level damage level, indicating a risk of local strength decline, with the risk trend showing an increasing state compared to the first warning instruction, a second warning instruction is generated, including: reducing the explosive filling rate by 11% - 15%, extending the detonation delay by 11% - 15%, and keeping the detonation point unchanged or slightly adjusting it 3 - 5 meters away from the edge of the high-pressure area; When Nn ≥ L and the damage level is the third-level damage level, indicating an obvious damage trend of "local strength decline", a third warning instruction is generated, including: reducing the explosive filling rate by 16% - 20%, extending the detonation delay by 16% - 20%, and keeping the detonation point unchanged or slightly adjusting it 6 - 8 meters away from the edge of the high-pressure area; When Nn ≥ L and the damage level is the fourth-level damage level, indicating that "concentrated damage" or "connected fracture penetration" has occurred in the structure, a fourth warning instruction is generated, including: stopping the blasting in this area or reducing the explosive filling rate by more than 40%, switching to an extended or segmented and zoned detonation mode, setting the detonation point 9 meters away from the high-pressure area, and adopting a multi-point dispersed detonation scheme.

[0044] In this embodiment, the warning unit establishes a four-level warning stratification strategy based on the preset damage quantity threshold L and the actual number Nn of high-pressure area shell units, combined with the level of the damage grade, which can accurately reflect different states from minor damage, risk enhancement to structural damage, and realize the dynamic evolution process of "no warning - weak warning - medium warning - strong warning - severe warning". This system not only considers the quantitative characteristics of the internal pressure field of the explosive exceeding the threshold (i.e., the number of high-pressure areas), but also superimposes the specific damage grade as a combined judgment basis, making the warning behavior based on the cross-analysis of multiple states, enhancing the accuracy and scientificity of the judgment of the damage evolution trend. When the number of high-pressure areas Nn < L, the system is identified as "within the structural bearing range", automatically enters the recording and monitoring state without triggering a warning, reducing the possibility of false alarms caused by local minor abnormalities, and improving the fault tolerance and decision-making robustness of the system in a complex explosion environment. Under the condition of Nn ≥ L, the warning unit not only makes a judgment of "whether to give a warning", but also makes a fine distinction of the risk trend based on the damage grade. For example, the first-level warning represents a minor brittle fracture omen, and the fourth-level warning represents the critical state of structural connected damage, realizing the two-dimensional judgment logic of "trend discrimination + grade response" for risks.

[0045] Warning instructions of different levels can be used as trigger signals for subsequent response processes, promoting the automatic initiation of intervention measures such as target structure reinforcement, evacuation suggestions, and subsequent structural integrity monitoring, providing accurate data support for rapid response, disaster mitigation, and protection optimization in an explosion environment.

[0046] It should be noted that all the calculation formulas in this application document adopt regression analysis including but not limited to machine learning algorithms to deeply analyze the relevant parameters collected, identify their natural trends and interrelationships. Using professional software such as the Scikit-learn library of Python or the R language, a mathematical model matching the data is automatically generated. Then, the performance of the model is objectively evaluated through methods such as cross-validation, and combined with continuous feedback and optimization to ensure that the created formula truly reflects the internal laws of the data, thereby ensuring its effectiveness and accuracy, and ensuring that the calculation process conforms to the constraints of natural laws rather than being based on artificially set rules.

[0047] Essentially, the technical solution of the present invention, or the part that contributes to the prior art, can be embodied in the form of a software product. This computer software product can be stored in a computer-readable storage medium, such as a floppy disk of a computer, a read-only memory (ROM), a random access memory (RAM), a flash memory (FLASH), a hard disk or an optical disc, etc., including several instructions for causing a computer device (which can be a personal computer, a server, or a network device, etc.) to execute the methods of various embodiments of the present invention.

[0048] The logic and / or steps represented in the flowchart or described in other ways herein, for example, can be considered as a definite sequence list of executable instructions for implementing logical functions, and can be specifically implemented in any computer-readable medium for use by an instruction execution system, apparatus or device (such as a computer-based system, a system including a processor, or other systems that can fetch instructions from the instruction execution system, apparatus or device and execute the instructions), or used in combination with these instruction execution systems, apparatus or devices. For the purposes of this specification, a "computer-readable medium" can be any device that can contain, store, communicate, propagate or transmit a program for use by or in connection with an instruction execution system, apparatus or device.

[0049] It should be noted that the above embodiments are only used to illustrate the technical solutions of the present invention and not to limit them. Although the present invention has been described in detail with reference to the preferred embodiments, those of ordinary skill in the art should understand that the technical solutions of the present invention can be modified or equivalently replaced without departing from the spirit and scope of the technical solutions of the present invention, and they should all be covered within the scope of the claims of the present invention.

[0050] It should be noted that the above embodiments are only used to illustrate the technical solutions of the present invention and not to limit them. Although the present invention has been described in detail with reference to the preferred embodiments, those of ordinary skill in the art should understand that the technical solutions of the present invention can be modified or equivalently replaced without departing from the spirit and scope of the technical solutions of the present invention, and they should all be covered within the scope of the claims of the present invention.

Claims

1. A damage effectiveness simulation and evaluation system based on model simulation, characterized in that: include: The target structure analysis module is used to obtain the 3D scanning data of the target simulated explosion structure, and reconstruct the physical field simulation model of the target based on the finite element mesh, establish the simulation domain and divide it into several shell units, and construct the internal pressure field of the explosive of the i-th shell unit in the process of simulating gunpowder explosion. and equivalent stress , and preset the pressure threshold Y, and the internal pressure field of the explosive of the i-th shell element Compare with the pressure threshold Y to obtain a first evaluation result, and construct an "initial distribution map of damage level" according to the first evaluation result; The first correction module is used to collect the particle size distribution parameters of explosive particles, obtain the explosive energy data set, and construct the energy release rate of the i-th shell unit The internal pressure field of the explosive for the i-th shell element Make corrections to obtain the corrected internal pressure field of the explosive for the i-th shell element ; and correct the "initial damage level distribution map" to obtain the "corrected damage level distribution map"; The second correction module is used to collect the terrain geometry information of the i-th shell unit, construct the environmental response data set, and construct the environmental impact index of the i-th shell unit. , and the equivalent stress of the i-th shell element Make corrections to obtain the corrected equivalent stress of the i-th shell element And evaluate, further identify the brittle fracture point area in the physical field simulation model, and dynamically link the brittle fracture point response with the damage level assessment standard and the corresponding early warning instructions during the damage simulation process.

2. The damage effectiveness simulation and evaluation system based on model simulation according to claim 1, characterized in that: The target structure analysis module includes a structure data acquisition unit and a finite element network construction unit: The structural data acquisition unit is used to use 3D-LiDAR or industrial CT scanning equipment to collect three-dimensional scanning data of the target simulated explosion physical structure, and after denoising and unified format conversion, a finite element network construction unit is used to establish a physical field simulation model based on the three-dimensional scanning data using Gmsh or AnsysMesh tools. After importing the three-dimensional scanning data, an initial three-dimensional grid is generated, and a simulation domain is established in the initial three-dimensional network, and the simulation domain includes: a target structure domain, an explosion source domain, an environmental medium domain, and a reflection boundary domain; and the simulation domain is divided into a plurality of shell units, which are marked in the physical field simulation model as: , to Represents the first shell element to the The status mark of each shell element at the time point sd.

3. The damage effectiveness simulation and evaluation system based on model simulation according to claim 2, characterized in that: The target structure analysis module also includes an explosion triggering process unit, a blasting pressure field construction unit and a blasting intensity division unit: The explosion triggering process unit is used to preset initial conditions, set the explosive heat value, define the detonation point and detonation speed in the simulation domain, and the detonation speed is set to 6930m / s; The blasting pressure field construction unit is used to construct the internal pressure field of the explosive of the i-th shell unit by applying the JWL state equation. : ； In the formula, and Respectively represent the pressure empirical constants at the initial and later stages of the explosion; obtained through the copper cylinder expansion test: is the first attenuation coefficient, which controls the rapid attenuation at the initial stage of the explosion. is the second attenuation coefficient, which controls the slow attenuation in the late stage of the explosion; represents the initial volume of the i-th shell element, represents the current volume of the i-th shell element, Represents the energy ratio of the explosive, represents the internal energy per unit mass of the ith shell element, which is calculated by the following formula: ； Where Q is the explosion value of all explosives in the i-th shell element, in J / kg, represents the explosive mass of the i-th shell element, Indicates the density of explosives; The blasting intensity division unit is used to preset the pressure threshold Y and divide the internal pressure field of the explosive of the i-th shell unit into Compare with the pressure threshold Y to obtain a first evaluation result, including: When the internal pressure field of the explosive of the i-th shell element ≥pressure threshold Y×120%, marked as high pressure area; When the pressure threshold Y≤the internal pressure field of the explosive of the i-th shell element <pressure threshold Y×120%, marked as medium pressure area; When the internal pressure field of the explosive of the i-th shell element <pressure threshold Y, marked as low pressure area; Traverse all shell elements to obtain the internal pressure field of the explosive , and construct an "initial distribution map of damage levels" by comparing the first assessment results.

4. The damage effectiveness simulation and evaluation system based on model simulation according to claim 3 is characterized in that: The target structure analysis module also includes an equivalent stress construction unit: The equivalent stress construction unit is used to construct the equivalent stress of the i-th shell element under the explosion load. : ； in, is the equivalent plastic strain value of the i-th shell element, represents the equivalent strain rate of the i-th shell element, represents the reference equivalent strain rate, set to , is the normalized temperature value of the ith shell element, defined as, T represents the current temperature, T0 represents the reference temperature, T m Indicates the melting point temperature of the material. represents the strain hardening exponent, It represents the thermal softening index. represents the initial yield strength of the material, represents the strain hardening coefficient, represents the strain rate sensitivity coefficient.

5. The damage effectiveness simulation and evaluation system based on model simulation according to claim 4 is characterized in that: The first correction module includes a first acquisition unit, a first correction unit and a pressure assessment unit; The first acquisition unit is used to collect the particle size distribution parameters of the explosive particles to obtain the explosive energy data set, which includes: energy release characteristics under different particle size combinations, dividing the value range of the particle size d into Z intervals, and the zth interval represents the average particle size of of explosive particles; For the i-th shell element, the proportion of each particle size interval sampled is , combined with the explosive mass of the i-th shell element , the energy release rate of the i-th shell element is constructed by the following formula : ； In the formula, The average particle size is The mass ratio of explosive particles to explosives satisfies ; The particle size per unit mass is The explosive release rate, in J / kg·s, is calculated using the following formula: ； Where q represents the total energy released per unit mass of explosive; The particle size is The energy release time of the explosive particles; The first correction unit is used to convert the energy release rate of the i-th shell element The internal pressure field of the explosive for the i-th shell element Make corrections to obtain the corrected internal pressure field of the explosive for the i-th shell element : ； In the formula, It represents the reference energy release rate under standard particle size conditions and is used for normalization.

6. The damage effectiveness simulation and evaluation system based on model simulation according to claim 5 is characterized in that: The pressure evaluation unit is used to convert the corrected internal pressure field of the explosive of the i-th shell element into Re-comparing and correcting with the pressure threshold value Y to obtain a corrected first evaluation result, including: When the corrected internal pressure field of the explosive of the i-th shell element ≥pressure threshold Y×120%, marked as high pressure area; When the pressure threshold Y≤the corrected internal pressure field of the explosive of the i-th shell element <pressure threshold Y×120%, marked as medium pressure area; When the corrected internal pressure field of the explosive of the i-th shell element <pressure threshold Y, marked as low pressure area; Traverse all shell elements to obtain the corrected internal pressure field of the explosive inside , and construct a "damage level correction distribution map" by comparing it with the corrected first assessment results.

7. The damage effectiveness simulation and evaluation system based on model simulation according to claim 6 is characterized in that: The second correction module includes a second acquisition unit and an associated unit; The second acquisition unit is used to acquire the terrain geometry information of the i-th shell element and construct an environmental response data set, wherein the environmental response data set includes the shock wave attenuation factor of the i-th shell element. , reflectivity factor and energy dissipation factor ; Shock wave attenuation factor of the i-th shell element , reflectivity factor and energy dissipation factor The specific way to obtain is: S21. The shock wave will attenuate due to the influence of the terrain during propagation. As the distance from the explosion center to the target location increases, the shock wave intensity gradually weakens. The shock wave attenuation of the i-th shell element is calculated by the following formula: : ； in, is the distance from the ith shell element to the explosion source, Represents the attenuation index, with a value range of 1-2; S22, reflectivity is affected by the terrain surface material and the shock wave incident angle. The reflectivity factor of the i-th shell element is calculated by the following formula: : ； In the formula, represents the incident angle of the shock wave with the surface of the i-th shell element; Represents the reflection coefficient of different surface materials; S23, characterizes the proportion of energy lost by shock waves due to absorption and scattering in different media, and calculates the energy dissipation factor of the i-th shell element using the following formula : ； In the formula, represents the basic medium absorption coefficient; Indicates the coverage density correction factor.

8. The damage effectiveness simulation and evaluation system based on model simulation according to claim 7 is characterized in that: The associated unit is used to extract the shock wave attenuation factor of the i-th shell element in the environmental response data set. , reflectivity factor and energy dissipation factor , after dimensionless processing, the environmental impact index of the i-th shell element is obtained by the following correlation formula : ； In the formula, , and are the shock wave attenuation factors of the i-th shell element. , reflectivity factor and energy dissipation factor The weight coefficient of , and the weight sum is 1.

9. The damage effectiveness simulation and evaluation system based on model simulation according to claim 8, characterized in that: The second correction module also includes a second correction unit, a simulated plastic deformation output unit and an early warning unit; The second correction unit is used to adjust the environmental impact index of the i-th shell element Equivalent stress on the i-th shell element Make corrections to obtain the corrected equivalent stress of the i-th shell element : ； The simulated plastic deformation output unit is used to identify the brittle fracture point area through the following steps, including: S11, preset fracture strength threshold X, when the corrected equivalent stress of the i-th shell element ≥ fracture strength threshold X, it is marked as a brittle fracture candidate unit and a "high stress fracture" label is added; S12, when When ≥1, it indicates that the explosion caused high strain rate fracture, and the unit is marked as a brittle fracture point candidate unit, and the "velocity instability fracture" label is added; S13, when 0.75≤normalized temperature value of the i-th shell element When <1, it means that the material of the i-th shell element has melt softening phenomenon, and it is marked as a brittle fracture candidate element, and the "thermal softening fracture" label is added; S14. When any of the conditions S11-S13 is met, the i-th shell element is marked as a "brittle fracture point area" and the information is recorded in the database, including: element number i, timestamp, spatial position coordinates and corresponding labels; S15. Count the total number of brittle fracture point areas and generate corresponding damage levels, including: When the total number of brittle fracture points in the three-dimensional internal structure model of the target object is less than 5, it means that tiny brittle fracture points appear, the overall structure is intact, and the spatial distribution of brittle fracture points is discrete, generating a first-level damage level; When 5≤the total number of brittle fracture point areas in the three-dimensional internal structure model of the target object is less than 15, it indicates local brittle fracture, local strength of the structure decreases, and there is a local clustering area of ​​brittle fracture points, generating a secondary damage level; When the total number of brittle fracture point areas in the three-dimensional internal structure model of the target object is ≥15, it indicates multiple concentrated brittle fractures, the main structure is damaged and concentrated on the main beam or supporting surface, and a third-level damage level is generated; The total number of brittle fracture point areas is counted. If the total number of brittle fracture point areas forms a continuous adjacent chain path in the longitudinal or transverse direction, and the total number of brittle fracture point areas in the three-dimensional internal structure model of the target object is ≥15, it means that the brittle fracture point paths are connected to form a chain path, there is a risk of fault or crack penetration, the structure has failed, and a fourth-level damage level is generated.

10. The damage effectiveness simulation and evaluation system based on model simulation according to claim 9, characterized in that: The early warning unit is used to generate a corresponding early warning instruction in combination with the corrected first assessment result and the corresponding damage level, including: The damage threshold L is preset, and the total number Nn of shell elements in the high-pressure area is counted. When Nn < L, it means that it is within the tolerance range, and continuous monitoring is performed without triggering an early warning, and only the status is recorded; When Nn≥L, and the damage level is level 1, it means that the overall damage of the target is still "minor brittle fracture", and there is a risk of local strength reduction, and the first warning instruction is generated; When Nn≥L, and the damage level is the second level, it means there is a risk of local strength reduction, and the risk trend is increasing compared to the first warning instruction, so a second warning instruction is generated; When Nn≥L, and the damage level is level 3, it indicates that there is an obvious damage trend of "local intensity reduction", and the third warning instruction is generated; When Nn≥L and the damage level is level 4, it means that the structure has experienced "concentrated damage" or "connected cracking", and the fourth warning instruction is generated.

Images