Damage effectiveness simulation and evaluation system based on model simulation

Abstract

The present invention provides a damage effectiveness simulation and evaluation system based on model simulation, which relates to the technical field of explosive action simulation analysis. The present invention obtains three-dimensional scanning data of the target simulated explosion physical structure and reconstructs the physical field simulation model of the target object based on the finite element grid, which can accurately restore the force behavior and stress distribution state of the structure in the actual blasting environment. The explosive particle size distribution parameters, terrain geometry factors, shock wave reflection and energy dissipation physical parameters, the first correction and the second correction modules are introduced to correct the explosive pressure field and equivalent stress, breaking through the limitation of traditional simulation models that ignore the complexity of the environment and avoiding evaluation distortion or misjudgment. During the simulation process, the fragile areas of stress concentration in the structure are automatically identified, and the damage level distribution map is dynamically calibrated in combination with the corrected pressure field and stress evaluation, thereby generating hierarchical early warning instructions, providing a scientific basis 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 in particular to a damage effectiveness simulation and evaluation system based on model simulation. Background Art

[0002] In scenarios such as mining, tunneling, foundation blasting, and military offensive and defensive operations, blasting operations often take place in complex environments, often characterized by rugged terrain, compact target structures, reflective obstacles, or inhomogeneous media. These environmental factors can significantly influence the propagation of explosive energy, thereby altering the stress state and damage pattern of the target structure.

[0003] Existing blasting simulation assessment techniques are typically based on uniformly simplified models, ignoring the complexity of the terrain and its impact on shock wave propagation paths, energy attenuation rates, reflection superposition effects, and the formation of localized overpressure zones. For example, mountains, walls, or underground structures can trigger secondary reflections and focusing of shock waves, resulting in localized structural stresses far exceeding theoretical estimates and causing "unintended damage." In summary, existing techniques neglect the reflection, diffraction, and focusing of shock waves after encountering obstacles (such as walls and rock masses), making it impossible to predict the localized high-pressure damage zones formed by energy accumulation, leading to insufficient accuracy.

[0004] The above information disclosed in this Background section is only for enhancement of understanding of the background of the present disclosure and therefore it may contain information that does not form the prior art that is already known to a person 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 technology.

[0006] To achieve the above object, the present invention provides the following technical solutions:

[0007] The damage effectiveness simulation and evaluation system based on model simulation includes:

[0008] 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 elements, and construct the internal pressure field of the explosive of the i-th shell element in the process of simulating the 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 assessment result, and construct an "initial damage level distribution map" based on the first assessment result;

[0009] 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";

[0010] 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 and damage level assessment standards, as well as the corresponding early warning instructions during the damage simulation process.

[0011] Furthermore, the target structure analysis module includes a structure data acquisition unit and a finite element network construction unit:

[0012] 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, use the finite element network construction unit 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. 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 elements, 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.

[0013] 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:

[0014] 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;

[0015] The blasting pressure field construction unit is used to construct the internal pressure field of the explosive of the i-th shell element by applying the JWL state equation. :

[0016] ;

[0017] Where, and Represent the empirical pressure constants at the initial and later stages of the explosion, respectively; 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;

[0018] For TNT nitroglycerin explosives, ;

[0019] For Comp-B explosives, ;

[0020] For RDX explosives, ;

[0021] 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, forming 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 caused by explosive explosion, 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, Expresses the energy ratio of explosives, represents the internal energy per unit mass of the i-th shell element, which is calculated by the following formula:

[0022] ;

[0023] 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;

[0024] 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:

[0025] When the internal pressure field of the explosive of the i-th shell element ≥pressure threshold Y×120%, marked as high pressure area;

[0026] 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;

[0027] When the internal pressure field of the explosive of the i-th shell element <pressure threshold Y, marked as low pressure area;

[0028] Traverse all shell elements to obtain the internal pressure field of the explosive , and construct an “initial distribution map of damage levels” by comparing it with the first assessment results.

[0029] Furthermore, the target structure analysis module also includes an equivalent stress construction unit:

[0030] The equivalent stress construction unit is used to construct the equivalent stress of the i-th shell element under the action of the explosion load. :

[0031] ;

[0032] 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, represents the thermal softening index, represents the initial yield strength of the material, represents the strain hardening coefficient, represents the strain rate sensitivity coefficient; obtained through experimental fitting, specifically:

[0033] When the material is AISI4340 high-strength steel, K1=792MPa; K2=510MPa; K3=0.014; n1=0.26; m1=1.03; T m =1793K, T0=300K

[0034] When the material is 6061-T6 aluminum alloy, K1=324MPa; K2=114MPa; K3=0.002; n1=0.42; m1=1.3403; T m =7751793K, T0=293K

[0035] When the material is Ti-6Al-4V titanium alloy, K1=1098MPa; K2=1092MPa; K3=0.014; n1=0.93; m1=1.1; T m=1878K, T0=298K.

[0036] Furthermore, the first correction module includes a first acquisition unit, a first correction unit and a pressure evaluation unit;

[0037] 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 explosive particles;

[0038] 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 :

[0039] ;

[0040] Where, Indicates the average particle size The mass ratio of explosive particles to explosives satisfies ; Indicates the particle size per unit mass is The explosive release rate, in J / kg·s, is calculated using the following formula:

[0041] ;

[0042] Where q represents the total energy released per unit mass of explosives. For TNT nitroglycerin explosives, 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; Indicates particle size The energy release time of the explosive particles;

[0043] 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 :

[0044] ;

[0045] Where, It represents the reference energy release rate under standard particle size conditions and is used for normalization.

[0046] Furthermore, the pressure evaluation unit is used to convert the corrected explosive internal pressure field of the i-th shell element into Re-comparing and correcting with the pressure threshold Y to obtain a corrected first assessment result, including:

[0047] When the corrected internal pressure field of the explosive of the i-th shell element ≥pressure threshold Y×120%, marked as high pressure area;

[0048] 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;

[0049] When the corrected internal pressure field of the explosive of the i-th shell element <pressure threshold Y, marked as low pressure area;

[0050] Traverse all shell elements to obtain the corrected internal pressure field of the explosive , and construct a "damage level correction distribution map" by comparing it with the corrected first assessment results.

[0051] Furthermore, the second correction module includes a second acquisition unit and an association unit;

[0052] The second acquisition unit is used to collect the terrain geometry 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 ;

[0053] Shock wave attenuation factor of the i-th shell element , reflectivity factor and energy dissipation factor The specific way to obtain it is:

[0054] S21. The shock wave will attenuate during its propagation due to the influence of the terrain. 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 using the following formula: :

[0055] ;

[0056] in, is the distance from the ith shell element to the explosion source, Represents the attenuation index, with a value range of 1-2;

[0057] 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 by the following formula: :

[0058] ;

[0059] Where, represents the incident angle of the shock wave on the surface of the i-th shell element; Represents the reflection coefficient of different surface materials, specifically: soil ; When concrete ; Vegetation ;

[0060] 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 :

[0061]

[0062] Where, Indicates the basic medium absorption coefficient, including: soil ; When concrete ; Vegetation ;

[0063] Indicates the coverage density correction factor, including: when bare ground Grassland ; Vegetation .

[0064] Furthermore, 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 :

[0065] ;

[0066] Where, 、 and are the shock wave attenuation factors of the i-th shell element. , reflectivity factor and energy dissipation factor The weight coefficients are , and the sum of the weights is 1.

[0067] Furthermore, the second correction module further includes a second correction unit, a simulated plastic deformation output unit and an early warning unit;

[0068] 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 :

[0069] ;

[0070] The simulated plastic deformation output unit is used to identify the brittle fracture point area through the following steps, including:

[0071] 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 the “high stress fracture” label is added;

[0072] The fracture strength threshold X is obtained from experiments or material manuals, including AISI4340, which is set to 1300–1500 MPa;

[0073] 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 label "velocity instability fracture" is added;

[0074] 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 added with the label of "thermal softening fracture";

[0075] S14. When any of the conditions S11-S13 are met, the i-th shell element is marked as a “brittle fracture point region” and the information is recorded in the database, including: element number i, timestamp, spatial position coordinates and corresponding label;

[0076] S15. Count the total number of brittle fracture point areas and generate corresponding damage levels, including:

[0077] 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 level 1 damage level;

[0078] When 5≤the total number of brittle fracture points in the three-dimensional internal structure model of the target object is less than 15, it indicates local brittle fracture, local strength reduction of the structure, and the presence of localized clustered areas of brittle fracture points, generating a secondary damage level;

[0079] When the total number of brittle fracture points in the target's 3D internal structure model is ≥15, it indicates multiple concentrated brittle fractures, and the main structure is damaged and concentrated on the main beam or supporting surface, generating a third level of damage level;

[0080] 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.

[0081] Furthermore, the warning unit is configured to generate a corresponding warning instruction based on the corrected first assessment result and the corresponding damage level, including:

[0082] A damage threshold L is preset, and the total number Nn of shell elements in the i-th high-pressure area is counted. When Nn < L, it means that the damage is within the tolerance range, and continuous monitoring is performed without triggering an early warning, and only the status is recorded.

[0083] When Nn ≥ L and the damage level is level 1, it means that the overall damage to the target is still classified as "minor brittle fracture" and there is a risk of local strength reduction, and the first warning command is generated;

[0084] When Nn ≥ L and the damage level is level 2, it indicates that there is a risk of local intensity reduction, and the risk trend is increasing compared to the first warning instruction, so a second warning instruction is generated;

[0085] 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;

[0086] When Nn≥L and the damage level is level four, it means that the structure has experienced "concentrated damage" or "connected crack penetration" and the fourth warning instruction is generated.

[0087] Compared with existing technologies, the present invention offers the following advantages: The system uses three-dimensional scanning and finite element meshing technology to acquire high-precision geometric information about the target structure and construct a detailed physical simulation domain, accurately reproducing the structure's mechanical behavior and stress distribution in an actual blasting environment. By incorporating physical parameters such as explosive particle size distribution, terrain geometry, shock wave reflection, and energy dissipation, a dual-layer correction module is constructed to correct the explosive pressure field and equivalent stress, overcoming the limitations of traditional simulation models that ignore environmental complexity and avoiding evaluation distortion or misjudgment. During the simulation process, the system automatically identifies vulnerable areas of stress concentration within the structure. Combining the corrected pressure field with stress assessments, it dynamically calibrates the damage level distribution map, generating hierarchical early warning instructions, providing a scientific basis for engineering protection design, emergency response, and risk control.

[0088] This method comprehensively considers equivalent plastic strain, normalized temperature, and strain rate variations based on equivalent stress. It characterizes the nonlinear yield behavior of materials at high temperatures and high strain rates through multi-parameter control factors (such as the thermal softening index, strain rate sensitivity coefficient, and strain hardening exponent). Compared to traditional ideal elastic-plastic models, this method significantly improves the accuracy of fitting the mechanical response of materials under extreme explosive loading environments. By fitting the experimental parameters of typical engineering materials such as AISI4340 high-strength steel, 6061-T6 aluminum alloy, and Ti-6Al-4V titanium alloy, the present invention clarifies the values ​​of key constants such as initial yield strength, strain hardening coefficient, and strain rate sensitivity coefficient, thereby enabling simultaneous simulation of explosion responses across multiple materials and regions.

[0089] The first correction unit normalizes the energy release rate of the i-th shell element with the reference rate of the standard particle size and corrects the original pressure field of the i-th shell element. This results in a corrected internal pressure field of the explosive that more accurately reflects the detonation pressure fluctuations caused by particle size refinement or coarsening, thereby improving the ability to depict the initial shock response of the explosion. The pressure assessment unit compares the corrected pressure field again with the pressure threshold Y to construct a "corrected first assessment result." Compared to the initial distribution assessment map, this eliminates the risk of misjudgment caused by differences in particle size distribution, effectively avoiding problems such as high-pressure mislabeling and low-pressure missed detection, thereby improving the credibility and sensitivity of the damage level assessment.

[0090] The present invention also uses the second correction module to collect the terrain geometry data of the shell unit's location in real time through the second acquisition unit, and constructs an environmental response data set including shock wave attenuation factor, reflectivity factor and energy dissipation factor. It can accurately simulate the complex behavior of shock waves during propagation in different landforms (such as mountains, hard ground, and grassland), and overcome the hypothetical error of ignoring environmental factors in traditional explosion models. BRIEF DESCRIPTION OF THE DRAWINGS

[0091] Figure 1 It is a schematic block diagram of the damage effectiveness simulation and evaluation system based on model simulation of the present invention. DETAILED DESCRIPTION

[0092] In order to make the objectives, technical solutions and advantages of the present invention more clearly understood, the present invention is further described in detail below with reference to specific embodiments.

[0093] It should be noted that, unless otherwise defined, the technical or scientific terms used in the present invention should have the usual meanings understood by people with ordinary skills in the field to which the present invention belongs. The "first", "second" and similar words used in the present invention do not indicate any order, quantity or importance, but are only used to distinguish different components. "Include" or "comprise" and similar words mean that the elements or objects appearing before the word include the elements or objects listed after the word and their equivalents, without excluding other elements or objects. "Connect" or "connected" and similar words are not limited to physical or mechanical connections, but may include electrical connections, whether direct or indirect. "Up", "down", "left", "right" and the like are only used to indicate relative position relationships. When the absolute position of the object being described changes, the relative position relationship may also change accordingly.

[0094] Example 1:

[0095] See also Figure 1 , the present invention provides a technical solution: including:

[0096] 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 elements, and construct the internal pressure field of the explosive of the i-th shell element in the process of simulating the 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 assessment result, and construct an "initial damage level distribution map" based on the first assessment result;

[0097] 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";

[0098] 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 and damage level assessment standards, as well as the corresponding early warning instructions during the damage simulation process.

[0099] In this embodiment, the system uses 3D scanning and finite element meshing technology to acquire high-precision geometric information of the target structure and construct a detailed physical simulation domain, accurately reproducing the structural force behavior and stress distribution in an actual blasting environment. Physical parameters such as explosive particle size distribution, terrain geometry, shock wave reflection, and energy dissipation are incorporated into a two-layer correction module to correct the explosive pressure field and equivalent stress. This overcomes the limitations of traditional simulation models that ignore environmental complexity and avoids evaluation distortion or misjudgment. During the simulation process, the system automatically identifies vulnerable areas of stress concentration within the structure. Combining the corrected pressure field with stress assessments, it dynamically calibrates the damage level distribution map, generating hierarchical early warning instructions, providing a scientific basis for engineering protection design, emergency response, and risk control.

[0100] Example 2

[0101] See also Figure 1 , the target structure analysis module includes a structure data acquisition unit and a finite element network construction unit:

[0102] 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, use the finite element network construction unit 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. 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 elements, 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.

[0103] The target structure analysis module also includes an explosion triggering process unit, a blasting pressure field construction unit, and a blasting intensity division unit:

[0104] 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;

[0105] The blasting pressure field construction unit is used to construct the internal pressure field of the explosive of the i-th shell element by applying the JWL state equation. :

[0106] ;

[0107] Where, and Represent the empirical pressure constants at the initial and later stages of the explosion, respectively; 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;

[0108] For TNT nitroglycerin explosives,

[0109] For Comp-B explosives,

[0110] For RDX explosives,

[0111] 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, forming 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 caused by explosive explosion, 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, Expresses the energy ratio of explosives, represents the internal energy per unit mass of the i-th shell element, which is calculated by the following formula:

[0112] ;

[0113] 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;

[0114] 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:

[0115] When the internal pressure field of the explosive of the i-th shell element ≥pressure threshold Y×120%, marked as high pressure area;

[0116] 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;

[0117] When the internal pressure field of the explosive of the i-th shell element <pressure threshold Y, marked as low pressure area;

[0118] Traverse all shell elements to obtain the internal pressure field of the explosive , and construct an “initial distribution map of damage levels” by comparing it with the first assessment results.

[0119] In this embodiment, the present invention uses 3D-LiDAR or industrial CT scanning equipment to acquire 3D point cloud data of real structures. After noise processing and unified format conversion, it is combined with Gmsh or Ansys Mesh modeling tools to construct a complete simulation domain consisting of the target structure, explosion source, ambient medium, and reflective boundary. This domain is then divided into multiple shell elements, significantly improving the realism and granularity of the structural model and avoiding the error accumulation problems associated with traditional idealized modeling. Based on the JWL (Jones-Wilkins-Lee) equation of state, the system models the evolution of the explosive pressure field within the shell elements for typical explosive materials such as TNT, Comp-B, and RDX. This model can reflect the initial violent impact of the explosion and the subsequent pressure decay process. This improves the reliability of the simulation of the target structure's response to the explosion physics field. By setting a pressure threshold Y and establishing a multi-segment pressure classification system (high, medium, and low pressure zones), the present invention accurately calibrates each shell element, generating an "initial damage distribution map." This lays the foundation for subsequent energy correction, environmental response correction, and dynamic tracking of brittle fracture points, achieving a comprehensive, hierarchical damage mapping process from explosion onset to structural response. The system supports customized parameter input for different explosive types, including key variables such as energy density, detonation velocity, explosive density, and detonation location. It is widely applicable to various fields, including geological blasting, military blasting testing, tunneling simulation, and building demolition safety assessment, demonstrating excellent engineering versatility and algorithm compatibility.

[0120] Example 3

[0121] See also Figure 1 , the target structure analysis module also includes an equivalent stress construction unit:

[0122] The equivalent stress construction unit is used to construct the equivalent stress of the i-th shell element under the action of the explosion load. :

[0123] ;

[0124] 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, with a value range of 0–1; represents the thermal softening index, represents the initial yield strength of the material, represents the strain hardening coefficient, represents the strain rate sensitivity coefficient; obtained through experimental fitting, specifically: when the material is AISI4340 high-strength steel, K1=792MPa; K2=510MPa; K3=0.014; n1=0.26; m1=1.03; T m =1793K, T0=300K;

[0125] When the material is 6061-T6 aluminum alloy, K1=324MPa; K2=114MPa; K3=0.002; n1=0.42; m1=1.3403; T m =7751793K, T0=293K;

[0126] When the material is Ti-6Al-4V titanium alloy, K1=1098MPa; K2=1092MPa; K3=0.014; n1=0.93; m1=1.1; T m =1878K, T0=298K.

[0127] In this embodiment, the equivalent stress comprehensively considers the equivalent plastic strain value, normalized temperature value, and strain rate variation. Multi-parameter control factors (such as the thermal softening index, strain rate sensitivity coefficient, and strain hardening exponent) are used to characterize the nonlinear yield behavior of the material under high temperature and high strain rate. Compared with the traditional ideal elastic-plastic model, this significantly improves the accuracy of fitting the mechanical response of the material under extreme explosive loading environments. By fitting the experimental parameters of typical engineering materials such as AISI4340 high-strength steel, 6061-T6 aluminum alloy, and Ti-6Al-4V titanium alloy, the present invention clarifies the values ​​of key constants such as initial yield strength, strain hardening coefficient, and strain rate sensitivity coefficient, thereby enabling simultaneous simulation of explosion response in multiple materials and regions.

[0128] Example 4

[0129] See also Figure 1 , the first correction module includes a first acquisition unit, a first correction unit and a pressure evaluation unit;

[0130] 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 explosive particles;

[0131] 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 :

[0132] ;

[0133] Where, Indicates the average particle size The mass ratio of explosive particles to explosives satisfies ; Indicates the particle size per unit mass is The explosive release rate, in J / kg·s, is calculated using the following formula:

[0134] ;

[0135] Where q represents the total energy released per unit mass of explosives. For TNT nitroglycerin explosives, 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; Indicates particle size The energy release time of the explosive particles;

[0136] 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 :

[0137] ;

[0138] Where, It represents the reference energy release rate under standard particle size conditions and is used for normalization.

[0139] Furthermore, the pressure evaluation unit is used to convert the corrected explosive internal pressure field of the i-th shell element into Re-comparing and correcting with the pressure threshold Y to obtain a corrected first assessment result, including:

[0140] When the corrected internal pressure field of the explosive of the i-th shell element ≥pressure threshold Y×120%, marked as high pressure area;

[0141] 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;

[0142] When the corrected internal pressure field of the explosive of the i-th shell element <pressure threshold Y, marked as low pressure area;

[0143] Traverse all shell elements to obtain the corrected internal pressure field of the explosive , and construct a "damage level correction distribution map" by comparing it with the corrected first assessment results.

[0144] In this embodiment, the first correction unit normalizes the energy release rate of the i-th shell element with the reference rate of the standard particle size and corrects the original pressure field of the i-th shell element. This results in a corrected internal pressure field of the explosive that more accurately reflects the detonation pressure fluctuations caused by particle size refinement or coarsening, thereby improving the ability to depict the initial shock response of the explosion. The pressure assessment unit compares the corrected pressure field again with the pressure threshold Y to construct a "corrected first assessment result." Compared to the initial distribution assessment map, this eliminates the risk of misjudgment caused by particle size distribution differences, effectively avoiding problems such as high-pressure mislabeling and low-pressure missed detection, thereby improving the credibility and sensitivity of the damage level assessment.

[0145] Example 5

[0146] See also Figure 1 , the second correction module includes a second acquisition unit and an associated unit;

[0147] The second acquisition unit is used to collect the terrain geometry 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 ;

[0148] Shock wave attenuation factor of the i-th shell element , reflectivity factor and energy dissipation factor The specific way to obtain it is:

[0149] S21. The shock wave will attenuate during its propagation due to the influence of the terrain. 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 using the following formula: :

[0150] ;

[0151] in, is the distance from the ith shell element to the explosion source, Represents the attenuation index, with a value range of 1-2;

[0152] 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 by the following formula: :

[0153] ;

[0154] Where, represents the incident angle of the shock wave on the surface of the i-th shell element; Represents the reflection coefficient of different surface materials, specifically: soil ; When concrete ; Vegetation ;

[0155] 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 :

[0156] ;

[0157] Where, Indicates the basic medium absorption coefficient, including: soil ; When concrete ; Vegetation ;

[0158] Indicates the coverage density correction factor, including: when bare ground Grassland ; Vegetation .

[0159] 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 :

[0160] ;

[0161] Where, 、 and are the shock wave attenuation factors of the i-th shell element. , reflectivity factor and energy dissipation factor The weight coefficients are , and the sum of the weights is 1.

[0162] The second correction module further includes a second correction unit, a simulated plastic deformation output unit and an early warning unit;

[0163] 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 :

[0164] ;

[0165] The simulated plastic deformation output unit is used to identify the brittle fracture point area through the following steps, including:

[0166] 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 the “high stress fracture” label is added;

[0167] The fracture strength threshold X is obtained from experiments or material manuals, including AISI4340, which is set to 1300–1500 MPa;

[0168] 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 label "velocity instability fracture" is added;

[0169] 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 added with the label "thermal softening fracture"; the normalized temperature value of the i-th shell element = 1, the stress returns to zero, but softening has already begun, so most materials =0.75, it enters the thermal softening sensitive zone, which is the key window before dynamic instability;

[0170] S14. When any of the conditions S11-S13 are met, the i-th shell element is marked as a “brittle fracture point region” and the information is recorded in the database, including: element number i, timestamp, spatial position coordinates and corresponding label;

[0171] S15. Count the total number of brittle fracture point areas and generate corresponding damage levels, including:

[0172] 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 level 1 damage level;

[0173] When 5≤the total number of brittle fracture points in the three-dimensional internal structure model of the target object is less than 15, it indicates local brittle fracture, local strength reduction of the structure, and the presence of localized clustered areas of brittle fracture points, generating a secondary damage level;

[0174] When the total number of brittle fracture points in the target's 3D internal structure model is ≥15, it indicates multiple concentrated brittle fractures, and the main structure is damaged and concentrated on the main beam or supporting surface, generating a third level of damage level;

[0175] 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.

[0176] In this embodiment, the second correction module collects the terrain geometry data of the location of the shell unit in real time through the second acquisition unit, and constructs an environmental response data set including the shock wave attenuation factor, reflectivity factor and energy dissipation factor. It can accurately simulate the complex behavior of the shock wave during the propagation process in different landforms (such as mountains, hard ground, and grassland), and overcome the hypothetical error of ignoring environmental factors in traditional explosion models.

[0177] The associated units process the terrain impact data through dimensionless normalization and calculate the "environmental impact index" based on preset weights (such as attenuation-dominant terrain, reflection-dominant surface, and dissipation-dominant medium) to achieve weighted modeling of the differential effects of different types of surface conditions on explosion effects, thereby making the damage analysis sensitive to the actual terrain conditions.

[0178] Compared to the assumption of wave propagation 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 element, making the actual impact value obtained closer to the on-site measurement results. This makes it particularly suitable for evaluating explosion effects in urban environments, forested areas, or military protection areas with complex terrain.

[0179] The module incorporates a combined correction method for the base medium absorption coefficient and the cover density correction factor. This method distinguishes natural conditions such as bare land, grassland, and vegetation, and assigns independent parameters to different material types (such as soil, concrete, and vegetation). This supports the coordinated prediction modeling of explosion effects in areas where natural surfaces and artificial structures are mixed. A second correction unit applies the environmental impact index to the explosion result of the i-th shell element, automatically correcting its damage level. This improves the smoothness and spatial resolution of the damage level distribution map at terrain boundaries, avoiding damage level errors or prediction bias in areas with drastic terrain changes.

[0180] Example 6

[0181] Please refer to Figure 1 The warning unit is configured to generate a corresponding warning instruction based on the corrected first assessment result and the corresponding damage level, including:

[0182] A damage threshold L is preset, and the total number Nn of shell elements in the i-th high-pressure area is counted. When Nn < L, it means that the damage is within the tolerance range, and continuous monitoring is performed without triggering an early warning, and only the status is recorded.

[0183] When Nn ≥ L and the damage level is level 1, it means that the overall damage to the target is still classified as "minor brittle fracture" and there is a risk of local strength reduction. The first warning instruction is generated, including: reducing the explosive filling rate by 5% to 10%, extending the detonation delay by 5% to 10%, and keeping the detonation point unchanged or slightly adjusting it to 1-2 meters away from the edge of the high-pressure zone;

[0184] When Nn ≥ L and the damage level is level 2, indicating a risk of local intensity reduction, the risk trend is increasing compared to the first warning instruction. A second warning instruction is generated, including: reducing the explosive filling rate by 11% to 15%, extending the detonation delay by 11% to 15%, and keeping the detonation point unchanged or slightly adjusting it to 3-5 meters away from the edge of the high-pressure area.

[0185] When Nn ≥ L and the damage level is level 3, indicating a significant "local intensity reduction" damage trend, a third warning instruction is generated, including: reducing the explosive filling rate by 16% to 20%, extending the detonation delay by 16% to 20%, and keeping the detonation point unchanged or slightly adjusting it to 6-8 meters away from the edge of the high-pressure area;

[0186] When Nn≥L and the damage level is level four, it means that the structure has experienced "concentrated damage" or "connected crack penetration" and a fourth warning instruction is generated, including: stopping blasting in the area or reducing the explosive filling rate by more than 40%, switching to extended or segmented detonation mode, keeping the detonation point 9 meters away from the high-voltage area, and adopting a multi-point dispersed detonation scheme.

[0187] In this embodiment, the early warning unit establishes a four-level early warning stratification strategy based on the preset damage quantity threshold L and the actual number Nn of high-pressure area shell units, and combines the severity of damage levels, which can accurately reflect different states from minor damage, enhanced risk to structural failure, and realize the dynamic evolution process of "no early warning - weak early warning - medium early warning - strong early warning - severe early warning". This system not only considers the quantity characteristic (i.e., the number of high-pressure areas) that the internal pressure field of the explosive exceeds the threshold, but also superimposes the specific damage level as the basis for joint judgment, making the early warning behavior based on the cross-analysis of multiple states, enhancing the accuracy and scientific nature of the judgment of the damage evolution trend. When the number of high-pressure areas Nn < L, the system identifies it as "within the structural bearing range", automatically enters the recording and monitoring state without triggering an early 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 early warning unit not only makes a judgment on "whether to give an early warning", but also makes a fine distinction of the risk trend based on the damage level. For example, a first-level early warning represents a minor brittle fracture omen, and a fourth-level early warning represents the critical state of structural connectivity damage, realizing the two-dimensional judgment logic of "trend discrimination + level response" for risks.

[0188] Early 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.

[0189] It should be noted that: all 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 based on artificially set rules.

[0190] Essentially, the technical solution of this 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, read-only memory (ROM), random access memory (RAM), flash memory (FLASH), hard disk, or optical disc of a computer, including several instructions for causing a computer device (which can be a personal computer, server, or network device, etc.) to execute the methods of various embodiments of this invention.

[0191] The logic and / or steps represented in the flowcharts or otherwise described herein, for example, can be considered as a sequenced list of executable instructions for implementing the logical functions, and can be embodied in any computer-readable medium for use by, or in conjunction with, an instruction execution system, apparatus, or device (e.g., a computer-based system, a system including a processor, or other system that can fetch and execute instructions from an instruction execution system, apparatus, or device). For purposes of this specification, a "computer-readable medium" can be any device that can contain, store, communicate, propagate, or transport a program for use by, or in conjunction with, an instruction execution system, apparatus, or device.

[0192] It should be noted that the above embodiments are only used to illustrate the technical solutions of the present invention and are not intended to limit the present invention. Although the present invention has been described in detail with reference to the preferred embodiments, those skilled in the art should understand that the technical solutions of the present invention may be modified or replaced by equivalents without departing from the spirit and scope of the technical solutions of the present invention, which should all be included in the scope of the claims of the present invention.

[0193] It should be noted that the above embodiments are only used to illustrate the technical solutions of the present invention and are not intended to limit the present invention. Although the present invention has been described in detail with reference to the preferred embodiments, those skilled in the art should understand that the technical solutions of the present invention may be modified or replaced by equivalents without departing from the spirit and scope of the technical solutions of the present invention, which should all be included in the scope of the claims of the present invention.

Claims

1. The damage effectiveness simulation and evaluation system based on model simulation is characterized by: 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 elements, and construct the internal pressure field of the explosive of the i-th shell element in the process of simulating the 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 this with the pressure threshold Y to obtain a first assessment result, and construct an "initial damage level distribution map" based on the first assessment 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 and damage level assessment standards, as well as 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, use the finite element network construction unit 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. 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 elements, 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 element by applying the JWL state equation. : ； Where, and Represent the empirical pressure constants at the initial and later stages of the explosion, respectively; 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, Expresses the energy ratio of explosives, represents the internal energy per unit mass of the i-th 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 it with the first assessment results.

4. The damage effectiveness simulation and evaluation system based on model simulation according to claim 3, 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 action of 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, 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, characterized in that: 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 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 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 : ； Where, Indicates the average particle size The mass ratio of explosive particles to explosives satisfies ; Indicates 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; Indicates particle size 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 : ； Where, 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, characterized in that: The pressure evaluation unit is used to convert the corrected explosive internal pressure field of the i-th shell element into Re-comparing and correcting with the pressure threshold Y to obtain a corrected first assessment 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 , 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, characterized in that: The second correction module includes a second acquisition unit and an associated unit; The second acquisition unit is used to collect the terrain geometry 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 ; Shock wave attenuation factor of the i-th shell element , reflectivity factor and energy dissipation factor The specific way to obtain it is: S21. The shock wave will attenuate during its propagation due to the influence of the terrain. 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 using 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: : ； Where, represents the incident angle of the shock wave on 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 : ； Where, 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, 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 : ； Where, 、 and are the shock wave attenuation factors of the i-th shell element. , reflectivity factor and energy dissipation factor The weight coefficients are , and the sum of the weights 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 further 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 the "high stress fracture" label is added; S12, when When ≥1, it indicates that the fracture is caused by explosive high strain rate, and the unit is marked as a brittle fracture point candidate 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 added with the "thermal softening fracture" label; S14. When any of the conditions S11-S13 are 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 label; 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 level 1 damage level; When 5≤the total number of brittle fracture points in the three-dimensional internal structure model of the target object is less than 15, it indicates local brittle fracture, local strength reduction of the structure, and the presence of localized clustered areas of brittle fracture points, generating a secondary damage level; When the total number of brittle fracture points in the target's 3D internal structure model is ≥15, it indicates multiple concentrated brittle fractures, and the main structure is damaged and concentrated on the main beam or supporting surface, generating a third level of damage level; 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 configured to generate a corresponding early warning instruction based on the corrected first assessment result and the corresponding damage level, including: A damage threshold L is preset, and the total number Nn of shell elements in the i-th high-pressure area is counted. When Nn < L, it means that the damage 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 to the target is still classified as "minor brittle fracture" and there is a risk of local strength reduction, and the first warning command is generated; When Nn ≥ L and the damage level is level 2, it indicates that there is a risk of local intensity 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 "local intensity reduction" damage trend, and the third warning instruction is generated; When Nn≥L and the damage level is level four, it means that the structure has experienced "concentrated damage" or "connected crack penetration" and the fourth warning instruction is generated.

Images