A method for constructing an ignition growth model of high-energy viscoelastic materials and related equipment

Abstract

The present application discloses a method for constructing an ignition growth model of a high-energy viscoelastic material and related equipment, the method comprising: constructing an initial model of the ignition growth of a high-energy viscoelastic material; obtaining a critical detonation probability interval of the high-energy viscoelastic material through a partition test; using a preset coupling algorithm to numerically simulate the minimum and maximum values ​​in the critical detonation probability interval, respectively, to generate simulation results; based on the simulation results, adjusting the model parameters of the initial model of the ignition growth of the high-energy viscoelastic material to obtain an ignition growth model of the high-energy viscoelastic material. The present application proposes a construction scheme for an ignition growth model of a high-energy viscoelastic material, obtains a variety of measured data based on a small partition test, and adjusts the model parameters in combination with numerical calculation and parameter fitting methods to improve the accuracy and applicability of the ignition growth model.

Description

Technical Field

[0001] The present application relates to the technical field of material testing, and in particular to a method for constructing an ignition growth model of a high-energy viscoelastic material and related equipment. Background Art

[0002] High-energy solid power devices are one of the main applications of high-energy viscoelastic materials. Under certain circumstances, such as operational errors and external shocks, the shock-to-detonation (XDT) phenomenon may occur, leading to serious consequences. As a mainstream method, the partition test is widely used in evaluating the impact sensitivity of solid viscoelastic materials. Its core is to explore its response ability to shock waves through small-charge tests, and to judge the impact effects it may suffer in the real power device environment. In the partition test, the dominant factor of shock initiation is the initial shock wave, and the shock wave reflected by the verification plate does not have a significant effect on shock initiation. Obviously, the shock initiation in the partition test is mainly determined by the initial shock wave.

[0003] The ignition growth model is the core issue of simulating the impact initiation process. It directly affects the accuracy and reliability of the simulation results and is crucial to simulating the impact initiation process. However, current research faces the challenge of the limited applicability of the ignition growth model. For many explosive materials, there is still a lack of a unified and accurate ignition growth model, especially for heterogeneous energetic materials such as solid viscoelastic materials, the research on the ignition growth model is still in an imperfect state. Summary of the invention

[0004] Based on this, it is necessary to address at least one of the problems mentioned above. The present invention provides a method for constructing an ignition growth model for high-energy viscoelastic materials and related equipment, and proposes a method for constructing an ignition growth model for high-energy viscoelastic materials. Based on a small partition test, a variety of measured data are obtained, and the model parameters are adjusted in combination with numerical calculation and parameter fitting methods to improve the accuracy and applicability of the ignition growth model.

[0005] In a first aspect, the present application provides a method for constructing an ignition growth model of a high-energy viscoelastic material, comprising the following steps:

[0006] Construct an initial model for the ignition growth of high-energy viscoelastic materials;

[0007] Obtain the critical detonation probability interval of high-energy viscoelastic materials through diaphragm tests;

[0008] Using a preset coupling algorithm to perform numerical simulation on the minimum value and the maximum value in the critical detonation probability interval respectively, and generate simulation results;

[0009] Based on the simulation results, the model parameters of the initial model of ignition growth of high-energy viscoelastic material are adjusted to obtain the ignition growth model of high-energy viscoelastic material.

[0010] Further, in some embodiments of the present application, the construction of the initial model of ignition growth of high-energy viscoelastic materials includes:

[0011] Calculating the number of hot spots and the ignition probability in the hot spot formation and ignition stages corresponding to the high-energy viscoelastic material to construct the first stage of the ignition growth initial model;

[0012] Calculating the number of hot spots and the reaction area in the hot spot growth stage to construct the second stage of the ignition growth initial model;

[0013] The energy data of the reaction completion stage is calculated to construct the third stage of the ignition growth initial model.

[0014] Further, in some embodiments of the present application, the step of obtaining the critical detonation probability interval of the high-energy viscoelastic material through a diaphragm test includes:

[0015] respectively setting corresponding partitions of different thicknesses for a plurality of high-energy viscoelastic materials of the same size and shape;

[0016] Conduct impact tests on each high-energy viscoelastic material using a standardized impact source to generate impact test results corresponding to partitions of different thicknesses;

[0017] Calculating the detonation probability at each partition thickness based on the impact test results;

[0018] Based on the impact test results and detonation probabilities corresponding to partitions of different thicknesses, a critical detonation probability interval of the high-energy viscoelastic material is determined.

[0019] Furthermore, in some embodiments of the present application, the step of obtaining the critical detonation probability interval of the high-energy viscoelastic material through the partition test further includes:

[0020] Constructing a partition test geometric model, wherein the partition test geometric model includes a main propellant, a partition material, a propellant column, a steel casing and a witness plate material;

[0021] Constructing a partition test material model, wherein the partition test material model includes a main charge material model, a partition material model, a steel casing material model and a witness plate material model;

[0022] The critical detonation probability interval of high-energy viscoelastic materials is simulated based on the partition test geometric model and the partition test material model.

[0023] Further, in some embodiments of the present application, the preset coupling algorithm is used to perform numerical simulation on the minimum value and the maximum value in the critical detonation probability interval respectively to generate simulation results, including:

[0024] Performing numerical simulation on a partition model corresponding to a minimum value in the critical detonation probability interval to generate a first simulation result in the simulation results, wherein the first simulation result is used to simulate an experimental result of the partition model corresponding to the minimum value;

[0025] A numerical simulation is performed on the partition model corresponding to the maximum value in the critical detonation probability interval to generate a second simulation result in the simulation results, wherein the second simulation result is used to simulate the experimental result of the partition model corresponding to the maximum value.

[0026] Furthermore, in some embodiments of the present application, the preset coupling algorithm is a smooth particle fluid dynamics-finite element coupling algorithm, and the simulation results include a numerical simulation pressure cloud map and pressure time history curves corresponding to multiple observation points.

[0027] Further, in some embodiments of the present application, after adjusting the model parameters of the initial model of ignition growth of the high-energy viscoelastic material based on the simulation results to obtain the ignition growth model of the high-energy viscoelastic material, the method further includes:

[0028] Calculate the corresponding theoretical analysis results through the pre-built shock wave parameter calculation model;

[0029] Obtaining test results corresponding to the diaphragm test;

[0030] Compare the experimental test results, the simulation results and the theoretical analysis results to obtain corresponding comparison results;

[0031] Verifying the ignition growth model of the high-energy viscoelastic material based on the comparison result to generate a verification result;

[0032] The ignition growth model of the high-energy viscoelastic material is adjusted secondary based on the verification results.

[0033] In a second aspect, the present application provides a device for constructing an ignition growth model of a high-energy viscoelastic material, comprising:

[0034] Building module for constructing the initial model of ignition growth of high-energy viscoelastic materials;

[0035] The test module is used to obtain the critical detonation probability interval of high-energy viscoelastic materials through partition tests;

[0036] A simulation module, used to perform numerical simulation on the minimum value and the maximum value in the critical detonation probability interval respectively by using a preset coupling algorithm to generate simulation results;

[0037] The adjustment module is used to adjust the model parameters of the initial model of the ignition growth of the high-energy viscoelastic material based on the simulation results to obtain the ignition growth model of the high-energy viscoelastic material.

[0038] In a third aspect, the present application also provides a computer device, comprising a memory, a processor, and a computer program stored in the memory and executable on the processor, wherein when the processor executes the computer program, the steps of the method for constructing an ignition growth model of a high-energy viscoelastic material as described in the first aspect are implemented.

[0039] In a fourth aspect, the present application further provides a computer-readable storage medium storing a computer program, which, when executed by a processor, implements the steps of the method for constructing an ignition growth model of a high-energy viscoelastic material as described in the first aspect.

[0040] As described above, the present application provides a method for constructing an ignition growth model of a high-energy viscoelastic material and related equipment, the method comprising: constructing an initial model of the ignition growth of a high-energy viscoelastic material; obtaining a critical detonation probability interval of the high-energy viscoelastic material through a partition test; using a preset coupling algorithm to numerically simulate the minimum and maximum values ​​in the critical detonation probability interval, respectively, to generate simulation results; based on the simulation results, adjusting the model parameters of the initial model of the ignition growth of the high-energy viscoelastic material to obtain an ignition growth model of the high-energy viscoelastic material. It can be seen that the present application proposes a construction scheme for an ignition growth model for high-energy viscoelastic materials, based on a small partition test, adding a variety of data collection methods to obtain a variety of measured data, combining numerical calculation and parameter fitting methods, adjusting the model parameters, improving the accuracy and applicability of the constructed ignition growth model, providing a physical model and numerical algorithm basis for the research on impact initiation of high-energy viscoelastic materials, and more economical, simple and fast than traditional methods. BRIEF DESCRIPTION OF THE DRAWINGS

[0041] In order to more clearly illustrate the embodiments of the present application or the technical solutions in the prior art, the drawings required for use in the embodiments or the description of the prior art will be briefly introduced below. Obviously, the drawings described below are only some embodiments of the present application. For ordinary technicians in this field, other drawings can be obtained based on these drawings without paying creative work.

[0042] Figure 1 It is a flow chart of a method for constructing an ignition growth model of a high-energy viscoelastic material provided in an embodiment of the present application;

[0043] Figure 2 It is a schematic diagram of the geometric model and mesh division of the small partition test provided in the embodiment of the present application;

[0044] Figure 3 is a schematic diagram of a numerical simulation process of a small partition test provided in an embodiment of the present application;

[0045] Figure 4 The model parameters of the 25.72 mm thick organic glass partition provided in the embodiment of this application are x =5.5 and the model is calculated to t= 30μs pressure cloud map;

[0046] Figure 5 The model parameters of the 25.72 mm thick organic glass partition provided in the embodiment of this application are x =5.5 when the pressure-time curve of the observation point AN;

[0047] Figure 6 The model parameters of the 25.72 mm thick organic glass partition provided in the embodiment of this application are x =5.6 and the model is calculated to t= 30μs pressure cloud map;

[0048] Figure 7 The model parameters of the 25.72 mm thick organic glass partition provided in the embodiment of this application are x =5.6 when the pressure-time curve of observation point AN;

[0049] Figure 8 The model parameters of the 28.4 mm thick organic glass partition provided in the embodiment of this application are x =5.5 and the model is calculated to t = 31μs pressure cloud map;

[0050] Fig. 9 The model parameters of the 28.4 mm thick organic glass partition provided in the embodiment of this application are x =5.5 when the pressure-time curve of the observation point AN;

[0051] Fig.10 The model parameters of the 28.4 mm thick organic glass partition provided in the embodiment of this application are x =5.4 and the model is calculated to t= 31μs pressure cloud map;

[0052] Fig.11 The model parameters of the 28.4 mm thick organic glass partition provided in the embodiment of this application are x =5.4 when the pressure-time curve of observation point AN;

[0053] Fig.12 It is a structural schematic diagram of a shock wave parameter calculation model provided in an embodiment of the present application;

[0054] Fig.13 is a comparison diagram of detonation pressures provided in the embodiments of the present application;

[0055] Fig.14 It is a structural schematic diagram of a device for constructing an ignition growth model of a high-energy viscoelastic material provided in an embodiment of the present application;

[0056] Fig.15 It is a schematic diagram of the structure of a computer device provided in an embodiment of the present application. DETAILED DESCRIPTION

[0057] Exemplary embodiments will be described in detail herein, examples of which are shown in the accompanying drawings. When the following description refers to the drawings, the same numbers in different drawings represent the same or similar elements unless otherwise indicated. The embodiments described in the following exemplary embodiments do not represent all embodiments consistent with the present application. Instead, they are merely examples of systems and methods consistent with the examples detailed in the attached claims, or with some aspects of the present application.

[0058] It should be noted that, in this article, the terms "include", "comprises" or any other variants thereof are intended to cover descriptions such as non-exclusive inclusion, so that a process, method, article or device including a series of elements includes not only those elements, but also other elements not explicitly listed, or also includes elements inherent to such process, method, article or device. In the absence of further restrictions, an element defined by the sentence "includes a ..." does not exclude the existence of other identical elements in the process, method, article or device including the element. In addition, components, features, and elements with the same name in different embodiments of the present application may have the same meaning or different meanings, and their specific meanings need to be determined by their explanation in the specific embodiment or further combined with the context of the specific embodiment.

[0059] It should be understood that the specific embodiments described herein are only used to explain the present application and are not used to limit the present application.

[0060] In the subsequent description, the suffixes such as "module", "component" or "unit" used to represent elements are only used to facilitate the description of the present application, and have no specific meanings. Therefore, "module", "component" or "unit" can be used in a mixed manner.

[0061] High-energy solid power devices are one of the main applications of high-energy viscoelastic materials. Under certain circumstances, such as operational errors and external impacts, shock-to-detonation phenomena may occur, leading to serious consequences. As a mainstream method, the partition test is widely used in evaluating the impact sensitivity of solid viscoelastic materials. Its core is to explore its response ability to shock waves through small-dose tests, and to judge the impact impact it may suffer in the real power device environment. At present, the classification of viscoelastic material hazard levels in China generally follows the "U.S. Naval Ordnance Laboratory (NOL) Partition Test Standard", which plays an important role in promoting the development of industry norms. It is found that in the partition test, the dominant factor of shock initiation is the initial shock wave, and the shock wave reflected by the verification plate has no significant effect on shock initiation. Based on the above research, it can be concluded that the shock initiation in the partition test is mainly determined by the initial shock wave.

[0062] The ignition growth model is the core issue in simulating the impact detonation process, which directly affects the accuracy and reliability of the simulation results. However, current research faces the challenge of the limited applicability of the ignition growth model. For many explosive materials, there is still a lack of a unified and accurate ignition growth model, especially for heterogeneous energetic materials such as solid viscoelastic materials, the research on the ignition growth model is still in an imperfect state.

[0063] In order to solve the above technical problems, the present application provides a method for constructing an ignition growth model of a high-energy viscoelastic material and related equipment, which can obtain a variety of measured data based on a small partition test, and adjust the model parameters in combination with numerical calculation and parameter fitting methods to improve the accuracy and applicability of the ignition growth model.

[0064] See also Figure 1 , Figure 1 1 is a flow chart of a method for constructing an ignition growth model of a high-energy viscoelastic material provided in an embodiment of the present application. The method for constructing an ignition growth model of a high-energy viscoelastic material may specifically include the following steps:

[0065] S1. Construct the initial model of ignition growth of high energy viscoelastic materials;

[0066] Specifically, for step S1, the initial ignition growth model of the high-energy viscoelastic material is first constructed; in order to accurately describe the physical process of the high-energy viscoelastic material during impact, the Lee-Tarver (trinomial) ignition and growth model is used:

[0067] In the formula, For time, For pressure, , is the density. is a constant, and F represents the reaction progress variable.

[0068] Furthermore, in some embodiments, step S1 of "constructing an initial model of ignition growth of high-energy viscoelastic material" may specifically include:

[0069] S11. Calculate the number of hot spots and the ignition probability during the hot spot formation and ignition stages corresponding to the high-energy viscoelastic material to construct the first stage of the initial model of ignition growth;

[0070] S12. Calculate the number of hot spots and the reaction area in the hot spot growth stage to construct the second stage of the ignition growth initial model;

[0071] S13. Calculate the energy data of the reaction completion stage to construct the third stage of the ignition growth initial model.

[0072] Specifically, for step S1, the three-term Lee-Tarver ignition growth model is a mathematical model used to describe the process of solid viscoelastic materials changing from a static state to a detonation state when subjected to external impact (such as collision or explosion). This model consists of three main parts, corresponding to different stages of the ignition growth process, each of which has its own unique characteristics and reaction laws:

[0073] (1) Hotspot formation and ignition. In this stage, the viscoelastic material is affected by the shock wave, and a high-temperature and high-pressure area is generated inside, which is the so-called "hotspot". The formation of these hotspots is closely related to the pressure of the shock wave. The greater the pressure, the more hotspots are formed and the higher the probability of ignition.

[0074] (2) Hotspots grow slowly. In this stage, the hotspots that have been formed begin to grow slowly, but the growth rate is relatively slow. This is because in this stage, the energy of the shock wave is mainly used to maintain the existence and growth of the hotspots, but fails to trigger large-scale chemical reactions. However, as the hotspots continue to increase and expand, they begin to gradually connect together to form a larger reaction area.

[0075] (3) Rapid reaction completion stage. When the hot spot grows to a certain extent, the chemical reaction of the viscoelastic material begins to proceed rapidly, and the unreacted viscoelastic material quickly decomposes, releasing a large amount of energy. This process is extremely rapid and is completed almost instantly. This stage marks the transition of the viscoelastic material from the detonation state to the detonation state, generating huge explosive force and destructive power.

[0076] Among them, the original values ​​of various parameters of high-energy viscoelastic materials in the ignition growth model are shown in Table 1.

[0077] Table 1 Parameters of ignition growth model for high energy viscoelastic materials

[0078]

[0079] S2. Obtain the critical detonation probability interval of high energy viscoelastic materials through partition test;

[0080] Specifically, for step S2, after using the original values ​​of the parameters in Table 1 to perform numerical simulation on the small partition test of high-energy viscoelastic material, the simulation results obtained are obviously inconsistent with the test results. It can be judged that the original values ​​of the parameters are not completely applicable to the high-energy viscoelastic material used, so it is necessary to re-fit these parameters to determine their accurate values, which is the main purpose of parameter fitting. The specific test results are shown in Table 2.

[0081] Table 2 Small partition test results

[0082]

[0083] At present, there are two methods to evaluate the impact sensitivity of solid viscoelastic materials through small partition tests. One method is to obtain the 50% detonation partition thickness of the solid viscoelastic material by the lifting method. , which means that the probability of high-energy viscoelastic material detonating is 50% when the partition thickness is 27mm. As the partition thickness becomes smaller than 27mm, the probability of viscoelastic material detonating becomes lower and lower until it drops to 0%; as the partition thickness becomes larger than 27mm, the probability of viscoelastic material detonating becomes higher and higher until it increases to 100%. This method is generally used to compare the sensitivity of two different viscoelastic materials.

[0084] Another method is to obtain the critical detonation probability interval of viscoelastic materials through the Langley method. The critical detonation probability interval of high-energy viscoelastic materials obtained in the experiment is 25.72mm-28.40mm, which means that when the thickness of the partition is less than or equal to 25.72, the probability of viscoelastic material detonation is 100%, when the thickness of the partition is greater than or equal to 28.4mm, the probability of viscoelastic material detonation is 0%, and when the thickness of the partition is between 25.72mm-28.40mm, the viscoelastic material may or may not detonate, and the detonation probability is determined by the thickness of the partition. This method is generally used to determine the impact initiation threshold of solid viscoelastic materials.

[0085] Since the ultimate goal of this application is to evaluate the impact initiation characteristics of solid viscoelastic materials and the impact safety of SRM fragments, the embodiments of this application refer to the critical detonation probability interval obtained by the experiment, use the SPH-FEM coupling algorithm, and perform numerical simulations of partition thicknesses of 25.72mm and 28.4mm to determine the parameters of the ignition growth model.

[0086] For the SPH-FEM coupling algorithm, the SPH-FEM coupling algorithm is a numerical calculation method that combines two different numerical discretization techniques: smoothed particle hydrodynamics (SPH) and finite element method (FEM). This coupling algorithm is designed to simulate complex physical phenomena, especially when it involves material failure, large deformation, and multi-physics field interactions. In the study of impact initiation of high-energy viscoelastic materials, the SPH-FEM coupling algorithm can help researchers more accurately simulate and analyze the behavior of viscoelastic materials when subjected to impact loads. The key points of the coupling algorithm include:

[0087] Data Exchange: Exchange boundary conditions and force information between SPH and FEM domains to ensure that the physics between the two methods is continuous and consistent.

[0088] Load transfer: Ensure that the loads transferred from the FEM region to the SPH region and the loads transferred from the SPH region to the FEM region match.

[0089] Stability and accuracy: The coupled algorithm needs to ensure the stability and accuracy of the numerical calculation, which may involve choosing an appropriate time step, particle spacing, and element size.

[0090] Computational efficiency: The coupling algorithm needs to balance computational efficiency and computational accuracy because SPH and FEM may differ in computational cost.

[0091] Simulation of physical phenomena: The coupled algorithm can simulate a variety of physical phenomena including shock wave propagation, material destruction, thermodynamic effects and chemical reactions.

[0092] In the study of impact initiation of high-energy viscoelastic materials, the SPH-FEM coupling algorithm can provide a detailed analysis of viscoelastic materials when subjected to impact loads, including ignition, propagation of detonation waves, dynamic response of materials, and destruction of structures. Through the coupling algorithm, researchers can better understand and predict the behavior of viscoelastic materials, thereby providing support for the design and safety assessment of viscoelastic materials.

[0093] Furthermore, in some embodiments, step S2 of “obtaining a critical detonation probability interval of a high-energy viscoelastic material through a partition test” may specifically include:

[0094] S21. respectively setting corresponding partitions of different thicknesses for a plurality of high-energy viscoelastic materials of the same size and shape;

[0095] S22. Performing an impact test on each high-energy viscoelastic material using a standardized impact source to generate impact test results corresponding to partitions of different thicknesses;

[0096] S23. Calculate the detonation probability for each partition thickness based on the impact test results;

[0097] S24. Based on the impact test results and detonation probabilities corresponding to partitions of different thicknesses, determine the critical detonation probability range of the high-energy viscoelastic material.

[0098] Specifically, for step S2, the critical detonation probability interval of the high-energy viscoelastic material is obtained through the partition test. The small partition test is an experimental method for evaluating the impact sensitivity of energetic materials (such as high-energy viscoelastic materials). The test adjusts the intensity of the shock wave by changing the thickness of the partition, so as to determine the critical condition for the high-energy viscoelastic material to detonate.

[0099] In this embodiment, the Langley method is mainly used to obtain the critical detonation probability interval of high-energy viscoelastic materials. The Langley method (Von Neumann's method), proposed by John von Neumann, is an experimental method for determining the impact initiation characteristics of energetic materials (such as high-energy viscoelastic materials). This method obtains the critical detonation probability interval of viscoelastic materials by changing the impact intensity in the partition test. The following are the steps to obtain the critical detonation probability interval using the Langley method:

[0100] Experimental design: Prepare multiple viscoelastic material samples of the same size and shape. Prepare a separator for each sample, and the thickness of the separator can be adjusted as needed.

[0101] Experimental setup: Place each viscoelastic material sample behind a baffle that attenuates the shock wave generated by the explosive source. Set up a standardized shock source such as an explosive charge or other device that can produce a consistent shock wave.

[0102] Experimental process: Each sample was subjected to impact test, and the thickness of the partition was changed to change the intensity of the shock wave.

[0103] Data Logging: For each different diaphragm thickness, record whether the viscoelastic material detonates. Typically, multiple tests are performed using a range of diaphragm thicknesses to collect sufficient statistical data.

[0104] Probability analysis: According to the ratio of the number of detonations of viscoelastic materials to the total number of tests, the probability of detonation at each partition thickness is calculated.

[0105] Determine the critical interval: By analyzing the relationship between the detonation probability and the partition thickness, the critical detonation probability interval is determined. The lower limit of this interval is the minimum partition thickness at which the viscoelastic material detonates 100%, and the upper limit is the maximum partition thickness at which no detonation occurs.

[0106] Statistical analysis: Use statistical methods (such as cumulative frequency distribution) to determine the confidence interval of the detonation probability, thereby obtaining a more accurate critical detonation probability interval.

[0107] Model fitting: Experimental data is used to fit the parameters of the ignition growth model, which determine the impact sensitivity and detonation characteristics of viscoelastic materials.

[0108] Application of results: The obtained critical detonation probability interval is very important for designing safer viscoelastic material formulations, improving protective measures and understanding the failure mechanism of materials.

[0109] The Langley method is an empirical method that relies on a large amount of experimental data and statistical analysis to determine the impact initiation characteristics of high-energy viscoelastic materials. This method is very useful in the research and application of high-energy viscoelastic materials, especially in evaluating the safety and sensitivity of viscoelastic materials. With this method, researchers can better predict and control the safety risks of high-energy viscoelastic materials in practical applications.

[0110] Furthermore, in some embodiments, step S2 of “obtaining a critical detonation probability interval of a high-energy viscoelastic material through a partition test” may specifically include:

[0111] Construct a geometric model of the partition test, which includes the main propellant, the partition material, the propellant column, the steel casing and the witness plate material;

[0112] Construct a partition test material model, which includes a main charge material model, a partition material model, a steel casing material model and a witness plate material model;

[0113] The critical detonation probability interval of high-energy viscoelastic materials is simulated based on the partition test geometry model and the partition test material model.

[0114] Specifically, the geometric model of the partition test includes: the main explosive is Pentolite explosive, with a size of φ50mm×50mm, and is divided by SPH particles. The average particle spacing is 5mm, and there are 800 SPH particles in total; the partition material is organic glass (PMMA), with a diameter of 5mm, and is divided by hexahedral grid; the fired charge is a high-energy viscoelastic material, with a size of φ36mm×140mm, and is divided by SPH particles. The particle spacing is 5mm, and there are 1836 particles in total; the steel sleeve is a No. 45 precision drawn tube, with a size of φ(48mm-36mm)×140mm, and is divided by hexahedral grid, with a grid size of 5mm, and there are 1216 grids in total; the witness plate material is No. 45 steel, with a size of φ150mm×10mm, and the grid is divided into two parts. The middle φ25mm area is divided by SPH particles, with a particle spacing of 50mm, and there are 160 particles in total. The outer area uses a hexahedral grid with a grid size of 50mm and 320 units.

[0115] like Figure 2 As shown, Figure 2The geometric model and mesh division of the small partition test are provided, and the initial position of the reference particles set by AN is convenient for data analysis.

[0116] Specifically, the material model for the partition test includes:

[0117] (1) Main grain material model

[0118] Pentolite explosives are described by the JWL (Jones-Wilkins-Lee) equation of state:

[0119] in It's pressure. is the relative volume, is the volume, is the main grain material density, is the specific internal energy, is a constant, e is the specific internal energy. The specific parameters are shown in Table 3.

[0120] Table 3 Explosive material parameters

[0121]

[0122] (2) Material models of diaphragms, steel casings and witness plates

[0123] Using the elastoplastic fluid dynamics material model:

[0124] In the formula, is the volumetric strain rate sum of the strain rates in all directions, is the stress deviator, is the deviatoric stress, p It's pressure. σ ij is the yield strength of the material. is the Kronecker delta function, which is 1 when i=j and 0 otherwise. It is used to represent the sum of diagonal elements in the equation, i.e., the calculation of volumetric strain rate; represents the strain rate tensor, which represents the strain rate components of the material in the i and j directions.

[0125] The von Mises yield criterion can be expressed as follows:

[0126] Using the Gruneisen equation of state (thermodynamic equation of the state of a material under high pressure):

[0127] in is the density, is a constant, E It is internal energy. p It's pressure. μ is the compression ratio. The Gruneisen parameters of PMMA and 45# steel are shown in Table 4.

[0128] Table 4 Material parameters of Gruneisen equation

[0129]

[0130] like Figure 3 As shown, Figure 3 Provides the numerical simulation process of small partition test. Figure 3 The test procedure for a PMMA separator with a thickness of 25.6 mm is described. t= 0μs, start to detonate; t= 5μs, PMMA produces projected shock waves; t= 10μs, the delivered drug receives the shock wave; t= 15μs, the viscoelastic material is detonated; t= 20μs, the detonation wave squeezes the steel casing outwards; t= 30μs, the detonation wave reaches the witness plate; t= 35μs, witnessing deformation and damage of the plate; t= At 55μs, the steel casing was completely ruptured and the steel witness plate was penetrated.

[0131] S3. using a preset coupling algorithm to perform numerical simulations on the minimum and maximum values ​​in the critical detonation probability interval, respectively, to generate simulation results;

[0132] Further, in some embodiments, step S3 of "using a preset coupling algorithm to perform numerical simulation on the minimum value and the maximum value in the critical detonation probability interval respectively to generate simulation results" may specifically include:

[0133] S31. Performing numerical simulation on the partition model corresponding to the minimum value in the critical detonation probability interval to generate a first simulation result in the simulation results, the first simulation result being used to simulate the experimental result of the partition model corresponding to the minimum value;

[0134] S32. Perform numerical simulation on the partition model corresponding to the maximum value in the critical detonation probability interval to generate a second simulation result in the simulation results, and the second simulation result is used to simulate the experimental result of the partition model corresponding to the maximum value.

[0135] Furthermore, in some embodiments, the preset coupling algorithm in step S3 is a smooth particle fluid dynamics-finite element coupling algorithm, and the simulation results include a numerical simulation pressure cloud map and pressure time history curves corresponding to multiple observation points.

[0136] Specifically, for step S3, the physical meaning of the three terms of the ignition growth model shows that the ignition term is the key term for controlling the detonation reaction process of high-energy viscoelastic materials under impact. I, b, a, x The specific method is to use the lifting method to fit these four parameters. I, b, a, x The four parameters are fitted one by one to ensure that the calculation process and results of the numerical simulation are consistent with the experiment, so as to identify the main parameters affecting the impact sensitivity of high-energy viscoelastic materials and determine the accurate values ​​of the main parameters.

[0137] After preliminary fitting, it is found that the parameters I, b, a The impact on the impact sensitivity of high-energy viscoelastic materials is small, so these three parameters can retain their original values ​​without modification, namely: I=45000, b= 0.667 ，a= 0.02. But the parameter x The impact sensitivity of high-energy viscoelastic materials is greatly affected, and the parameters x The value is negatively correlated with the impact sensitivity of high-energy viscoelastic materials, that is, the parameter x The smaller the value, the higher the impact sensitivity of the high-energy viscoelastic material, and thus the parameter x It is the main parameter that affects the impact sensitivity of high-energy viscoelastic materials. Therefore, the fitting parameters x The accurate value of is the focus of this parameter fitting work, and is also the basic work for the subsequent construction of the high-energy viscoelastic material injection detonation model. x The specific fitting work is as follows:

[0138] 1. Perform numerical simulation on the partition model corresponding to the minimum value in the critical detonation probability interval, that is, the numerical simulation of 25.72 mm thick organic glass partition:

[0139] (1) Parameters x Value: 5.5

[0140] For the 25.72 mm thick plexiglass partition model, when the parameter x= 5.5, model operation to t= 30μs, through Figure 4 It can be seen from the pressure cloud diagram that the shock wave caused by the detonation of the main propellant has caused the detonation of the secondary propellant, which is consistent with the test results.

[0141] like Figure 5As shown in the figure, the pressure time history curve AN of the 14 observation points is analyzed again. It can be seen from the pressure time history curve that t= At 12 μs, the pressure value at observation point A suddenly increased to about 6.02 GPa, which indicates that an explosion occurred in the area near the observation point under the impact of the shock wave; t= 12μs-17μs, the shock wave passes through points B, C, D, and E in turn, and the pressure increases from 6.02GPA to 20.6GPA, which explains the growth process of the detonation; t= Starting from 17μs, the detonation pressure was maintained at about 20.6GPA, and the viscoelastic material underwent a stable detonation process; t= At 26μs, the observation point N closest to the witness plate has a steep pressure jump, reaching 27GPa, which is the result of the superposition of the shock wave reflected by the witness plate and the detonation wave. From the pressure time history curve of the observation point, the SPH-FEM coupling algorithm can better reflect the ignition growth process of solid viscoelastic materials, and the simulation results are consistent with the experimental results.

[0142] In summary, the information displayed by the pressure time history curve is basically consistent with that of the pressure cloud map, and both can accurately reflect the entire process of detonation and continuous explosion transmission of high-energy viscoelastic materials under impact.

[0143] (2) Parameters x Value: 5.6

[0144] For the 25.72 mm thick plexiglass partition model, when the parameter x= 5.6, model operation to t= 30μs, through Figure 6 It can be seen from the pressure cloud diagram that the shock wave caused by the detonation of the main propellant has reached the end of the propellant column, but has not yet caused the propellant to detonate. There is no deformation or damage trend in the center of the witness plate, which is obviously inconsistent with the test results.

[0145] like Figure 7 As shown in Figure 2, the pressure-time curve AN of the 14 observation points is analyzed again. It can be seen from the pressure-time curve that t= At 12μs, the pressure value at observation point A increased to about 6.09GPa, and then dropped rapidly, which indicates that the explosive was not detonated by the incident shock wave. t= 14μs to t= 50μs, the shock wave passes through the BE point in turn, and the pressure peak gradually decreases, which shows that the explosive has not produced a detonation growth process. From the pressure-time curve of the observation point, the high-energy viscoelastic material was not detonated when the partition thickness was less than 25.72mm. The simulation results are completely inconsistent with the test results, which proves that the parameters x The value of cannot be greater than or equal to 5.6.

[0146] 2. Perform numerical simulation on the partition model corresponding to the maximum value in the critical detonation probability interval, i.e., numerical simulation on a 28.4 mm thick plexiglass partition:

[0147] (1) Parameters x Value: 5.5

[0148] For the 28.4 mm thick plexiglass partition model, when the parameter x= 5.5, model operation to t= At 31μs, through Figure 8 It can be seen from the pressure cloud diagram that the shock wave caused by the detonation of the main propellant has reached the end of the propellant column, but has not yet caused the propellant to detonate. There is no deformation or damage trend in the center of the witness plate, which is consistent with the test results.

[0149] like Fig. 9 As shown in Figure 2, the pressure time history curve AN of the 14 observation points is analyzed again. It can be seen from the pressure time history curve that t= At 12μs, the pressure value at observation point A increased to about 4.17GPa, and then dropped rapidly, which indicated that the explosive was not detonated by the incident shock wave. t= 14μs to t=50 μs, the shock wave passes through the BE point in turn, and the pressure peak gradually decreases, which indicates that the explosive has not produced a detonation growth process. From the pressure time history curve of the observation point, the SPH-FEM coupling algorithm can better reflect the ignition growth process of solid viscoelastic materials, and the simulation results are consistent with the experimental results.

[0150] (2) Parameters x Value: 5.4

[0151] For the 28.4 mm thick plexiglass partition model, when the parameter x= 5.4, ​​model operation to t= At 31μs, through Fig.10 From the pressure cloud diagram shown, it can be seen that the shock wave caused by the detonation of the main propellant has caused the detonation of the secondary propellant, which is obviously inconsistent with the test results.

[0152] like Fig.11 As shown in Figure 1, the pressure time history curve AN of the 14 observation points (A is the first observation point, N is the last observation point) is analyzed again. It can be seen from the pressure time history curve that t= At 13μs, the pressure value at observation point A suddenly increased to about 4.17GPa; t= 13μs to t= 20μs, the shock wave passes through points BH in turn (B is the second observation point, H is the eighth observation point), and the pressure increases from 4.17GPa to 20.5GPa; t=Starting from 20μs, the detonation wave pressure is maintained at about 20.6GPa for a stable detonation; t = At 27μs, the observation point N closest to the witness plate has a steep pressure jump, reaching 27.2GPa. From the pressure-time curve of the observation point, the high-energy viscoelastic material completely detonates when the partition thickness is less than 28.4mm. The simulation results are completely inconsistent with the experimental results, which proves that the parameters x The value of cannot be less than or equal to 5.4.

[0153] S4. adjusting the model parameters of the initial model of the ignition growth of the high-energy viscoelastic material based on the simulation results to obtain the ignition growth model of the high-energy viscoelastic material;

[0154] Specifically, for step S4, after multiple rounds of numerical simulation, it is found that the parameter x The impact sensitivity of high-energy viscoelastic materials is greatly affected, and the parameters x The value of is negatively correlated with the impact sensitivity of the high-energy viscoelastic material, that is, the smaller the value of parameter x, the higher the impact sensitivity of the high-energy viscoelastic material. x It is the main parameter that affects the impact sensitivity of high-energy viscoelastic materials. The comparison between numerical simulation and experiment is shown in Table 5. L =25.72mm is the thickness of the complete detonation partition of the high-energy viscoelastic material. Under this partition thickness, when the parameter x= 5.5, the numerical simulation results show complete detonation. x= 5.6, the numerical simulation results show that no detonation occurs, which indicates that the parameters x Should be less than 5.6. L =28.4mm is the thickness of the high-energy viscoelastic material that does not cause detonation. Under this thickness, when the parameter x= 5.5, the numerical simulation results show that no detonation occurs. x= 5.4, ​​the numerical simulation results show complete detonation, which indicates that the parameters x Should be greater than 5.4. From this, the parameters can be determined x= 5.5.

[0155] Table 5 Ignition growth model parameter fitting key operating conditions and results

[0156]

[0157] After the above parameter fitting work, the precise values ​​of various parameters of high-energy viscoelastic materials in the ignition growth model can be obtained, as shown in Table 6.

[0158] Table 6 Accurate values ​​of parameters of the ignition growth model for high-energy viscoelastic materials

[0159]

[0160] Furthermore, in some embodiments, after step S4 "adjusting the model parameters of the initial model of the ignition growth of the high-energy viscoelastic material based on the simulation results to obtain the ignition growth model of the high-energy viscoelastic material", the method may specifically include: S5. Verifying and adjusting the ignition growth model of the high-energy viscoelastic material.

[0161] Wherein, step S5 may specifically include:

[0162] S51. Calculate the corresponding theoretical analysis results through the pre-built shock wave parameter calculation model;

[0163] S52. Obtain the test results corresponding to the partition test;

[0164] S53. Compare the experimental test results, simulation results and theoretical analysis results to obtain corresponding comparison results;

[0165] S54. Verify the ignition growth model of high energy viscoelastic materials based on the comparison results and generate verification results;

[0166] S55. Make secondary adjustments to the ignition growth model of high-energy viscoelastic materials based on the verification results.

[0167] In this embodiment, a verification method for the ignition growth model of high-energy viscoelastic materials is also provided. After the main propellant column is detonated, the generated shock wave undergoes a series of complex physical processes and finally achieves the impact detonation of the propellant. Specifically, this process can be divided into the following four steps:

[0168] (1) By detonating the main propellant column, huge energy is released, thus generating a detonation wave.

[0169] (2) The detonation wave propagates forward and passes through the contact interface between the main charge and the partition. During this process, part of the energy is reflected by the partition, and the remaining part forms the partition-incident shock wave.

[0170] (3) The incident shock wave propagates inside the partition. Due to the damping and attenuation of the partition material, the energy of the shock wave gradually weakens during the propagation process. When the shock wave reaches the bottom of the partition, its intensity has been significantly reduced, forming a partition outgoing shock wave.

[0171] (4) The outgoing shock wave of the partition passes through the interface between the partition and the explosive, and energy transfer and conversion occur again. In this process, the shock wave forms the incident shock wave of the explosive, which in turn triggers the shock detonation of the explosive.

[0172] After four steps, the numerical solution of the incident shock wave of the explosive can be obtained. The calculation process is as follows: Fig.12 As shown in the figure P is the shock wave pressure, D is the shock wave velocity.

[0173] like Fig.12 As shown, Fig.12 A shock wave parameter calculation model is provided. When the main charge type and the parameters of the partition material are known, the above four physical processes can be solved by theoretical calculation according to the theory of explosive detonation and shock wave propagation. Therefore, this section calculates the theoretical value of the incident shock wave pressure of the charged charge, and compares the experimental data, numerical simulation results and analytical solution of the incident shock wave pressure of the charged charge in the small partition test to verify the accuracy of the shock initiation criterion and ignition growth model of solid viscoelastic materials constructed in this chapter, and verify the applicability of the SPH-FEM coupling algorithm in the shock initiation problem of solid viscoelastic materials.

[0174] 1. Calculation of detonation pressure and detonation velocity of Pentiolite explosive

[0175] Pentiolite explosives are made from a mixture of two components, PETN and TNT. The method can calculate the detonation pressure and velocity of Pentiolite explosives.

[0176] (1) Pentiolite CJ Explosion Pressure

[0177] Among them, the initial density of explosives , It's pressure. is the adiabatic index, It's explosive speed, , is the specific heat at constant pressure, is the specific heat at constant volume, is the overall adiabatic index, is the amount of substance, is the mass ratio, is the adiabatic index, is the mole fraction.

[0178] (2) Detonation velocity of pentiolite explosive

[0179] in, Pentiolite explosives Characteristic calorific value, potential energy factor and mass fraction of components.

[0180] The material parameters of each component of Pentiolite explosive are shown in Table 7.

[0181] Table 7 Material parameters of TNT and PETN

[0182]

[0183] According to the formula of detonation velocity and detonation pressure, combined with the material parameters in Table 3.12, the CJ state of Pentiolite explosive can be calculated. , .

[0184] 2. Parameters of the incident shock wave on the partition

[0185] The propagation law of the explosion shock wave in the partition is affected by the relationship between the impact impedance of the explosive and the partition. When the impact impedance of the explosive is different from that of the partition, the shock wave will be attenuated; when the impact impedance of the explosive is consistent with that of the partition, the shock wave passing through the interface will not be attenuated.

[0186] The impact impedance of Pentolite explosive is greater than the shock wave impedance of PMMA (Polymethyl Methacrylate), that is, At this time, the detonation product is expressed as:

[0187] In the formula, It's explosive pressure. It's explosive speed, is the particle velocity, is the density, It is the speed of sound.

[0188] Interface movement speed It can be expressed as:

[0189] From the isentropic equation of the product:

[0190] Combining the above formulas, we can get:

[0191] Combining the formulas, we get:

[0192] Substituting into the above formula again, we can get:

[0193] Then there is

[0194] The pressure equation of the medium behind the wave is:

[0195]

[0196] The velocity of the shock wave is related to the velocity of the particle as follows:

[0197]

[0198] in, is the speed of sound at zero pressure, is a constant. The material parameters on both sides of the interface are shown in Table 8.

[0199] Table 8 Hugoniot parameters of organic glass and high-energy viscoelastic materials

[0200]

[0201] Combining the above formulas, the calculation formula for the incident pressure of the partition shock wave can be obtained:

[0202] From the above calculation, we can know that the detonation velocity of the main propellant is , explosion pressure , the multi-party index of mixed explosives is , we can get:

[0203]

[0204] From the above two formulas, the incident shock wave parameters of the PMMA partition can be obtained: , , as shown in Table 9.

[0205] Table 9 Parameters of incident shock wave on organic glass partition

[0206]

[0207] 3. Parameters of the shock wave emitted by the partition

[0208] The shock wave attenuation law in the partition is exponential, which can be known from the literature and experimental laws:

[0209] in The thickness is The explosion pressure, explosion velocity, particle velocity and local sound velocity of the shock wave emitted by the plexiglass partition.

[0210] 4. The incident pressure of the fired grain

[0211] The exit pressure of the organic glass partition is not the incident pressure of the discharged medicine. Due to the shock wave impedance, the shock wave still needs to be attenuated at the interface between the partition and the discharged medicine column before reaching the discharged medicine column. The shock impedance of PMMA is smaller than the shock impedance of high-energy viscoelastic materials, that is, , at this time the medicine is given to gain extra speed :

[0212]

[0213] According to the polynomial equation of gas, we can get:

[0214] Then we have:

[0215] Combining the above three formulas, we can get the expression of the detonation pressure of the explosive:

[0216]

[0217] 5. Comparative analysis of incident pressure

[0218] The measured results, theoretical analysis results and numerical simulation results of the incident pressure test of the high-energy viscoelastic material small diaphragm test were compared. Fig.13 shown.

[0219] from Fig.13 It can be seen that the results of the test, theoretical analysis and numerical simulation are more consistent when the thickness of the partition is greater than 27mm. This shows that the SPH-FEM method can accurately describe the physical process of shock wave propagation across the medium in solving the shock initiation problem, and verifies the advancement of the algorithm. When the thickness is less than 27mm, the test results and numerical simulation results are significantly improved. The reason for this is that during the small partition test, when the thickness of the partition is less than 27mm, the explosive is in a completely detonated state, that is, Fig.14 In the four working conditions marked by dotted lines, the explosive has been detonated by shock. The pressure values ​​of the test and numerical simulation are the result of the superposition of the incident shock wave and the detonation wave of the high-energy viscoelastic material. Therefore, there will be an obvious rise in the position where the thickness of the partition is less than 27 mm in the figure. The theoretical analysis only considers the physical transmission process of the shock wave, but does not consider the shock detonation effect of the explosive, so the curve is relatively smooth. The rise in the test results and numerical simulation results here illustrates the accuracy of the solid viscoelastic material shock detonation criterion and the fitted ignition growth model constructed in the embodiment of the present application.

[0220] In summary, the method for constructing the ignition growth model of high-energy viscoelastic materials provided in the embodiment of the present application is as follows: first, construct the initial model of ignition growth of high-energy viscoelastic materials; then, obtain the critical detonation probability interval of high-energy viscoelastic materials through the partition test; then, use the preset coupling algorithm to numerically simulate the minimum and maximum values ​​in the critical detonation probability interval respectively, and generate simulation results; finally, adjust the model parameters of the initial model of ignition growth of high-energy viscoelastic materials based on the simulation results to obtain the ignition growth model of high-energy viscoelastic materials. It can be seen that this embodiment fits the parameters of the ignition growth model of solid viscoelastic materials through the information such as the shock wave incident pressure and the detonation velocity of solid viscoelastic materials obtained during the test, that is, establishes a numerical model of the detonation growth stage of the solid power device. The model is added to the SPH-FEM coupling algorithm, and the accuracy of the fitting results of the ignition growth model and the advancement of the SPH-FEM coupling algorithm for the impact initiation problem are verified through the numerical calculation of the test conditions and the comparison of the critical partition thickness and the numerical calculation results obtained in the test, so as to provide a physical model and numerical algorithm basis for the research on the impact initiation of solid viscoelastic materials.

[0221] In order to facilitate the implementation of the high-energy viscoelastic material ignition growth model construction method of the embodiment of the present application, the embodiment of the present application also provides a construction device based on the above-mentioned high-energy viscoelastic material ignition growth model construction method. The meanings of the terms are the same as those in the above-mentioned high-energy viscoelastic material ignition growth model construction method, and the specific implementation details can refer to the description in the method embodiment.

[0222] See also Fig.14 , Fig.14 The schematic diagram of the structure of the device for constructing the ignition growth model of a high-energy viscoelastic material provided in the embodiment of the present application includes a construction module 201, an experimental module 202, a simulation module 203 and an adjustment module 204, which can be specifically as follows:

[0223] Construction module 201, used to construct an initial model of ignition growth of high energy viscoelastic materials;

[0224] The test module 202 is used to obtain the critical detonation probability interval of the high-energy viscoelastic material through a partition test;

[0225] A simulation module 203 is used to perform numerical simulation on the minimum value and the maximum value in the critical detonation probability interval respectively using a preset coupling algorithm to generate simulation results;

[0226] The adjustment module 204 is used to adjust the model parameters of the initial model of the ignition growth of high-energy viscoelastic materials based on the simulation results to obtain the ignition growth model of high-energy viscoelastic materials.

[0227] Furthermore, in some embodiments, the construction module 201 may specifically include:

[0228] A first calculation unit is used to calculate the number of hot spots and the ignition probability in the hot spot formation and ignition stage;

[0229] A second calculation unit is used to calculate the number of hot spots and the reaction area in the hot spot growth stage;

[0230] The third calculation unit is used to calculate the energy data of the reaction completion stage.

[0231] Furthermore, in some embodiments, the test module 202 may be specifically used to:

[0232] respectively setting corresponding partitions of different thicknesses for a plurality of high-energy viscoelastic materials of the same size and shape;

[0233] Conduct impact tests on each high-energy viscoelastic material using a standardized impact source to generate impact test results corresponding to partitions of different thicknesses;

[0234] Calculating the detonation probability at each partition thickness based on the impact test results;

[0235] Based on the impact test results and detonation probabilities corresponding to partitions of different thicknesses, a critical detonation probability interval of the high-energy viscoelastic material is determined.

[0236] Furthermore, in some embodiments, the test module 202 may be specifically used for:

[0237] Constructing a partition test geometric model, wherein the partition test geometric model includes a main propellant, a partition material, a propellant column, a steel casing and a witness plate material;

[0238] A partition test material model is constructed, wherein the partition test material model includes a main charge material model, a partition material model, a steel casing material model and a witness plate material model.

[0239] Furthermore, in some embodiments, the simulation module 203 may be specifically used for:

[0240] Performing numerical simulation on a partition model corresponding to a minimum value in the critical detonation probability interval to generate a first simulation result in the simulation results;

[0241] A numerical simulation is performed on a partition model corresponding to the maximum value in the critical detonation probability interval to generate a second simulation result in the simulation results.

[0242] Furthermore, in some embodiments, the preset coupling algorithm is a smooth particle fluid dynamics-finite element coupling algorithm, and the simulation results include a numerical simulation pressure cloud map and pressure time history curves corresponding to multiple observation points.

[0243] Furthermore, in some embodiments, the construction device may further include a verification module, and the verification module is specifically used to:

[0244] Calculate the corresponding theoretical analysis results through the pre-built shock wave parameter calculation model;

[0245] Obtaining test results corresponding to the diaphragm test;

[0246] Compare the experimental test results, the simulation results and the theoretical analysis results to obtain corresponding comparison results;

[0247] Verifying the ignition growth model of the high-energy viscoelastic material based on the comparison result to generate a verification result;

[0248] The ignition growth model of the high-energy viscoelastic material is adjusted secondary based on the verification results.

[0249] As described above, an ignition growth model construction device for high-energy viscoelastic materials provided in an embodiment of the present application adds a variety of data collection means to obtain a variety of measured data based on a small partition test, and adjusts the model parameters in combination with numerical calculation and parameter fitting methods to improve the accuracy and applicability of the constructed ignition growth model, providing a physical model and numerical algorithm basis for the research on impact initiation of high-energy viscoelastic materials, and is more economical, simple and quicker than traditional methods.

[0250] The process of each module in the high-energy viscoelastic material ignition growth model construction device provided in this embodiment realizing its own function can be specifically referred to the aforementioned Figure 1 The description of the illustrated embodiment and other related method embodiments will not be repeated here.

[0251] It should be noted that the information interaction, execution process, etc. between the above-mentioned devices / units are based on the same concept as the method embodiment of the present application. Their specific functions and technical effects can be found in the method embodiment part and will not be repeated here.

[0252] It should be noted that in the embodiments of the present application, if the above test method is implemented in the form of a software function module and sold or used as an independent product, it can also be stored in a computer-readable storage medium. Based on this understanding, the technical solution of the embodiment of the present application can be essentially or partly reflected in the form of a software product that contributes to the prior art. The computer software product is stored in a storage medium, including several instructions to enable a computer device (which can be a personal computer, a server, or a network device, etc.) to execute all or part of the methods described in each embodiment of the present application. The aforementioned storage medium includes: various media that can store program codes, such as a USB flash drive, a mobile hard disk, a read-only memory (ROM, Read Only Memory), a magnetic disk or an optical disk. In this way, the embodiments of the present application are not limited to any specific combination of hardware and software.

[0253] Correspondingly, an embodiment of the present application further discloses a storage medium storing a computer program that can be loaded by a processor and execute the above-mentioned method for constructing an ignition growth model of a high-energy viscoelastic material.

[0254] The present application also discloses a computer device, such as Fig.15 As shown, it includes a processor 100, at least one communication bus 200, a user interface 300, at least one external communication interface 400 and a memory 500. The processor 100 is used to provide computing and control capabilities, the communication bus 200 is configured to achieve connection and communication between these components, the user interface 300 may include a display screen, the external communication interface 400 may include a standard wired interface and a wireless interface, and the memory 500 stores a method for rapid cold impact test of solid power devices. The processor 100 is also used to execute the steps of the method for constructing an ignition growth model of a high-energy viscoelastic material stored in the memory 500.

[0255] Optionally, in some embodiments, the processor implements the following steps when executing the computer program:

[0256] Construct an initial model of the ignition growth of high-energy viscoelastic materials; obtain the critical detonation probability interval of high-energy viscoelastic materials through partition tests; use a preset coupling algorithm to numerically simulate the minimum and maximum values ​​in the critical detonation probability interval to generate simulation results; adjust the model parameters of the initial model of the ignition growth of high-energy viscoelastic materials based on the simulation results to obtain the ignition growth model of high-energy viscoelastic materials.

[0257] Optionally, in some embodiments, a computer-readable storage medium is provided, wherein the computer-readable storage medium stores a computer program, and when the computer program is executed by a processor, the steps of the method for constructing an ignition growth model of a high-energy viscoelastic material are implemented in any one of the above-mentioned steps.

[0258] This embodiment adds a variety of data collection methods based on the small partition test to obtain a variety of measured data, combines numerical calculation and parameter fitting methods to adjust the model parameters, improves the accuracy and applicability of the constructed ignition growth model, and provides a physical model and numerical algorithm basis for the research on impact initiation of high-energy viscoelastic materials. Compared with traditional methods, it is more economical, simple and quick.

[0259] It should be understood that "one embodiment" or "some embodiments" mentioned throughout the specification means that specific features, structures or characteristics related to the embodiment are included in at least one embodiment of the present application. Therefore, "in one embodiment" or "in an embodiment" appearing throughout the specification does not necessarily refer to the same embodiment. In addition, these specific features, structures or characteristics can be combined in one or more embodiments in any suitable manner. It should be understood that in various embodiments of the present application, the size of the sequence number of the above-mentioned processes does not mean the order of execution. The execution order of each process should be determined by its function and internal logic, and should not constitute any limitation on the implementation process of the embodiments of the present application. The above-mentioned sequence numbers of the embodiments of the present application are for description only and do not represent the advantages and disadvantages of the embodiments.

[0260] It should be noted that, in this article, the terms "include", "comprises" or any other variations thereof are intended to cover non-exclusive inclusion, so that a process, method, article or device including a series of elements includes not only those elements, but also other elements not explicitly listed, or also includes elements inherent to such process, method, article or device. In the absence of further restrictions, an element defined by the sentence "comprises a ..." does not exclude the existence of other identical elements in the process, method, article or device including the element.

[0261] In the several embodiments provided in the present application, it should be understood that the disclosed devices and methods can be implemented in other ways. The device embodiments described above are only schematic. For example, the division of the units is only a logical function division. There may be other division methods in actual implementation, such as: multiple units or components can be combined, or can be integrated into another system, or some features can be ignored or not executed. In addition, the coupling, direct coupling, or communication connection between the components shown or discussed can be through some interfaces, and the indirect coupling or communication connection of the devices or units can be electrical, mechanical or other forms.

[0262] The units described above as separate components may or may not be physically separated, and the components displayed as units may or may not be physical units; they may be located in one place or distributed on multiple network units; some or all of the units may be selected according to actual needs to achieve the purpose of the present embodiment.

[0263] In addition, all functional units in the embodiments of the present application may be integrated into one processing unit, or each unit may be a separate unit, or two or more units may be integrated into one unit; the above-mentioned integrated units may be implemented in the form of hardware or in the form of hardware plus software functional units.

[0264] A person of ordinary skill in the art can understand that: all or part of the steps of implementing the above-mentioned method embodiment can be completed by hardware related to program instructions, and the aforementioned program can be stored in a computer-readable storage medium. When the program is executed, it executes the steps of the above-mentioned method embodiment; and the aforementioned storage medium includes: various media that can store program codes, such as mobile storage devices, ROMs, magnetic disks or optical disks.

[0265] Alternatively, if the above-mentioned integrated unit of the present application is implemented in the form of a software function module and sold or used as an independent product, it can also be stored in a computer-readable storage medium. Based on this understanding, the technical solution of the embodiment of the present application can essentially or in other words, the part that contributes to the prior art can be embodied in the form of a software product, which is stored in a storage medium and includes several instructions for a device to execute all or part of the methods described in each embodiment of the present application. The aforementioned storage medium includes: various media that can store program codes, such as mobile storage devices, ROMs, magnetic disks or optical disks.

[0266] The above disclosure is only the preferred embodiment of the present application, which certainly cannot be used to limit the scope of rights of the present application. Therefore, equivalent changes made according to the claims of the present application are still within the scope covered by the present application.

Claims

1. A method for constructing an ignition growth model of a high-energy viscoelastic material, characterized in that: include: Construct an initial model for the ignition growth of high-energy viscoelastic materials; The critical detonation probability interval of the high-energy viscoelastic material is obtained through a partition test, including: constructing a partition test geometric model, the partition test geometric model includes a main propellant, a partition material, a propellant column, a steel casing and a witness plate material; constructing a partition test material model, the partition test material model includes a main propellant column material model, a partition material model, a steel casing material model and a witness plate material model; based on the partition test geometric model and the partition test material model, the critical detonation probability interval of the high-energy viscoelastic material is simulated; the lower limit of the critical detonation probability interval is the minimum partition thickness at which the viscoelastic material detonates 100%, and the upper limit is the maximum partition thickness at which the viscoelastic material does not detonate; A preset coupling algorithm is used to numerically simulate the minimum and maximum values ​​in the critical detonation probability interval, generate simulation results, and fit the ignition item parameters in the ignition growth initial model one by one to ensure that the calculation process and calculation results of the numerical simulation are consistent with the experiment, identify the main parameters affecting the impact sensitivity of the high-energy viscoelastic material, and determine the accurate values ​​of the main parameters; the preset coupling algorithm is a smooth particle fluid dynamics-finite element coupling algorithm, and the simulation results include a numerical simulation pressure cloud map and a pressure time history curve corresponding to multiple observation points; Based on the simulation results, the model parameters of the initial model of ignition growth of high-energy viscoelastic material are adjusted to obtain the ignition growth model of high-energy viscoelastic material.

2. The method for constructing an ignition growth model of a high-energy viscoelastic material according to claim 1, characterized in that: The construction of the initial model of ignition growth of high-energy viscoelastic materials comprises: Calculating the number of hot spots and the ignition probability in the hot spot formation and ignition stages corresponding to the high-energy viscoelastic material to construct the first stage of the ignition growth initial model; Calculating the number of hot spots and the reaction area in the hot spot growth stage to construct the second stage of the ignition growth initial model; The energy data of the reaction completion stage is calculated to construct the third stage of the ignition growth initial model.

3. The method for constructing an ignition growth model of a high-energy viscoelastic material according to claim 1, characterized in that: The method of obtaining the critical detonation probability interval of the high-energy viscoelastic material through the partition test includes: respectively setting corresponding partitions of different thicknesses for a plurality of high-energy viscoelastic materials of the same size and shape; Conduct impact tests on each high-energy viscoelastic material using a standardized impact source to generate impact test results corresponding to partitions of different thicknesses; Calculating the detonation probability at each partition thickness based on the impact test results; Based on the impact test results and detonation probabilities corresponding to partitions of different thicknesses, a critical detonation probability interval of the high-energy viscoelastic material is determined.

4. The method for constructing an ignition growth model of a high-energy viscoelastic material according to claim 1, characterized in that: The method of using a preset coupling algorithm to perform numerical simulation on the minimum value and the maximum value in the critical detonation probability interval to generate simulation results includes: Performing numerical simulation on a partition model corresponding to a minimum value in the critical detonation probability interval to generate a first simulation result in the simulation results, wherein the first simulation result is used to simulate an experimental result of the partition model corresponding to the minimum value; A numerical simulation is performed on the partition model corresponding to the maximum value in the critical detonation probability interval to generate a second simulation result in the simulation results, wherein the second simulation result is used to simulate the experimental result of the partition model corresponding to the maximum value.

5. The method for constructing an ignition growth model of a high-energy viscoelastic material according to claim 1, characterized in that: After adjusting the model parameters of the initial model of the ignition growth of the high-energy viscoelastic material based on the simulation results to obtain the ignition growth model of the high-energy viscoelastic material, the method further includes: Calculate the corresponding theoretical analysis results through the pre-built shock wave parameter calculation model; Obtaining test results corresponding to the diaphragm test; Compare the experimental test results, the simulation results and the theoretical analysis results to obtain corresponding comparison results; Verifying the ignition growth model of the high-energy viscoelastic material based on the comparison result to generate a verification result; The ignition growth model of the high-energy viscoelastic material is adjusted secondary based on the verification results.

6. A device for constructing an ignition growth model of a high-energy viscoelastic material, characterized in that: include: Building module for constructing the initial model of ignition growth of high-energy viscoelastic materials; A test module, for obtaining a critical detonation probability interval of a high-energy viscoelastic material through a partition test, comprising: constructing a partition test geometric model, the partition test geometric model comprising a main propellant, a partition material, a propellant column, a steel casing and a witness plate material; constructing a partition test material model, the partition test material model comprising a main propellant column material model, a partition material model, a steel casing material model and a witness plate material model; simulating a critical detonation probability interval of a high-energy viscoelastic material based on the partition test geometric model and the partition test material model; the lower limit of the critical detonation probability interval is a minimum partition thickness at which 100% of the viscoelastic material detonates, and the upper limit is a maximum partition thickness at which no detonation occurs in the viscoelastic material; A simulation module, used to use a preset coupling algorithm to numerically simulate the minimum and maximum values ​​in the critical detonation probability interval, generate simulation results, and fit the ignition item parameters in the ignition growth initial model one by one to ensure that the calculation process and calculation results of the numerical simulation are consistent with the experiment, identify the main parameters affecting the impact sensitivity of the high-energy viscoelastic material, and determine the accurate values ​​of the main parameters; the preset coupling algorithm is a smooth particle fluid dynamics-finite element coupling algorithm, and the simulation results include a numerical simulation pressure cloud map and a pressure time history curve corresponding to multiple observation points; The adjustment module is used to adjust the model parameters of the initial model of the ignition growth of the high-energy viscoelastic material based on the simulation results to obtain the ignition growth model of the high-energy viscoelastic material.

7. A computer device comprising a memory, a processor, and a computer program stored in the memory and executable on the processor, characterized in that: When the processor executes the computer program, the steps of the method for constructing an ignition growth model of a high-energy viscoelastic material as claimed in any one of claims 1 to 5 are implemented.

8. A computer-readable storage medium storing a computer program, characterized in that: When the computer program is executed by a processor, the steps of the method for constructing an ignition growth model of a high-energy viscoelastic material as claimed in any one of claims 1 to 5 are implemented.

Images