Damage assessment system and method for internal explosion of building

Abstract

The invention belongs to the technical field of target damage assessment, and relates to a damage assessment system and method for internal explosion of a building. The system comprises a target modeling module, an ammunition library module, a material library module, a damage calculation module, a damage tree module and a result display module. According to the selected warhead parameter in the ammunition library and the target characteristic parameter of the target modeling module, the damage calculation module can obtain the overpressure, impulse and fragment distribution of different positions in the building, and sends the calculation result data to the damage tree module so as to obtain the damage condition of the building. And finally, the damage effect is displayed through a result display module.

Description

Technical Field

[0001] The invention belongs to the technical field of target damage assessment and relates to a damage assessment system and method for explosion inside a building. Background Art

[0002] In modern warfare, precise air strikes on enemy building targets have become the preferred option for combat. Various types of buildings, such as command and communication buildings, ammunition depots, office buildings, communication buildings, etc., are the priority targets. At present, there are three main methods for damage assessment of building target implosion: numerical simulation, experimental simulation and theoretical derivation. Among them, numerical simulation mainly uses finite element simulation software to establish a finite element model of the corresponding building, and obtains the degree of damage to the building after the explosion through simulation software; this method requires detailed modeling and simulation analysis of the building, consumes a lot of computing resources, takes a long time to analyze, and can only analyze a single explosion condition for a specific target, and its versatility is poor. Experimental simulation usually uses a scaled building model to conduct explosion experiments to evaluate the damage of the explosion to the building target, but the experimental simulation takes a long time and is costly. Theoretical derivation is usually based on empirical formulas to estimate the shock wave overpressure of the building components, so as to obtain the dynamic response of a certain component of the building. This derivation method is mainly suitable for evaluating the damage of the external explosion of the building, as well as the dynamic response and damage of a single component. At present, there is no timely and efficient damage assessment method that can accurately evaluate the propagation law of shock waves inside buildings and their damaging effects. Therefore, it is difficult to meet the urgent need for damage assessment of buildings in the event of an explosion. Summary of the invention

[0003] The technical problem solved by the present invention is: to overcome the shortcomings of the prior art and to propose a damage assessment system and method for explosions inside buildings. The method can be used to assess the probability of damage caused by weapon strikes to buildings, the number of strikes, and to optimize strike plans, thereby providing an effective platform and technical support for the protective design of buildings, the research and development of warhead power, and firepower control.

[0004] The objective of the present invention is achieved through the following technical solutions:

[0005] A damage assessment system for explosion inside a building, comprising a target modeling module, an ammunition depot module, a material depot module, a damage calculation module, a damage tree module and a result display module;

[0006] The target modeling module is used to input the building target characteristic parameters and upload the characteristic parameter data to the damage calculation module;

[0007] The ammunition depot module is used to determine the warhead parameters used to strike the building, and transmit the warhead parameters to the damage calculation module;

[0008] The material library module is used to input relevant material parameters of different components of the building, and transmit the material parameters of different components of the building to the damage calculation module;

[0009] The damage calculation module is used to receive characteristic parameters transmitted by the target modeling module, and is also used to receive warhead parameters transmitted by the ammunition depot module, and is also used to receive material parameters of different components transmitted by the material depot module, and calculates the coordinates of the warhead penetrating the building according to the received warhead parameters, characteristic parameters and material parameters of the building target, and calculates the overpressure value inside the building according to the coordinates using a CFD solver, and calculates the coordinates of the warhead fragments at the same time, and uploads the calculated coordinates of the warhead penetrating the building, the overpressure value inside the building, and the coordinates of the warhead fragments to the damage tree module and the result display module;

[0010] The damage tree module includes a component damage criterion, and the damage tree module is used to receive the coordinates of the warhead penetrating the building, the overpressure value inside the building, and the coordinates of the warhead fragments output by the damage calculation module, and calculate the damage probability of different components of the building according to the received coordinates of the warhead penetrating the building, the overpressure value inside the building, the coordinates of the warhead fragments, and the component damage criterion, to obtain the overall damage probability of the building, and upload the overall damage probability of the building to the result display module;

[0011] The result display module is used to receive the coordinates of the warhead penetrating the building, the overpressure value inside the building, and the coordinates of the warhead fragments output by the damage calculation module, and is also used to receive the overall damage probability of the building output by the damage tree module, and based on the received coordinates of the warhead penetrating the building, the overpressure value inside the building, the coordinates of the warhead fragments and the damage probability, the process of the warhead penetrating the building and the explosion scene of the warhead are displayed in a three-dimensional visual form.

[0012] Further, the target characteristic parameters include the size characteristic parameters, structural form and basic components of the building, wherein the size characteristic parameters of the building include the x, y and z coordinates of the reference point;

[0013] The structural form includes the number of rows, columns and floors of the building, and the length of the room in the x, y and z directions;

[0014] The basic components include the lengths of concrete columns in x, y and z directions, the lengths of concrete slabs in x, y and z directions, the lengths of concrete beams in x, y and z directions, the lengths of walls in x, y and z directions, and the lengths of doors in x, y and z directions;

[0015] Furthermore, the warhead parameters include geometric parameters, physical parameters and projectile-target intersection parameters of the warhead;

[0016] The geometric parameters of the warhead include the radius of the warhead, the length of the cylindrical part of the projectile, the length of the charge, and the diameter of the charge;

[0017] The physical parameters of the warhead include warhead weight, charge type, and equivalent TNT mass;

[0018] The missile-target intersection parameters include warhead velocity, starting position, CEP radius, pitch angle, yaw angle, roll angle, trajectory inclination angle and trajectory deviation angle;

[0019] Furthermore, the material parameters include parameters of concrete and parameters of metal, the parameters of concrete include axial compressive strength, axial tensile strength, elastic modulus, stirrup ratio and longitudinal reinforcement ratio; the parameters of metal include density, shear modulus, yield stress and melting temperature;

[0020] The coordinates of the warhead when it touches the building are: (x, y, z);

[0021] The coordinates of the explosion after the warhead penetrates the building are: (x, y', z); where y' = yh s ,h s is the penetration depth of the warhead;

[0022] The formula used to obtain the penetration depth of the warhead is the ACE formula;

[0023]

[0024] Where d w is the warhead diameter; D w is the warhead caliber density; f c ′ is the axial compressive strength of concrete; v 0 Warhead speed; m p is the warhead mass;

[0025] The following formula is solved by the Roe approximate Riemann solver to obtain the overpressure value Q inside the building:

[0026]

[0027] Where ρ is the fluid density, u is the velocity vector, E is the specific total energy, is the partial derivative of time, is the divergence, represents the Kronecker product operation, and γ=1.4 represents the adiabatic index of air.

[0028] Furthermore, the component damage criteria include penetration damage criteria, fragment damage criteria and shock wave damage criteria;

[0029] The penetration damage criteria include a 0-1 damage criterion and a linear damage criterion; the 0-1 damage criterion is

[0030]

[0031] The linear damage criterion is

[0032]

[0033] Where P 1 is the damage probability of 0-1 damage criterion; P 2 is the linear damage criterion damage probability; h s is the penetration depth; D is the floor thickness;

[0034] The fragment damage criteria include the limit penetration speed damage criterion and the hit number damage criterion;

[0035] The limit penetration speed damage criterion is:

[0036]

[0037] Where P 3 is the damage probability of the limit penetration speed damage criterion; V 破片速度 is the fragment velocity; V 极限穿透速度 is the ultimate penetration speed;

[0038] The damage criterion of the number of hits is

[0039]

[0040] Where P 4 N is the damage probability of the hit number damage criterion; 构件击中 is the number of components hit by fragments; N component damage is the number of fragments that hit the component.

[0041] The shock wave damage criteria include overpressure damage criteria, impulse damage criteria and overpressure-impulse damage criteria;

[0042] The overpressure damage criterion includes personnel overpressure damage criterion and component overpressure damage criterion;

[0043] The personnel overpressure injury criterion is:

[0044]

[0045] Among them, P 5 P is the injury probability of personnel overpressure injury criterion; 人员 is the overpressure value to which personnel are subjected; P 轻伤 P is the overpressure threshold when personnel are slightly injured; 中等伤P is the overpressure threshold when personnel are moderately injured; 重伤 P is the overpressure threshold when personnel are seriously injured; 危重伤害 is the overpressure threshold when personnel are seriously injured; P 死亡 It is the overpressure threshold at which people die.

[0046] The component overpressure damage criterion is

[0047]

[0048] Where P 6 is the probability of damage; P 构件 is the overpressure value of the component; P 轻微毁伤 is the overpressure threshold when the component is slightly damaged; P 中等毁伤 is the overpressure threshold when the component is moderately damaged; P 严重毁伤 It is the overpressure threshold when the component is seriously damaged.

[0049] The impulse damage criterion is

[0050]

[0051] Where P 7 I is the impulse damage criterion damage probability; 构件 I is the impulse received by the component; 构件阈值 is the impulse threshold of the component.

[0052] The overpressure impulse damage criterion is

[0053]

[0054] Where: P 0.2 , P 0.5 and P 0.8 is the overpressure asymptote of the PI curve corresponding to the damage degree; I 0.2 , I 0.5 and I 0.8 is the impulse asymptote of the PI curve corresponding to the degree of damage;

[0055]

[0056] f c ′ is the axial compressive strength of concrete; ρ s is the stirrup ratio; ρ is the longitudinal reinforcement ratio; H is the length of the column in the y direction, b is the length of the column in the x direction, and h is the length of the column in the z direction.

[0057] The overall damage probability of the building is calculated by the damage probability of each component through a logical relationship. After obtaining the damage probability of each component of the building, the damage tree can be used to analyze according to the following method to finally obtain the overall damage probability of the building.

[0058] The damage probability of the lowest level component i in the M-level damage tree is C i is the influence weight of component i, which has the following definition:

[0059] (1) If a sub-function n of the mth layer is composed of T lower-level components in parallel, that is, when all the lower-level components are destroyed, the upper level can be destroyed. At this time, its destruction probability is

[0060]

[0061] (2) If the T next-level components of a sub-function n of the mth layer are independent and connected in series, that is, if one of the next-level components is damaged, then the previous-level function is considered damaged. At this time, its damage probability is

[0062]

[0063] Therefore, according to the damage probability of the bottom layer components and the logical connection relationship of the damage tree of this layer, the damage probability of the event on the previous layer is calculated according to the above two formulas, and finally the damage probability of the entire building is obtained: 1 .

[0064] A method for assessing damage caused by explosion inside a building comprises the following steps:

[0065] The first step is to input the building target characteristic parameters through the target modeling module and upload the characteristic parameter data to the damage calculation module;

[0066] The second step is to determine the warhead parameters used to attack the building through the ammunition depot module and transmit the warhead parameters to the damage calculation module;

[0067] The third step is to input the relevant material parameters of different components of the building through the material library module, and transmit the material parameters of different components of the building to the damage calculation module;

[0068] The fourth step is to calculate the end point coordinates of the building after the warhead penetrates the building through the warhead parameters received by the damage calculation module, the characteristic parameters of the building target and the material parameters. The CFD solver is called at this end point coordinate to calculate the overpressure value inside the building, and the distribution position of the warhead fragments at this end point coordinate is calculated. The calculated damage results are uploaded to the damage tree module and the result display module.

[0069] Step 5: Calculate the damage probability of building components through the shock wave overpressure distribution and fragment distribution received by the damage tree module and the component damage criterion to obtain the overall damage probability of the building, and upload the damage probability to the result display module;

[0070] In the sixth step, the result display module displays the process of the weapon striking the building target and the explosion scene in a three-dimensional visual form based on the received overpressure data, fragmentation data and damage probability.

[0071] The beneficial effects of the present invention are:

[0072] The CFD solver can accurately calculate the overpressure distribution inside the building, evaluate the weak points of the building, estimate the amount of ammunition consumed by the building under different damage probabilities, evaluate the damage of different weapons to the building, and generate corresponding damage result reports, providing a decision-making basis for the development of my country's military damage assessment technology and the precise strike of building targets.

[0073] The present invention discloses a damage assessment system and method for explosions inside buildings, belonging to the technical field of target damage assessment. The system includes a target modeling module, an ammunition depot module, a material library module, a damage calculation module, a damage tree module and a result display module. According to the warhead parameters in the selected ammunition depot and the target characteristic parameters of the target modeling module, the damage calculation module can obtain the overpressure, impulse and fragment distribution at different positions in the building, and send the calculation result data to the damage tree module to obtain the damage of the building, and finally display the damage effect through the result display module. The present invention can estimate the ammunition consumption of the building under different damage probabilities; can compare the damage probabilities of different warheads attacking the same building target; can evaluate the weak positions of the building, generate a damage assessment report, and provide support for building target protection, weapon and ammunition development and combat command decision-making. BRIEF DESCRIPTION OF THE DRAWINGS

[0074] Figure 1 It is a schematic diagram of the connection of each module of the present invention;

[0075] Figure 2 It is a schematic diagram of CFD solver data exchange;

[0076] Figure 3 It is a flow chart of the damage assessment method for internal explosion in a building;

[0077] Figure 4 It is a schematic diagram of parameterized modeling in Example 1;

[0078] Figure 5 is a pressure cloud diagram of the building in Example 1;

[0079] Figure 6 is the building damage probability cloud map in Example 1;

[0080] Figure 7 This is a simulation diagram of building damage assessment in Example 1. DETAILED DESCRIPTION

[0081] The following will be combined with the drawings in the embodiments of the present invention to clearly and completely describe the technical solutions in the embodiments of the present invention. It is obvious that the described embodiments are only part of the embodiments of the present invention, not all of the embodiments. Based on the embodiments of the present invention, all other embodiments obtained by ordinary technicians in this field without creative work are within the scope of protection of the present invention.

[0082] like Figure 1 As shown, a damage assessment system for explosion inside a building includes a target modeling module, an ammunition depot module, a material library module, a damage calculation module, a damage tree module and a result display module;

[0083] The target modeling module is used to input the target characteristic parameters of the building and upload the characteristic parameter data to the damage calculation module; the target characteristic parameters include the size characteristic parameters, structural form and basic components of the building, wherein the size characteristic parameters of the building include the x, y and z coordinates of the reference point; the structural form includes the number of rows, columns and floors of the building, and the length of the room in the x, y and z directions; the basic components include the length of the concrete column in the x, y and z directions, the length of the concrete slab in the x, y and z directions, the length of the concrete beam in the x, y and z directions, the length of the wall in the x, y and z directions, and the length of the door in the x, y and z directions;

[0084] The ammunition depot module is used to determine the warhead parameters used to strike the building, and transmit the warhead parameters to the damage calculation module; the warhead parameters include the warhead's geometric parameters, physical parameters and projectile-target intersection parameters; the warhead's geometric parameters include the warhead radius, the length of the cylindrical part of the projectile body, the charge length, and the charge diameter; the warhead's physical parameters include the warhead weight, the charge type, and the equivalent TNT mass; the projectile-target intersection parameters include the warhead speed, starting position, CEP radius, pitch angle, yaw angle, roll angle, trajectory inclination angle, and trajectory deviation angle;

[0085] The material library module is used to input relevant material parameters of different components of the building, and transmit the material parameters of different components of the building to the damage calculation module; the material parameters include parameters of concrete and parameters of metal, the parameters of concrete include axial compressive strength, axial tensile strength, elastic modulus and fatigue deformation modulus; the parameters of metal include density, shear modulus, yield stress and melting temperature;

[0086] The damage calculation module is used to receive characteristic parameters transmitted by the target modeling module, and is also used to receive warhead parameters transmitted by the ammunition depot module, and is also used to receive material parameters of different components transmitted by the material depot module, and calculates the coordinates of the warhead penetrating the building according to the received warhead parameters, characteristic parameters and material parameters of the building target, and calculates the overpressure value inside the building according to the coordinates, and calculates the coordinates of the warhead fragments at the same time, and uploads the calculated coordinates of the warhead penetrating the building, the overpressure value inside the building and the coordinates of the warhead fragments to the damage tree module and the result display module; the coordinates when the warhead contacts the building are: (x, y, z);

[0087] The coordinates of the explosion after the warhead penetrates the building are: (x, y', z); where y' = yh, h is the penetration depth of the warhead;

[0088] The formula used to obtain the penetration depth of the warhead is the ACE formula;

[0089]

[0090] Where d w is the warhead diameter; D w is the warhead caliber density; f c ′ is the axial compressive strength of concrete; v 0 Warhead speed; m p is the warhead mass;

[0091] The overpressure value Q inside the building is calculated by the following formula:

[0092]

[0093] Where ρ is the fluid density, u is the velocity vector, E is the specific total energy, is the partial derivative of time, is the divergence, represents the Kronecker product operation, and γ=1.4 represents the adiabatic index of air.

[0094] The above formula is solved using the Roe approximate Riemann solver. During the solution process, the level set function is used to map the boundary geometry to a Cartesian grid.

[0095] The damage tree module includes a component damage criterion, which is used to receive the coordinates of the warhead penetrating the building, the overpressure value inside the building, and the coordinates of the warhead fragments output by the damage calculation module, and calculate the damage probability of different components of the building according to the received coordinates of the warhead penetrating the building, the overpressure value inside the building, the coordinates of the warhead fragments, and the component damage criterion to obtain the overall damage probability of the building, and upload the overall damage probability of the building to the result display module; further, the component damage criterion includes a penetration damage criterion, a fragment damage criterion, and a shock wave damage criterion;

[0096] The penetration damage criteria include a 0-1 damage criterion and a linear damage criterion; the 0-1 damage criterion is

[0097]

[0098] The linear damage criterion is

[0099]

[0100] Where P 1 is the damage probability of 0-1 damage criterion; P 2 is the linear damage criterion damage probability; h is the penetration depth; D is the floor thickness;

[0101] The fragment damage criteria include the limit penetration speed damage criterion and the hit number damage criterion;

[0102] The limit penetration speed damage criterion is:

[0103]

[0104] Where P 3 is the damage probability of the limit penetration speed damage criterion; V 破片速度 is the fragment velocity; V 极限穿透速度 is the ultimate penetration speed;

[0105] The damage criterion of the number of hits is

[0106]

[0107] Where P 4 N is the damage probability of the hit number damage criterion; 构件击中 is the number of components hit by fragments; N component damage is the number of fragments that hit the component.

[0108] The shock wave damage criteria include overpressure damage criteria, impulse damage criteria and overpressure-impulse damage criteria;

[0109] The overpressure damage criterion includes personnel overpressure damage criterion and component overpressure damage criterion;

[0110] The personnel overpressure injury criterion is:

[0111]

[0112] Among them, P 5 P is the injury probability of personnel overpressure injury criterion; 人员 is the overpressure value to which personnel are subjected; P 轻伤 P is the overpressure threshold when personnel are slightly injured; 中等伤 P is the overpressure threshold when personnel are moderately injured; 重伤 P is the overpressure threshold when personnel are seriously injured; 危重伤害 is the overpressure threshold when personnel are seriously injured; P 死亡 It is the overpressure threshold at which people die.

[0113] The component overpressure damage criterion is

[0114]

[0115] Where P 6 is the probability of damage; P 构件 is the overpressure value of the component; P 轻微毁伤 is the overpressure threshold when the component is slightly damaged; P 中等毁伤 is the overpressure threshold when the component is moderately damaged; P 严重毁伤 It is the overpressure threshold when the component is seriously damaged.

[0116] The impulse damage criterion is

[0117]

[0118] Where P 7 I is the impulse damage criterion damage probability; 构件 I is the impulse received by the component; 构件阈值 is the impulse threshold of the component.

[0119] The overpressure impulse damage criterion is

[0120]

[0121] Where: P 0.2 , P 0.5 and P 0.8 is the overpressure asymptote of the PI curve corresponding to the damage degree; I 0.2 , I 0.5 and I 0.8 is the impulse asymptote of the PI curve corresponding to the degree of damage;

[0122]

[0123]

[0124] The overall damage probability of the building is calculated by the damage probability of each component through a logical relationship. After obtaining the damage probability of each component of the building, the damage tree can be used to analyze according to the following method to finally obtain the overall damage probability of the building.

[0125] The overall damage probability of the building is calculated by the damage probability of each component through a logical relationship. After obtaining the damage probability of each component of the building, the damage tree can be used to analyze according to the following method to finally obtain the overall damage probability of the building.

[0126] The damage probability of the lowest level component i in the M-level damage tree is C i is the influence weight of component i, which has the following definition:

[0127] (1) If a sub-function n of the mth layer is composed of T lower-level components in parallel, that is, when all the lower-level components are destroyed, the upper level can be destroyed. At this time, its destruction probability is

[0128]

[0129] (2) If the T next-level components of a sub-function n of the mth layer are independent and connected in series, that is, if one of the next-level components is damaged, then the previous-level function is considered damaged. At this time, its damage probability is

[0130]

[0131] Therefore, according to the damage probability of the bottom layer components and the logical connection relationship of the damage tree of this layer, the damage probability of the event on the previous layer is calculated according to the above two formulas, and finally the damage probability of the entire building is obtained: 1 .

[0132] The result display module is used to receive the coordinates of the warhead penetrating the building, the overpressure value inside the building, and the coordinates of the warhead fragments output by the damage calculation module, and is also used to receive the overall damage probability of the building output by the damage tree module, and based on the received coordinates of the warhead penetrating the building, the overpressure value inside the building, the coordinates of the warhead fragments and the damage probability, the process of the warhead penetrating the building and the explosion scene of the warhead are displayed in a three-dimensional visual form.

[0133] Example

[0134] There is a warhead a, the target is a building on the ground. This system conducts damage assessment research on the warhead to the target, providing certain guidance and assistance to both the attacker and the defender.

[0135] like Figure 3As shown, damage assessment for an internal explosion in a building requires the following steps:

[0136] Step 1. In the target modeling module of the internal explosion damage assessment system in the building, input the parameters related to the target to determine the structure of the building and the materials of the corresponding components; the detailed parameters of the headquarters building are: a 6-story 8-span frame structure building; the x coordinate of the reference point is 0.25; the y coordinate of the reference point is 2.5; the z coordinate of the reference point is 0; the row is 9; the room length x is 5; the column is 4; the room length y is 5; the layer is 4; the room length z is 5; the concrete column x is 0.5; y is 5; z is 0.5; the concrete slab x is 4.5; y is 0.2; z is 4.5; the concrete beam dx is 4.5; y is 0.2; z is 0.5; the wall dx is 4.5; y is 4.8; z is 0.2; the door dx is 4.5; y is 1.8; z is 0.2, Figure 4 Schematic diagram for parametric modeling.

[0137] Step 2: In the warhead module of the internal explosion damage assessment system in the building, select the corresponding warhead and fill in the missile-target intersection parameters; the specific parameters are shown in the following table:

[0138] Table 1: Warhead parameters

[0139]

[0140] Table 2: Projectile-target intersection parameters

[0141]

[0142] Step 3: Select material parameters for different components; select C30 concrete for concrete beams and columns; select C20 concrete for concrete plates; and select C15 concrete for concrete slabs.

[0143] Step 4: The explosion position of the warhead in the building is obtained through penetration calculation, the distribution of shock wave overpressure and the loads on building components, equipment and personnel are calculated through the CFD solver, and the damage of the projectile fragments to the building components, equipment and personnel is determined through fragment calculation; wherein, the penetration effect algorithm mainly includes the ACE formula widely used in the current penetration effect calculation and analysis, which will not be repeated here; in this embodiment, the penetration route of the warhead is calculated through the penetration effect algorithm, and the specific position of the final explosion point is given in combination with the delayed explosion time set by the missile fuze, wherein the final explosion point of the warhead a is calculated to be located in the corner room on the 4th floor;

[0144] Step 5: Select the corresponding damage criterion and damage criteria according to the damage element type to carry out damage analysis of building explosion, among which the damage criterion for personnel is shown in Table 3:

[0145] Table 3: Criteria for personnel damage

[0146]

[0147] The damage criteria for building components are shown in Table 4:

[0148] Table 4: Criteria for damage to building components

[0149]

[0150] According to the results calculated by CFD solver and combined with Table 3 and Table 4, single component damage analysis of personnel and building components is carried out; the pressure cloud map of the building is as follows Figure 5 As shown in the figure; the damage of the entire building is analyzed through the damage tree module, and the damage effect of the building is finally obtained; the damage probability cloud map of the building is shown in the figure Figure 6 shown.

[0151] Step 6: Visualize the fragment dispersion process, shock wave overpressure distribution, damage level and damage probability through a three-dimensional visualization platform.

[0152] The simulation demonstration screen is as follows Figure 7 shown.

[0153] In summary, the above are only preferred embodiments of the present invention and are not intended to limit the protection scope of the present invention. Any modifications, equivalent substitutions, improvements, etc. made within the spirit and principles of the present invention should be included in the protection scope of the present invention.

Claims

1. A damage assessment system for explosions inside buildings, characterized in that: The damage assessment system includes a target modeling module, an ammunition depot module, a material depot module, a damage calculation module and a damage tree module; The target modeling module is used to input the building target characteristic parameters and upload the characteristic parameter data to the damage calculation module; The ammunition depot module is used to determine the warhead parameters used to strike the building, and transmit the warhead parameters to the damage calculation module; The material library module is used to input relevant material parameters of different components of the building, and transmit the material parameters of different components of the building to the damage calculation module; The damage calculation module is used to receive characteristic parameters transmitted by the target modeling module, and is also used to receive warhead parameters transmitted by the ammunition depot module, and is also used to receive material parameters of different components transmitted by the material depot module, and calculates the coordinates of the warhead penetrating the building according to the received warhead parameters, characteristic parameters and material parameters of the building target, and calculates the overpressure value inside the building according to the coordinates using a CFD solver, and calculates the coordinates of the warhead fragments at the same time, and uploads the calculated coordinates of the warhead penetrating the building, the overpressure value inside the building, and the coordinates of the warhead fragments to the damage tree module; The damage tree module includes a component damage criterion, and the damage tree module is used to receive the coordinates of the warhead penetrating the building, the overpressure value inside the building, and the coordinates of the warhead fragments output by the damage calculation module, and calculate the damage probability of different components of the building based on the received coordinates of the warhead penetrating the building, the overpressure value inside the building, the coordinates of the warhead fragments, and the component damage criterion to obtain the overall damage probability of the building.

2. A damage assessment system for explosion inside a building according to claim 1, characterized in that: The damage assessment system also includes a result display module, which is used to receive the coordinates of the warhead penetrating the building, the overpressure value inside the building, and the coordinates of the warhead fragments output by the damage calculation module, and is also used to receive the overall damage probability of the building output by the damage tree module, and based on the received coordinates of the warhead penetrating the building, the overpressure value inside the building, the coordinates of the warhead fragments and the damage probability, the process of the warhead penetrating the building and the explosion scene of the warhead are displayed in a three-dimensional visual form.

3. A damage assessment system for explosion inside a building according to claim 1 or 2, characterized in that: The target characteristic parameters include the size characteristic parameters, structural form and basic components of the building, wherein the size characteristic parameters of the building include the x, y and z coordinates of the reference point; The structural form includes the number of rows, columns and floors of the building, and the length of the room in the x, y and z directions; The basic components include the lengths of concrete columns in x, y and z directions, the lengths of concrete slabs in x, y and z directions, the lengths of concrete beams in x, y and z directions, the lengths of walls in x, y and z directions, and the lengths of doors in x, y and z directions.

4. The damage assessment system for explosion inside a building according to claim 1, characterized in that: The parameters of the warhead include geometric parameters, physical parameters and projectile-target intersection parameters of the warhead; The geometric parameters of the warhead include the radius of the warhead, the length of the cylindrical part of the projectile, the length of the charge and the diameter of the charge; The physical parameters of the warhead include warhead weight, charge type and equivalent TNT mass; The missile-target intersection parameters include warhead velocity, starting position, CEP radius, pitch angle, yaw angle, roll angle, trajectory inclination angle and trajectory deviation angle.

5. The damage assessment system for explosion inside a building according to claim 1, characterized in that: The material parameters include parameters of concrete and parameters of metal. The parameters of concrete include axial compressive strength, axial tensile strength, elastic modulus, stirrup ratio and longitudinal reinforcement ratio; the parameters of metal include density, shear modulus, yield stress and melting temperature.

6. The damage assessment system for explosion inside a building according to claim 1, characterized in that: The coordinates of the warhead penetrating the building include: The coordinates of the warhead when it touches the building are: (x, y, z); The coordinates of the explosion after the warhead penetrates the building are: (x, y', z); where y' = yh s ,h s is the penetration depth of the warhead; The formula used to obtain the penetration depth of the warhead is the ACE formula; Where d w is the warhead diameter; D w is the warhead caliber density; f c ′ is the axial compressive strength of concrete; v0 is the impact velocity of the warhead; m p is the warhead mass; The following formula is solved by the Roe approximate Riemann solver to obtain the overpressure value Q inside the building: Where ρ is the fluid density, u is the velocity vector, E is the specific total energy, is the partial derivative of time, is the divergence, represents the Kronecker product operation, and γ=1.4 represents the adiabatic index of air.

7. The damage assessment system for explosion inside a building according to claim 1, characterized in that: The component damage criteria include penetration damage criteria, fragment damage criteria and shock wave damage criteria; The penetration damage criteria include 0-1 damage criterion and linear damage criterion; The 0-1 damage criterion is The linear damage criterion is Where P1 is the damage probability of 0-1 damage criterion; P2 is the damage probability of linear damage criterion; h s is the penetration depth; D is the floor thickness; The fragment damage criteria include the limit penetration speed damage criterion and the hit number damage criterion; The limit penetration speed damage criterion is: Where P3 is the damage probability of the limit penetration speed damage criterion; V 破片速度 is the fragment velocity; V 极限穿透速度 is the ultimate penetration speed; The damage criterion of the number of hits is Where P4 is the damage probability of the number of hits damage criterion; N 构件击中 is the number of components hit by fragments; N component damage is the number of fragments that hit the component.

8. The damage assessment system for explosion inside a building according to claim 1, characterized in that: The shock wave damage criteria include overpressure damage criteria, impulse damage criteria and overpressure-impulse damage criteria; The overpressure damage criterion includes personnel overpressure damage criterion and component overpressure damage criterion; The personnel overpressure injury criterion is: Among them, P5 is the injury probability of personnel overpressure injury criterion; P 人员 is the overpressure value to which personnel are subjected; P 轻伤 P is the overpressure threshold when personnel are slightly injured; 中等伤 P is the overpressure threshold when personnel are moderately injured; 重伤 is the overpressure threshold when personnel are seriously injured; P 危重伤害 is the overpressure threshold when personnel are seriously injured; P 死亡 is the overpressure threshold at which people die.

9. The damage assessment system for explosion inside a building according to claim 1, characterized in that: The component overpressure damage criterion is Where P6 is the probability of damage; P 构件 is the overpressure value of the component; P 轻微毁伤 is the overpressure threshold when the component is slightly damaged; P 中等毁伤 is the overpressure threshold when the component is moderately damaged; P 严重毁伤 is the overpressure threshold when the component is seriously damaged; The impulse damage criterion is Where P7 is the impulse damage criterion damage probability; I 构件 I is the impulse received by the component; 构件阈值 is the impulse threshold of the component; The overpressure impulse damage criterion is Where: P 0.2 , P 0.5 and P 0.8 is the overpressure asymptote of the PI curve corresponding to the damage degree; I 0.2 , I 0.5 and I 0.8 is the impulse asymptote of the PI curve corresponding to the degree of damage; Among them, f c ′ is the axial compressive strength of concrete; ρ s is the stirrup ratio; ρ is the longitudinal reinforcement ratio; H is the length of the column in the y direction, b is the length of the column in the x direction, and h is the length of the column in the z direction; The damage probability of the lowest level component i in the M-level damage tree is C i is the influence weight of component i, If a sub-function n of the mth layer is composed of T lower-level components in parallel, that is, when all the lower-level components are damaged, the upper level can be damaged. At this time, its damage probability is If the T next-level components of a sub-function n of the mth layer are independent and connected in series, that is, if one of the next-level components is damaged, then the previous level function is considered damaged. At this time, its damage probability is 10. A method for assessing damage caused by explosion inside a building, characterized in that The steps of the method include: The first step is to input the building target characteristic parameters through the target modeling module and upload the characteristic parameter data to the damage calculation module; The second step is to determine the warhead parameters used to attack the building through the ammunition depot module and transmit the warhead parameters to the damage calculation module; The third step is to input the relevant material parameters of different components of the building through the material library module, and transmit the material parameters of different components of the building to the damage calculation module; The fourth step is to calculate the end point coordinates of the building after the warhead penetrates the building through the warhead parameters received by the damage calculation module, the characteristic parameters of the building target and the material parameters. The CFD solver is called at this end point coordinate to calculate the overpressure value inside the building, and the distribution position of the warhead fragments at this end point coordinate is calculated. The calculated damage results are uploaded to the damage tree module and the result display module. Step 5: Calculate the damage probability of building components through the shock wave overpressure distribution and fragment distribution received by the damage tree module and the component damage criterion to obtain the overall damage probability of the building, and upload the damage probability to the result display module; In the sixth step, the result display module displays the process of the weapon striking the building target and the explosion scene in a three-dimensional visual form based on the received overpressure data, fragmentation data and damage probability.

Images