A missile weapon damage simulation evaluation method

Abstract

The present invention discloses a missile weapon damage simulation evaluation method, which specifically relates to the technical field of weapon equipment development. The damage process calculation and analysis module selects calculation targets and configures simulation parameters in turn according to the working condition data and target type obtained by the blasting warhead parameter input module, and establishes a numerical model based on the fluid, metal, and concrete constitutive equations; the numerical model is used for calculation, and the damage process calculation and analysis module retrieves and organizes the calculation results of the user-specified working condition through a program calling the simulation result database, and then visualizes the damage data of the blasting warhead on the target through the damage effect visualization module until the calculation is completed. The present invention obtains a large amount of damage effect data at a low cost through numerical simulation, effectively makes up for the lack of small sample test data, analyzes and evaluates the damage of the blasting warhead to typical targets, and visualizes the killing effect of the blasting warhead on typical targets.

Description

Technical Field

[0001] The invention relates to the technical field of weapon equipment development, and in particular to a missile weapon damage simulation evaluation method. Background Art

[0002] In recent years, with the rapid development of missile-type high-tech weapons and equipment, the design technology and level of explosive warheads have been rapidly improved, a large number of new explosive warheads have been developed, and tactical performance has reached a new height. However, due to the inadequate analysis of the damage effect of the target during the design of the explosive warhead, the lack of attention to the damage efficiency of the explosive warhead on the target, and the lack of in-depth research on the vulnerability of the target, damage criteria and damage effects, the design of the damage element of the explosive warhead is somewhat blind, and the actual combat damage performance of the explosive warhead is poor, which is reflected in the poor performance of the weapon in some military exercises and flight target tests in non-predetermined environments. Summary of the invention

[0003] To this end, the present invention provides a missile weapon damage simulation evaluation method to solve the problems raised in the background technology.

[0004] In order to achieve the above object, the present invention provides the following technical solution: a missile weapon damage simulation evaluation method, wherein a blasting warhead parameter input module reads relevant parameters input by a user through an input interface, and reads a typical target selected by the user, and the method comprises the following steps:

[0005] Step 1: The damage process calculation and analysis module selects calculation targets and configures simulation parameters in turn according to the working condition data and target type obtained by the blasting warhead parameter input module, and establishes a numerical model based on the constitutive equations of fluid, metal, and concrete;

[0006] Step 2: The user chooses whether to use the numerical model. If the numerical model is used for calculation, the damage process calculation and analysis module calls the simulation result database through the program to retrieve and organize the calculation results of the user-specified working conditions, and then visualizes the damage data of the blasting warhead on the target through the damage effect visualization module until the calculation is completed. If it is not used, the relevant parameters can be re-entered.

[0007] Preferably, the relevant parameters input by the user include the terminal velocity of the blasting warhead, the incident angle, the charge equivalent of the blasting warhead, the filling coefficient of the blasting warhead and the distance to the target point.

[0008] Preferably, typical targets selected by the user include ammunition depots, fire shelters, beachheads, and bridges.

[0009] Preferably, the visualization processing includes outputting dynamic cloud maps, data curves and target damage degree analysis.

[0010] Preferably, in step 2, if the numerical model is not used for calculation, return to step 2 to configure the simulation parameters.

[0011] Preferably, in step 1, the fluid constitutive equation constructed is: the pressure p' of the explosion product is determined by the product of the pressure given by the state equation and the combustion fraction F, that is,

[0012] p'=F·p EOS ;

[0013] The combustion fraction F controls the chemical energy release of the explosive.

[0014] F = max(F1, F2);

[0015]

[0016] F2 = β(1-V);

[0017]

[0018] Where: t L is the detonation time of the explosive particle, which is equal to the distance between the particle and the detonation point divided by the detonation velocity D; h is the characteristic size of the particle, P CJ is the Chapman-Jouguet pressure; V CJ is the Chapman-Jouguet relative volume; F1 is used to simulate shock detonation, and F2 is used to simulate the propagation of detonation waves;

[0019] The JWL equation is used to describe the relationship between the pressure, volume and energy of detonation products, and its expression is:

[0020]

[0021] Where: A, B, R1, R2 and ω are material parameters; E0 is the initial internal energy of the material per unit volume.

[0022] Preferably, the metal constitutive equation is constructed as follows: Considering stress hardening, the influence of strain rate and thermal effect, the yield stress of the metal material is expressed as:

[0023]

[0024] Where: ε is the equivalent plastic strain; is the dimensionless plastic strain rate ( T is the reference strain rate); *is the dimensionless temperature; A, B, n, C, m are material constants, which are obtained through torsion tests at different strain rates, Hopkinson bar tests at different temperatures, and quasi-static tensile tests, or are determined by the simpler and more economical Taylor bar collision test. For specific values, please refer to relevant literature;

[0025] The failure strain considering the damage is expressed as:

[0026]

[0027] Where: * is the average stress; D1, D2, D3, D4, D5 are the input damage constants;

[0028]

[0029] When the damage D reaches 1, the material is destroyed; after the material is destroyed, the deviatoric stress will be set to zero, and the material cannot be stretched, that is, if the pressure is less than zero, it will be set to zero.

[0030] Preferably, the constitutive equation of concrete is constructed as follows: in the elastic stage, the tension and compression of the material are linear elastic; in the plastic stage, the numerical model uses the non-associated plastic flow law; the plastic potential G takes the form of a Drucker-Prager hyperbolic function:

[0031]

[0032] in, is the expansion angle, ∈ is the eccentricity of the flow potential, σ t0 is the tensile failure stress;

[0033] The yield function when expressed using effective stress is:

[0034]

[0035] in,

[0036]

[0037]

[0038]

[0039] in, is the maximum principal stress effective stress, σ b0 / σ c0 is the ratio of the biaxial compressive yield stress to the uniaxial compressive yield stress, K c is the constant stress ratio, and are the effective tensile cohesive stress and the effective compressive cohesive stress, respectively.

[0040] Preferably, the damage process calculation and analysis module adopts an explicit dynamic finite element method, which does not form an overall stiffness matrix, avoids matrix inversion, and reduces the amount of calculation. The explicit dynamic finite element method adopts central difference time integration, and the time step adopted is less than the critical value;

[0041] Solution steps: establish integral equation; divide unit; determine unit shape function; unit analysis; overall synthesis; boundary condition processing; solve equation group; output calculation results.

[0042] Preferably, a non-uniform grid method is used in modeling the damage process calculation and analysis module. For areas where physical quantities change dramatically or local areas of detailed parts of the structure, centimeter-level grids are used; meter-level grids are used in distant areas, and transition grids are used in between.

[0043] Preferably, based on the numerical calculation results of multiple working conditions, a parametric evaluation of the destructive effectiveness of the blasting warhead on different targets is achieved, and the simulation calculation process is visualized, and different damage levels are given to the damaged area of ​​the target in the range of 0 to 1 (0 is no damage, 1 is complete damage). The calculation results are analyzed and a description of the destructive effect evaluation is given, and the conditions for complete destruction of the target are given, such as parameters such as warhead speed and explosive equivalent).

[0044] The present invention has the following advantages:

[0045] The present invention obtains a large amount of damage effect data at low cost through numerical simulation, effectively makes up for the shortage of small sample test data, analyzes and evaluates the damage of blasting warheads to typical targets, and visualizes the killing effect of blasting warheads on typical targets. BRIEF DESCRIPTION OF THE DRAWINGS

[0046] Figure 1 A diagram showing the relationship between components of the damage assessment system provided by the present invention;

[0047] Figure 2 A simulation flow chart of blasting warhead damage assessment provided by the present invention;

[0048] Figure 3 An example diagram of concrete material structure damage provided by the present invention;

[0049] Figure 4 This is an example diagram of structural damage to metal materials provided by the present invention. DETAILED DESCRIPTION

[0050] The following is a description of the implementation of the present invention by specific embodiments. People familiar with the art can easily understand other advantages and effects of the present invention from the contents disclosed in this specification. Obviously, the described embodiments are 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.

[0051] Example 1

[0052] This embodiment discloses a missile weapon damage simulation evaluation method, wherein the blasting warhead parameter input module reads the relevant parameters input by the user through the input interface, and reads the typical targets selected by the user, the typical targets selected by the user include ammunition depots, fire shelters, beachheads and bridges; the relevant parameters input by the user include the terminal velocity of the blasting warhead, the incident angle, the charge equivalent of the blasting warhead, the filling coefficient of the blasting warhead and the distance to the target point;

[0053] The method comprises the following steps:

[0054] Step 1: The damage process calculation and analysis module selects calculation targets and configures simulation parameters in turn according to the working condition data and target type obtained by the blasting warhead parameter input module, and establishes a numerical model based on the constitutive equations of fluid, metal, and concrete;

[0055] In the process of processing fluid materials, it is necessary to use both constitutive equations and state equations to describe the characteristics of a material: use the constitutive equation to describe the deviatoric stress, and use the state equation to describe the relationship between pressure and volume deformation.

[0056] Taking explosives as an example, after high-energy explosives are detonated, a detonation wave will be formed inside the explosive that propagates stably at the speed of detonation. This is a strong shock wave followed by a high-speed chemical reaction zone.

[0057] The fluid constitutive equation constructed is: The pressure p' of the explosion product is determined by the product of the pressure given by the state equation and the combustion fraction F, that is,

[0058] p'=F·p EOS ;

[0059] The combustion fraction F controls the chemical energy release of the explosive.

[0060] F = max(F1, F2);

[0061]

[0062] F2 = β(1-V);

[0063]

[0064] Where: t L is the detonation time of the explosive particle, which is equal to the distance between the particle and the detonation point divided by the detonation velocity D; h is the characteristic size of the particle, P CJ is the Chapman-Jouguet pressure; V CJ is the Chapman-Jouguet relative volume; F1 is used to simulate impact initiation, and F2 is used to simulate the propagation of detonation waves; after the explosive is detonated, the combustion fraction of the explosive particles generally takes several time steps to reach 1, so this processing method expands the discontinuous detonation wave front into a narrow area, smoothing it into a rapidly changing but continuous wave front, effectively avoiding the numerical oscillation introduced by the discontinuity.

[0065] The JWL equation is used to describe the relationship between the pressure, volume and energy of detonation products, and its expression is:

[0066]

[0067] Where: A, B, R1, R2 and ω are material parameters; E0 is the initial internal energy of the material per unit volume.

[0068] The constitutive relationship of metal materials can be selected from the elastic-plastic model commonly used in the field of explosion mechanics, that is, the state equation that describes the relationship between pressure and volume deformation. The stress deviator and shape distortion are related by the elastic-plastic relationship, which includes the elastic distortion law, the plastic yield criterion and the flow law of plastic strain.

[0069] Preferably, the metal constitutive equation is constructed as follows: Considering stress hardening, the influence of strain rate and thermal effect, the yield stress of the metal material is expressed as:

[0070]

[0071] Where: ε is the equivalent plastic strain; is the dimensionless plastic strain rate ( T is the reference strain rate); * is the dimensionless temperature; A, B, n, C, m are material constants, which are obtained through torsion tests at different strain rates, Hopkinson bar tests at different temperatures, and quasi-static tensile tests, or are determined by the simpler and more economical Taylor bar collision test. For specific values, please refer to relevant literature;

[0072] The failure strain considering damage can be expressed as:

[0073]

[0074] Where: * is the average stress; D1, D2, D3, D4, D5 are the input damage constants;

[0075]

[0076] When the damage D reaches 1, the material is destroyed; after the material is destroyed, the deviatoric stress will be set to zero, and the material cannot be stretched, that is, if the pressure is less than zero, it will be set to zero.

[0077] Concrete is a mixture of cement (usually silicate cement) and aggregate. Any significant tension and bending will cause the microstructure of the concrete to crack and separate, leading to structural damage. For concrete with complex mechanical properties, the finite element analysis mainly uses the plastic damage model (Concrete Damaged Plasticity), which can be used in the explicit analysis module.

[0078] The concrete damage plasticity model (CDP) is determined based on the damage plasticity model proposed by Lubliner, Lee and Fenves. Its purpose is to provide a universal material model for analyzing the mechanical response of concrete structures under cyclic loading and dynamic loading conditions. It is a continuous medium damage model based on plasticity. The material hardens first and then softens after compression yielding, and softens after tension yielding. Compared with the diffuse cracking model and the brittle cracking model, the damage factors used for tension and compression damage are different, taking into account the differences in the tensile and compressive properties of the material. At the same time, it can be applied to unidirectional loading, cyclic loading, dynamic loading and other occasions. It assumes that the two failure mechanisms of concrete are: tensile cracking and compressive crushing. The evolution of the yield surface and the failure surface is controlled by two hardening parameters: and is the tensile effective plastic strain, is the compressive effective plastic strain. Damage plasticity needs to define three aspects, namely plastic behavior, tensile damage behavior and compressive damage behavior.

[0079] Preferably, the constitutive equation of concrete is constructed as follows: in the elastic stage, the tension and compression of the material are linear elastic; in the plastic stage, the numerical model uses the non-associated plastic flow law; the plastic potential G takes the form of a Drucker-Prager hyperbolic function:

[0080]

[0081] in, is the expansion angle, ∈ is the eccentricity of the flow potential, σ t0 is the tensile failure stress;

[0082] The yield function when expressed using effective stress is:

[0083]

[0084] in,

[0085]

[0086]

[0087]

[0088] in, is the maximum principal stress effective stress, σ b0 / σ c0 is the ratio of the biaxial compressive yield stress to the uniaxial compressive yield stress, K c is the constant stress ratio, and are the effective tensile cohesive stress and the effective compressive cohesive stress, respectively.

[0089] Step 2: The user chooses whether to use the numerical model. If the numerical model is used for calculation, the damage process calculation and analysis module calls the simulation result database through the program to retrieve and organize the calculation results of the user-specified working conditions (read the simulation result data from the database, and process the working condition data into a unified format that is easy to analyze). Then, the damage data of the blasting warhead on the target is visualized through the damage effect visualization module. The visualization includes outputting dynamic cloud maps, data curves, and target damage degree analysis until the calculation is completed. If it is not used, the relevant parameters can be re-entered. Preferably, in step 2, if the numerical model is not used for calculation, return to step 2 to configure the simulation parameters.

[0090] Among them, the calculation boundary processing method is to impose pressure outflow boundary conditions on the calculation boundary of the model to avoid the reflection of shock waves on the boundary and affect the calculation results.

[0091] Preferably, during data simulation calculations, the simulation data is used to train a deep learning model, and then the trained model is applied to the simulation result database.

[0092] Preferably, the damage process calculation and analysis module adopts an explicit dynamic finite element method, the explicit dynamic finite element method adopts central difference time integration, and the time step adopted is smaller than a critical value;

[0093] Solution steps: establish integral equation; divide unit; determine unit shape function; unit analysis; overall synthesis; boundary condition processing; solve equation group; output calculation results.

[0094] In numerical simulation, the density and division method of the grid directly affect the accuracy of the calculation and whether the calculation can be carried out smoothly. The overall density of the grid division affects the size of the model. In three-dimensional problems, the influence of the grid size on the calculation scale is more prominent. If uniform grids are used for all, the total number of grid cells will increase greatly, and the time required for numerical simulation will increase dramatically.

[0095] Preferably, a non-uniform grid method is used in modeling the damage process calculation and analysis module. A finer grid is used for areas where physical quantities change dramatically or for local areas of detailed parts of the structure; a coarser grid is used in farther areas, with a transition grid in between.

[0096] Preferably, based on the numerical calculation results of multiple working conditions, a parametric evaluation of the destructive effectiveness of the blasting warhead on different targets is achieved, and the simulation calculation process is visualized, and different damage levels are given to the damaged area of ​​the target in the range of 0 to 1 (0 is no damage, 1 is complete damage). The calculation results are analyzed and a description of the destructive effect evaluation is given, and the conditions for complete destruction of the target are given, such as parameters such as warhead speed and explosive equivalent).

[0097] Although the present invention has been described in detail above by general description and specific embodiments, it is obvious to those skilled in the art that some modifications or improvements can be made to the present invention. Therefore, these modifications or improvements made without departing from the spirit of the present invention all belong to the scope of protection claimed by the present invention.

Claims

1. A missile weapon damage simulation evaluation method, wherein the blasting warhead parameter input module reads the relevant parameters input by the user through the input interface and reads the typical target selected by the user, characterized in that: The method comprises the following steps: Step 1: The damage process calculation and analysis module selects calculation targets and configures simulation parameters in turn according to the working condition data and target type obtained by the blasting warhead parameter input module, and establishes a numerical model based on the constitutive equations of fluid, metal, and concrete; The fluid constitutive equation constructed is: The pressure p' of the explosion product is determined by the product of the pressure given by the state equation and the combustion fraction F, that is, p'=F·p EOS ; The combustion fraction F controls the chemical energy release of the explosive. F = max(F1, F2); F2 = β(1-V); Where: t L is the detonation time of the explosive particle, which is equal to the distance between the particle and the detonation point divided by the detonation velocity D; h is the characteristic size of the particle, P CJ is the Chapman-Jouguet pressure; V CJ is the Chapman-Jouguet relative volume; F1 is used to simulate shock detonation, and F2 is used to simulate the propagation of detonation waves; The JWL equation is used to describe the relationship between the pressure, volume and energy of detonation products, and its expression is: Where: A, B, R1, R2 and ω are material parameters; E0 is the initial internal energy of the material per unit volume; The metal constitutive equation constructed is: Considering stress hardening, the influence of strain rate and thermal effect, the yield stress of metal material is expressed as: Where: ε is the equivalent plastic strain; is the dimensionless plastic strain rate; T * is the dimensionless temperature; A, B, n, C, m are material constants, obtained by torsion tests at different strain rates, Hopkinson bar tests at different temperatures, and quasi-static tensile tests, or determined by simple and economical Taylor bar impact tests; The failure strain considering the damage is expressed as Where: * is the average stress; D1, D2, D3, D4, D5 are the input damage constants; When the damage D reaches 1, the material is destroyed; after the material is destroyed, the deviatoric stress will be set to zero, and the material cannot be stretched, that is, if the pressure is less than zero, it will be set to zero; The constructed constitutive equation of concrete is as follows: in the elastic stage, the tension and compression of the material are linear elastic; in the plastic stage, the numerical model uses the non-associated plastic flow law; the plastic potential G takes the form of a Drucker-Prager hyperbolic function: in, is the expansion angle, ∈ is the eccentricity of the flow potential, σ t0 is the tensile failure stress; The yield function when expressed using effective stress is: in, in, is the maximum principal stress effective stress, σ b0 / σ c0 is the ratio of the biaxial compressive yield stress to the uniaxial compressive yield stress, K c is the constant stress ratio, and are the effective tensile cohesive stress and the effective compressive cohesive stress respectively; Step 2: The numerical model is used for calculation. The damage process calculation and analysis module retrieves and organizes the calculation results of the user-specified working conditions through the simulation result database through the program, and then visualizes the damage data of the blasting warhead on the target through the damage effect visualization module until the calculation is completed; The damage process calculation and analysis module adopts the explicit dynamic finite element method, which uses central difference time integration and the time step adopted is less than the critical value. The damage process calculation and analysis module adopts the non-uniform grid method in modeling. For areas where physical quantities change drastically or local areas of detailed parts of the structure, centimeter-level grids are used; meter-level grids are used in far areas, and transition grids are used in the middle.

2. A missile weapon damage simulation evaluation method according to claim 1, characterized in that: The relevant parameters input by the user include the terminal velocity of the blasting warhead, the incident angle, the charge equivalent of the blasting warhead, the filling coefficient of the blasting warhead and the distance to the target.

3. A missile weapon damage simulation evaluation method according to claim 1, characterized in that: Typical targets selected by the user include ammunition depots, fire bunkers, beachheads, and bridges.

4. A missile weapon damage simulation evaluation method according to claim 1, characterized in that: Visualization processing includes outputting dynamic cloud maps, data curves and target damage analysis.

5. The missile weapon damage simulation evaluation method according to claim 1 is characterized in that: In step 2, if the numerical model is not used for calculation, return to step 2 to configure the simulation parameters.

6. A missile weapon damage simulation evaluation method according to claim 1, characterized in that: Based on the numerical calculation results of multiple working conditions, the parametric evaluation of the destructive effectiveness of the blasting warhead on different targets is realized, and the simulation calculation process is visualized. Different damage levels are given to the damaged area of ​​the target. The calculation results are analyzed and an evaluation description of the destructive effect is given, and the conditions for complete destruction of the target are given.

Images