Construction method, device and equipment of viscous-elastic material detonation model and medium

Abstract

The invention discloses a viscous-elastic material detonation model construction method and device, equipment and a medium, and the method comprises the steps: obtaining formula data and density data of a target viscous-elastic material, and obtaining at least two detonation velocity data of the target viscous-elastic material during a small partition experiment; determining a detonation reaction equation of the target viscoelastic material according to the formula data, and determining a local isentropic index of the target viscoelastic material; constructing an initial detonation model according to a JWL equation; obtaining a parameter group of a JWL equation according to the density data, the local isentropic index and each detonation velocity data; target detonation velocity data are selected from all the detonation velocity data, and target parameters are selected from the parameter set; and updating the initial detonation model based on the target detonation velocity data and the target parameters to obtain a detonation model of the target viscous-elastic material. The detonation model of the viscous-elastic material is constructed on the basis of measured data of a small partition plate test in combination with numerical calculation, and the method is more economical, simpler, more convenient and quicker.

Description

Technical Field

[0001] The invention relates to the technical field of viscoelastic materials, and in particular to a method, device, equipment and medium for constructing a detonation model of a viscoelastic material. Background Art

[0002] The detonation product state equation can describe the relationship between various physical quantities in the detonation products of high-energy viscoelastic materials, reflect the detonation ability of high-energy viscoelastic materials, and is also the basis for solving explosion mechanics problems by numerical means. The construction process of the viscoelastic material detonation model can be understood as the parameter calculation process of the detonation product state equation. Due to the wide application of the Jones-Wilkins-Lee (JWL) state equation, the viscoelastic material detonation model in the prior art often adopts the JWL state equation, that is, the construction process of the viscoelastic material detonation model can be understood as the parameter calculation process of the JWL state equation. The Jones-Wilkins-Lee (JWL) state equation is a semi-empirical state equation proposed by Lee based on the research of Jones and Wilkins, which does not explicitly contain chemical reactions and has parameters determined by experimental methods. It can accurately describe the expansion-driven work process of the detonation products.

[0003] In the prior art, when calculating the parameters of the JWL state equation, experimental data of a specific experiment (cylinder test) and two-dimensional fluid dynamics calculations are required. However, the experimental cost of the specific experiment (cylinder test) is high, and the experimental process is cumbersome, which is not conducive to the construction of the detonation model of viscoelastic materials. At the same time, the calculation process of two-dimensional fluid dynamics is cumbersome and requires a large computing power, which is not conducive to the construction of the detonation model of viscoelastic materials. Therefore, when constructing the detonation model of viscoelastic materials, the efficiency is low and the process is complicated. Summary of the invention

[0004] Based on this, it is necessary to propose a method, device, equipment and medium for constructing a viscoelastic material detonation model in order to address the technical problems of low efficiency and complicated process in the existing technology when constructing a viscoelastic material detonation model.

[0005] In a first aspect, a method for constructing a detonation model of a viscoelastic material is provided, the method comprising:

[0006] Obtaining formula data and density data of a target viscoelastic material, and obtaining at least two detonation velocity data of the target viscoelastic material when performing a small partition experiment;

[0007] Determining a detonation reaction equation of the target viscoelastic material according to the formula data and the Binkley-Wilson method, and determining a local isentropic index of the target viscoelastic material based on the detonation reaction equation;

[0008] Construct the initial detonation model according to the JWL equation;

[0009] Substituting the density data, the local isentropic index and each detonation velocity data into the initial detonation model to obtain a parameter group of the JWL equation;

[0010] Selecting target detonation velocity data from among the detonation velocity data based on the actual pressure data, and selecting a target parameter from among the parameter group;

[0011] The initial detonation model is updated based on the target detonation velocity data and the target parameters to obtain a detonation model of the target viscoelastic material.

[0012] Optionally, the step of determining the local isentropic index of the target viscoelastic material based on the detonation reaction equation comprises:

[0013] Determine various detonation products of the target viscoelastic material and the proportion data of each detonation product based on the detonation reaction equation;

[0014] Determining a sub-local isentropic index of each detonation product according to preset local isentropic index data;

[0015] The local isentropic index of the target viscoelastic material is calculated by using the sub-local isentropic index of each detonation product and the proportion data of each detonation product.

[0016] Optionally, the step of calculating the local isentropic index of the target viscoelastic material by using the sub-local isentropic index of each detonation product and the proportion data of each detonation product comprises:

[0017] Calculate the product of the sub-local isentropic index and the proportion data of the same detonation product to obtain the sub-value of the detonation product;

[0018] The sub-values ​​of the detonation products are accumulated to obtain the local isentropic index of the target viscoelastic material.

[0019] Optionally, the step of constructing an initial detonation model according to the JWL equation comprises:

[0020] Determine the relationship equations of various physical quantities of the target viscoelastic material according to the JWL equation to obtain a relationship equation group;

[0021] Each equation in the relationship equation group is combined to obtain the initial detonation model.

[0022] Optionally, the step of substituting the density data, the local isentropic index and each detonation velocity data into the initial detonation model to obtain the parameter group of the JWL equation includes:

[0023] Substituting the density data, the local isentropic index and each detonation velocity data into the initial detonation model to obtain a γ-law state equation;

[0024] Draw the isentropic expansion curve of the viscoelastic material under different detonation velocity data according to the γ-law state equation;

[0025] Fitting the isentropic expansion curve and the detonation velocity data corresponding to the isentropic expansion curve to obtain parameters of the JWL equation under different detonation velocity data;

[0026] The parameters of the JWL equation under different detonation velocity data are summarized to form the parameter group.

[0027] Optionally, the step of selecting target detonation velocity data from various detonation velocity data based on actual pressure data comprises:

[0028] A preselected JWL equation is established according to the parameters of the JWL equation under different detonation velocity data to obtain a preselected JWL equation group;

[0029] Based on each preselected JWL equation in the preselected JWL equation group, a SPH-FEM coupling algorithm is used to simulate an internal explosion test to obtain a simulation data group;

[0030] Acquiring actual pressure data of the target viscoelastic material in an actual explosion experiment, and matching the actual pressure data with each simulation data in the simulation data group;

[0031] using simulation data matching the actual pressure data as target simulation data;

[0032] The detonation velocity data corresponding to the target simulation data among the detonation velocity data is used as the target detonation velocity data.

[0033] Optionally, the step of selecting a target parameter in the parameter group includes:

[0034] The preselected JWL equation corresponding to the target simulation data in the preselected JWL equation group is used as a target JWL equation, and the parameters corresponding to the target JWL equation in the parameter group are used as target parameters.

[0035] In a second aspect, a device for constructing a detonation model of a viscoelastic material is provided, the device comprising:

[0036] A data acquisition unit, used to acquire formula data and density data of a target viscoelastic material, and to acquire at least two detonation velocity data of the target viscoelastic material when performing a small partition experiment;

[0037] A first calculation unit, configured to determine a detonation reaction equation of the target viscoelastic material according to the formula data and the Binkley-Wilson method, and determine a local isentropic index of the target viscoelastic material based on the detonation reaction equation;

[0038] A model building unit, used for building an initial detonation model according to the JWL equation;

[0039] A second calculation unit is used to substitute the density data, the local isentropic index and each detonation velocity data into the initial detonation model to obtain a parameter group of the JWL equation;

[0040] a detonation velocity selection unit, for selecting target detonation velocity data from among the detonation velocity data based on the actual pressure data, and selecting a target parameter from among the parameter group;

[0041] A model updating unit is used to update the initial detonation model based on the target detonation velocity data and the target parameters to obtain a detonation model of the target viscoelastic material.

[0042] In a third aspect, an electronic device is provided, comprising a memory, a processor, and a computer program stored in the memory and executable on the processor, wherein the processor implements the steps of the method for constructing the above-mentioned viscoelastic material detonation model when executing the computer program.

[0043] In a fourth aspect, 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 a viscoelastic material detonation model are implemented.

[0044] The present application obtains the formula data and density data of the target viscoelastic material, and obtains at least two detonation velocity data of the target viscoelastic material when conducting a small partition experiment; determines the detonation reaction equation of the target viscoelastic material according to the formula data and the Binkley-Wilson method, and determines the local isentropic index of the target viscoelastic material based on the detonation reaction equation; constructs an initial detonation model according to the JWL equation; substitutes the density data, the local isentropic index and each detonation velocity data into the initial detonation model to obtain a parameter group of the JWL equation; selects target detonation velocity data from each detonation velocity data based on actual pressure data, and selects target parameters from the parameter group; updates the initial detonation model based on the target detonation velocity data and the target parameters to obtain a detonation model of the target viscoelastic material. By obtaining the measured data of the small partition test and combining it with numerical calculation (calculation of the JWL equation parameters), a detonation model of the viscoelastic material is constructed. Compared with the specific test (cylinder test) in the prior art, the small partition test is more economical, simpler and faster. Compared with the calculation process of the two-dimensional fluid dynamics calculation in the prior art, the numerical calculation (calculation of the JWL equation parameters) has fewer calculation steps and parameters involved in the calculation, and the calculation process is more convenient, which improves the efficiency of constructing the detonation model of the target viscoelastic material. BRIEF DESCRIPTION OF THE DRAWINGS

[0045] In order to more clearly illustrate the embodiments of the present invention 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 invention. For ordinary technicians in this field, other drawings can be obtained based on these drawings without paying creative work.

[0046] in:

[0047] Figure 1 is an application environment diagram of a method for constructing a viscoelastic material detonation model in one embodiment;

[0048] Figure 2 is a flow chart of a method for constructing a detonation model of a viscoelastic material in one embodiment;

[0049] Figure 3 is a structural block diagram of a device for constructing a detonation model of a viscoelastic material in one embodiment;

[0050] Figure 4 is a structural block diagram of a computer device in one embodiment;

[0051] Figure 5 It is a structural block diagram of a computer device in another embodiment. DETAILED DESCRIPTION

[0052] The following will be combined with the drawings in the embodiments of the present invention to clearly and completely describe the technical solutions in the embodiments of the present invention. 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.

[0053] The method for constructing a viscoelastic material detonation model provided by the embodiment of the present invention can be applied in the following aspects: Figure 1In an application environment, an experimental data measuring device (sensor) 110 communicates with a server 120, the server 120 obtains the formula data and density data of a target viscoelastic material, and obtains at least two detonation velocity data of the target viscoelastic material during a small partition experiment through the experimental data measuring device (sensor) 110; a detonation reaction equation of the target viscoelastic material is determined according to the formula data and the Binkley-Wilson method, and a local isentropic index of the target viscoelastic material is determined based on the detonation reaction equation; an initial detonation model is constructed according to a JWL equation; the density data, the local isentropic index and each detonation velocity data are substituted into the initial detonation model to obtain a parameter group of the JWL equation; target detonation velocity data is selected from each detonation velocity data based on actual pressure data, and a target parameter is selected from the parameter group; the initial detonation model is updated based on the target detonation velocity data and the target parameter to obtain a detonation model of the target viscoelastic material. By obtaining the measured data of the small partition test and combining it with numerical calculation (calculation of the JWL equation parameters), a detonation model of the viscoelastic material is constructed. Compared with the specific test (cylinder test) in the prior art, the small partition test is more economical, simpler and faster. Compared with the calculation process of the two-dimensional fluid dynamics calculation in the prior art, the numerical calculation (calculation of the JWL equation parameters) has fewer calculation steps and parameters involved in the calculation, and the calculation process is more convenient, which improves the efficiency of constructing the detonation model of the target viscoelastic material.

[0054] The present invention is described in detail below through specific embodiments.

[0055] See also Figure 2 As shown, Figure 2 A flow chart of a method for constructing a viscoelastic material detonation model provided in an embodiment of the present invention includes the following steps:

[0056] S101, obtaining formula data and density data of a target viscoelastic material, and obtaining at least two detonation velocity data of the target viscoelastic material when performing a small partition experiment;

[0057] Exemplarily, the formula data includes the components of the target viscoelastic material and the proportion of each component. For example, the target viscoelastic material per unit volume includes component A and component B, component A accounts for 40%, and component B accounts for 60%. The detonation velocity data is the speed at which the detonation wave propagates stably in the explosive, and is an important indicative number for measuring the explosive performance of the explosive. The detonation velocity data can be estimated using the Martin formula, the Kamlet formula, or can be measured using the Dautriche method, the timing method, and the optics method.

[0058] Exemplarily, the small diaphragm experiment accurately measures the shock wave pressure of NEPE (nitrate plasticized polyether) propellant detonation products at different positions in the internal explosion model tunnel, and obtains at least two detonation velocity data of the target viscoelastic material when performing the small diaphragm experiment through at least two groups of small diaphragm experiments.

[0059] S102, determining a detonation reaction equation of the target viscoelastic material according to the formula data and the Binkley-Wilson method, and determining a local isentropic index of the target viscoelastic material based on the detonation reaction equation;

[0060] Exemplarily, the chemical formula of the target viscoelastic material is determined according to the formula data, wherein when the proportion of a component of the target viscoelastic material is less than 5% and substantially does not participate in the exothermic reaction, the component is ignored.

[0061] Exemplarily, the detonation reaction equation of the target viscoelastic material is determined according to the chemical equation of the target viscoelastic material and the Binkley-Wilson method.

[0062] Exemplarily, the Binkley-Wilson method includes detonation products of various components of the target viscoelastic material after detonation, for example, aluminum powder of the target viscoelastic material undergoes oxidation reaction during detonation to produce aluminum oxide.

[0063] S103, constructing an initial detonation model according to the JWL equation;

[0064] Exemplarily, the Jones-Wilkins-Lee (JWL) state equation is a semi-empirical state equation proposed by Lee based on the research of Jones and Wilkins. It does not explicitly contain chemical reactions and has parameters determined by experimental methods. It can accurately describe the expansion-driven work process of detonation products.

[0065] S104, substituting the density data, the local isentropic index and each detonation velocity data into the initial detonation model to obtain a parameter group of the JWL equation;

[0066] For example, if the fluid does not exchange heat with the outside world during the state change (from one state to another), the process is called an adiabatic process. If the adiabatic process is a reversible process, then in this process, the product of the pressure P of the fluid and the κ-th power of the specific volume V is a constant, and κ is the local isentropic index.

[0067] Exemplarily, at least two detonation velocity data are obtained from the small partition experiment. It is obviously unscientific to select a single detonation velocity data or calculate the average value of multiple detonation velocity data to perform JWL equation parameter fitting. Therefore, each detonation velocity data is substituted into the initial detonation model respectively to obtain the parameters of at least two JWL equations, and then the JWL equation parameters are fitted according to the parameters of at least two JWL equations.

[0068] Exemplarily, after substituting each detonation velocity data into the initial detonation model, different JWL equations can be obtained, and different JWL equations have different parameters. The parameters of each JWL equation are summarized to obtain the parameter group. For example, after substituting the detonation velocity data A into the initial detonation model, the JWL equation N can be obtained, and the parameter of the JWL equation N is Q. After substituting the detonation velocity data B into the initial detonation model, the JWL equation M can be obtained, and the parameter of the JWL equation N is P. Then the parameters of each JWL equation are summarized to obtain the parameter group (Q, P).

[0069] S105, selecting target detonation velocity data from each detonation velocity data based on the actual pressure data, and selecting a target parameter from the parameter group;

[0070] Exemplarily, the actual pressure data is pressure data obtained through an actual detonation experiment, for example, actual pressure data of a target point in an actual detonation experiment is measured, and initial pressure data of the target point under different detonation velocity data are calculated based on the JWL equation under different detonation velocity data to obtain an initial pressure data set;

[0071] The initial pressure data with the smallest error with the actual pressure data is selected from the initial pressure data set as the target pressure data, the detonation velocity data corresponding to the target pressure data is used as the target detonation velocity data, and the target detonation velocity data is substituted into the initial detonation model to obtain the model to be updated.

[0072] S106: Update the initial detonation model based on the target detonation velocity data and the target parameters to obtain a detonation model of the target viscoelastic material.

[0073] The present application obtains the formula data and density data of the target viscoelastic material, and obtains at least two detonation velocity data of the target viscoelastic material when conducting a small partition experiment; determines the detonation reaction equation of the target viscoelastic material according to the formula data and the Binkley-Wilson method, and determines the local isentropic index of the target viscoelastic material based on the detonation reaction equation; constructs an initial detonation model according to the JWL equation; substitutes the density data, the local isentropic index and each detonation velocity data into the initial detonation model to obtain the parameters of the JWL equation; selects the target detonation velocity data from each detonation velocity data based on the actual pressure data, and substitutes the target detonation velocity data into the initial detonation model to obtain the model to be updated; updates the model to be updated based on the parameters to obtain the detonation model of the target viscoelastic material. By obtaining the measured data of the small partition test and combining it with numerical calculation (calculation of the JWL equation parameters), a detonation model of the viscoelastic material is constructed. Compared with the specific test (cylinder test) in the prior art, the small partition test is more economical, simpler and faster. Compared with the calculation process of the two-dimensional fluid dynamics calculation in the prior art, the numerical calculation (calculation of the JWL equation parameters) has fewer calculation steps and parameters involved in the calculation, and the calculation process is more convenient, which improves the efficiency of constructing the detonation model of the target viscoelastic material.

[0074] In a possible implementation, the step of determining the local isentropic index of the target viscoelastic material based on the detonation reaction equation includes:

[0075] Determine various detonation products of the target viscoelastic material and the proportion data of each detonation product based on the detonation reaction equation;

[0076] Determining a sub-local isentropic index of each detonation product according to preset local isentropic index data;

[0077] The local isentropic index of the target viscoelastic material is calculated by using the sub-local isentropic index of each detonation product and the proportion data of each detonation product.

[0078] For example, taking 1 kg of nitrate ester plasticized polyether (NEPE) propellant as an example, when the mass percentage of aluminum powder is 5% to 30%, only 33% to 38% of the aluminum powder participates in the reaction, the detonation reaction equation is:

[0079] 0.15N 8.53 H 34.13 Cl 8.53 O 34.13 +0.17Al 37.05 +0.42C 13.52 H 27.05 N 27.05 O 27.05 +0.18C 16.67 H29.17 N 12.5 O 37.5 →0.76Al 2 O 3 +1.28HCl+10.22H 2 O+6.63CO+2.04CO 2 +7.45N 2

[0080] Among them, N 8.53 The number of nitrogen atoms is 8.53, H 34.13 The number of atoms representing hydrogen is 34.13, Cl 8.53 The number of atoms representing chlorine is 8.53, Al 37.05 The number of atoms representing aluminum is 37.05, O 34.13 The number of atoms representing oxygen is 34.13, C 13.52 The number of atoms representing the element carbon is 13.52.

[0081] Exemplarily, according to the detonation reaction equation, it is determined that the detonation product of the target viscoelastic material includes aluminum oxide (Al 2 O 3 ), hydrogen chloride (HCl), water (H 2 O), carbon monoxide (CO), carbon dioxide (CO 2 ) and nitrogen (N 2 ), among which aluminum oxide (Al 2 O 3 ) accounted for 0.76, hydrogen chloride (HCl) accounted for 1.28, and water (H 2 O) accounted for 10.22, carbon monoxide (CO) accounted for 6.63, and carbon dioxide (CO 2 ) accounted for 2.04, nitrogen (N 2 ) accounts for 7.45.

[0082] Exemplarily, the preset local isentropic index data includes the sub-local isentropic index of each detonation product, for example, aluminum oxide (Al 2 O 3 ) has a sub-local isentropic index of 1.145, that of hydrogen chloride (HCl) is 1.39, and that of water (H 2 O) is 1.9, carbon monoxide (CO) is 2.85, and carbon dioxide (CO 2 ) has a sub-local isentropic index of 4.5, and nitrogen (N 2 ) has a sub-local isentropic index of 3.70.

[0083] In a possible implementation, the step of calculating the local isentropic index of the target viscoelastic material by using the sub-local isentropic index of each detonation product and the proportion data of each detonation product includes:

[0084] Calculate the product of the sub-local isentropic index and the proportion data of the same detonation product to obtain the sub-value of the detonation product;

[0085] The sub-values ​​of the detonation products are accumulated to obtain the local isentropic index of the target viscoelastic material.

[0086] Exemplarily, the local isentropic index calculation formula is:

[0087]

[0088] The product of the sub-local isentropic index of each detonation product and the proportion data is calculated to obtain the sub-value of each detonation product, and the sub-values ​​of each detonation product are accumulated to obtain the local isentropic index of the target viscoelastic material, wherein γ is the local isentropic index of the target viscoelastic material, γ i is the sublocal isentropic index of the ith detonation product, x i is the amount of substance in the i-th detonation product, and ∑ is the calculation symbol for cumulative summation.

[0089] In a possible implementation manner, the step of constructing an initial detonation model according to the JWL equation includes:

[0090] Determine the relationship equations of various physical quantities of the target viscoelastic material according to the JWL equation to obtain a relationship equation group;

[0091] Each equation in the relationship equation group is combined to obtain the initial detonation model.

[0092] Exemplarily, the JWL equation is:

[0093]

[0094] Where p is the pressure, V = V e / V 0 =ρ 0 / ρ e is the relative volume, V e is the volume of the detonation products, ρ e is the density of the detonation product, E is the specific internal energy, A, B, R 1 , R 2 , ω is a constant, ρ 0 is the initial density of the target viscoelastic material, V 0 is the initial volume of the target viscoelastic material, and e is an exponential constant.

[0095] The JWL equation can be simplified to:

[0096]

[0097] Where C is the pressure constant.

[0098] Exemplarily, based on the simplified JWL equation, the relationship equations of various physical quantities of the target viscoelastic material are determined to obtain a relationship equation group:

[0099] p=Av -γ =Aρ γ

[0100]

[0101] Wherein, ρ is the density of the target viscoelastic material after detonation, ρ 0 is the initial density of the target viscoelastic material, v is the volume of the target viscoelastic material after detonation, and v 0 is the initial volume of the target viscoelastic material, D is the detonation velocity data, A is the pressure constant, and γ is the local isentropic index of the target viscoelastic material.

[0102] Exemplarily, each equation in the relationship equation group is combined to obtain the initial detonation model:

[0103]

[0104] In a possible implementation manner, the step of substituting the density data, the local isentropic index and each detonation velocity data into the initial detonation model to obtain the parameters of the JWL equation includes:

[0105] Substituting the density data, the local isentropic index and each detonation velocity data into the initial detonation model to obtain a γ-law state equation;

[0106] Draw the isentropic expansion curve of the viscoelastic material under different detonation velocity data according to the γ-law state equation;

[0107] Fitting the isentropic expansion curve and the detonation velocity data corresponding to the isentropic expansion curve to obtain parameters of the JWL equation under different detonation velocity data;

[0108] The parameters of the JWL equation under different detonation velocity data are summarized to form the parameter group.

[0109] Exemplarily, the density data, the local isentropic index and each detonation velocity data are substituted into the initial detonation model, and four detonation velocity values: 6866.7 m / s, 6876.7 m / s, 7191.8 m / s and 7357.9 m / s are taken as examples to obtain the γ-law state equation:

[0110]

[0111] Exemplarily, an isentropic expansion curve (PV curve) of a high-energy viscoelastic material is plotted according to the γ-law state equation, and the isentropic expansion curve and detonation velocity data corresponding to the isentropic expansion curve are fitted to obtain parameters of the JWL equation under different detonation velocity data.

[0112] In a possible implementation, the step of selecting target detonation velocity data from various detonation velocity data based on actual pressure data includes:

[0113] A preselected JWL equation is established according to the parameters of the JWL equation under different detonation velocity data to obtain a preselected JWL equation group;

[0114] Based on each preselected JWL equation in the preselected JWL equation group, a SPH-FEM coupling algorithm is used to simulate an internal explosion test to obtain a simulation data group;

[0115] Acquiring actual pressure data of the target viscoelastic material in an actual explosion experiment, and matching the actual pressure data with each simulation data in the simulation data group;

[0116] using simulation data matching the actual pressure data as target simulation data;

[0117] The detonation velocity data corresponding to the target simulation data among the detonation velocity data is used as the target detonation velocity data.

[0118] Exemplarily, 6866.7m / s, 6876.7m / s, 7191.8m / s and 7357.9m / s are respectively substituted into the initial detonation model to obtain four preselected JWL equations (preselected JWL equation group), and the SPH-FEM coupling algorithm is used to simulate the internal explosion test based on each preselected JWL equation to obtain a simulation data group. The finite element-smoothed particle fluid dynamics coupling (SPH-FEM coupling) algorithm is an algorithm in which the control quantity and the controlled quantity influence and are interrelated with each other, and a simulation data group is obtained.

[0119] Exemplarily, the actual pressure data can be understood as the pressure data actually exerted on the target point in an actual explosion experiment. The simulated pressure of the target point is predicted based on each preselected JWL equation to obtain a simulation data group. The actual pressure data is compared with the simulated pressure in the simulation data group, and the simulated pressure (simulated data) with the smallest error with the actual pressure data is selected as the target simulation data, and the detonation velocity data corresponding to the target simulation data is used as the target detonation velocity data.

[0120] In a possible implementation manner, the step of selecting target detonation velocity data from various detonation velocity data based on actual pressure data, and substituting the target detonation velocity data into the initial detonation model to obtain the model to be updated includes:

[0121] The preselected JWL equation corresponding to the target simulation data in the preselected JWL equation group is used as a target JWL equation, and the parameters corresponding to the target JWL equation in the parameter group are used as target parameters.

[0122] Exemplarily, 6866.7m / s, 6876.7m / s, 7191.8m / s and 7357.9m / s are respectively substituted into the initial detonation model to obtain four preselected JWL equations (preselected JWL equation group), and the SPH-FEM coupling algorithm is used to simulate the internal explosion test based on each preselected JWL equation to obtain a simulation data group. The finite element-smoothed particle fluid dynamics coupling (SPH-FEM coupling) algorithm is an algorithm in which the control quantity and the controlled quantity influence and are interrelated with each other, and a simulation data group is obtained.

[0123] Exemplarily, the actual pressure data is used as a criterion, and more accurate JWL equation parameters are selected as target parameters in each simulation data.

[0124] For example, the simulation results of the internal explosion test simulated by the SPH-FEM coupling algorithm are as follows: At t = 10 μs, the propellant detonation products begin to expand outward. At t = 50 μs, the shock wave is transmitted to the tank wall and the pressure relief port. At t = 100 μs, the shock wave at the pressure relief port begins to propagate outward, and the shock wave at the explosion tank wall is reflected inward due to contact with the wall. At t = 200 μs, due to the reflection of the shock wave, an irregular local high-pressure area appears in the explosion tank, and the particles propagate outward along the pressure relief pipe.

[0125] Exemplarily, the relative errors between the simulated values ​​of the parameters of each pre-selected JWL equation in the pre-selected JWL equation group and the experimental data are shown in Table 1:

[0126] Table 1: Relative error between parameter simulation values ​​and test data

[0127] Serial number D(m / s) P(Mpa) Relative error 1 6866.7 4.45 -6.5％ 2 6876.7 4.46 -6.3％ 3 7191.8 4.62 -2.9％ 4 7357.9 4.79 0.63％

[0128] It can be seen that the JWL equation obtained by fitting with D = 7357.9 m / s can most accurately reflect the state of detonation products of the high-energy viscoelastic material used in the experiment, with a peak pressure error of 0.63%. That is, 7357.9 m / s is used as the target detonation velocity data, and the parameters of the JWL equation obtained by fitting with D = 7357.9 m / s are used as the target parameters.

[0129] See also Figure 3As shown, in one embodiment, a device for constructing a viscoelastic material detonation model is provided, the device comprising:

[0130] The data acquisition unit 201 is used to acquire the formula data and density data of the target viscoelastic material, and acquire at least two detonation velocity data of the target viscoelastic material when performing a small partition experiment;

[0131] A first calculation unit 202 is used to determine a detonation reaction equation of the target viscoelastic material according to the formula data and the Binkley-Wilson method, and determine a local isentropic index of the target viscoelastic material based on the detonation reaction equation;

[0132] A model building unit 203 is used to build an initial detonation model according to the JWL equation;

[0133] A second calculation unit 204 is used to substitute the density data, the local isentropic index and each detonation velocity data into the initial detonation model to obtain a parameter group of the JWL equation;

[0134] A detonation velocity selection unit 205, for selecting target detonation velocity data from various detonation velocity data based on actual pressure data, and selecting target parameters from the parameter group;

[0135] The model updating unit 206 is used to update the initial detonation model based on the target detonation velocity data and the target parameters to obtain the detonation model of the target viscoelastic material.

[0136] In a possible implementation manner, the first calculation unit 202 is configured to:

[0137] Determine various detonation products of the target viscoelastic material and the proportion data of each detonation product based on the detonation reaction equation;

[0138] Determining a sub-local isentropic index of each detonation product according to preset local isentropic index data;

[0139] The local isentropic index of the target viscoelastic material is calculated by using the sub-local isentropic index of each detonation product and the proportion data of each detonation product.

[0140] In a possible implementation manner, the first calculation unit 202 is configured to:

[0141] Calculate the product of the sub-local isentropic index and the proportion data of the same detonation product to obtain the sub-value of the detonation product;

[0142] The sub-values ​​of the detonation products are accumulated to obtain the local isentropic index of the target viscoelastic material.

[0143] In a possible implementation, the model building unit 203 is used to:

[0144] Determine the relationship equations of various physical quantities of the target viscoelastic material according to the JWL equation to obtain a relationship equation group;

[0145] Each equation in the relationship equation group is combined to obtain the initial detonation model.

[0146] In a possible implementation manner, the second calculation unit 204 is configured to:

[0147] Substituting the density data, the local isentropic index and each detonation velocity data into the initial detonation model to obtain a γ-law state equation;

[0148] Draw the isentropic expansion curve of the viscoelastic material under different detonation velocity data according to the γ-law state equation;

[0149] Fitting the isentropic expansion curve and the detonation velocity data corresponding to the isentropic expansion curve to obtain parameters of the JWL equation under different detonation velocity data;

[0150] The parameters of the JWL equation under different detonation velocity data are summarized to form the parameter group.

[0151] In a possible implementation manner, the detonation velocity selection unit 205 is used to:

[0152] A preselected JWL equation is established according to the parameters of the JWL equation under different detonation velocity data to obtain a preselected JWL equation group;

[0153] Based on each preselected JWL equation in the preselected JWL equation group, a SPH-FEM coupling algorithm is used to simulate an internal explosion test to obtain a simulation data group;

[0154] Acquiring actual pressure data of the target viscoelastic material in an actual explosion experiment, and matching the actual pressure data with each simulation data in the simulation data group;

[0155] using simulation data matching the actual pressure data as target simulation data;

[0156] The detonation velocity data corresponding to the target simulation data among the detonation velocity data is used as the target detonation velocity data.

[0157] In a possible implementation manner, the detonation velocity selection unit 205 is used to:

[0158] The preselected JWL equation corresponding to the target simulation data in the preselected JWL equation group is used as a target JWL equation, and the parameters corresponding to the target JWL equation in the parameter group are used as target parameters.

[0159] The present application obtains the formula data and density data of the target viscoelastic material, and obtains at least two detonation velocity data of the target viscoelastic material when conducting a small partition experiment; determines the detonation reaction equation of the target viscoelastic material according to the formula data and the Binkley-Wilson method, and determines the local isentropic index of the target viscoelastic material based on the detonation reaction equation; constructs an initial detonation model according to the JWL equation; substitutes the density data, the local isentropic index and each detonation velocity data into the initial detonation model to obtain a parameter group of the JWL equation; selects target detonation velocity data from each detonation velocity data based on actual pressure data, and selects target parameters from the parameter group; updates the initial detonation model based on the target detonation velocity data and the target parameters to obtain a detonation model of the target viscoelastic material. By obtaining the measured data of the small partition test and combining it with numerical calculation (calculation of the JWL equation parameters), a detonation model of the viscoelastic material is constructed. Compared with the specific test (cylinder test) in the prior art, the small partition test is more economical, simpler and faster. Compared with the calculation process of the two-dimensional fluid dynamics calculation in the prior art, the numerical calculation (calculation of the JWL equation parameters) has fewer calculation steps and parameters involved in the calculation, and the calculation process is more convenient, which improves the efficiency of constructing the detonation model of the target viscoelastic material.

[0160] In one embodiment, a computer device is provided. The computer device may be a server, and its internal structure diagram may be as follows: Figure 4 As shown. The computer device includes a processor, a memory, a network interface and a database connected via a system bus. Among them, the processor of the computer device is used to provide computing and control capabilities. The memory of the computer device includes a non-volatile and / or volatile storage medium and an internal memory. The non-volatile storage medium stores an operating system, a computer program and a database. The internal memory provides an environment for the operation of the operating system and the computer program in the non-volatile storage medium. The network interface of the computer device is used to communicate with an external client via a network connection. When the computer program is executed by the processor, the functions or steps on the server side of a method for constructing a viscoelastic material detonation model are realized.

[0161] In one embodiment, a computer device is provided. The computer device may be a client, and its internal structure diagram may be as follows: Figure 5As shown. The computer device includes a processor, a memory, a network interface, a display screen and an input device connected through a system bus. Among them, the processor of the computer device is used to provide computing and control capabilities. The memory of the computer device includes a non-volatile storage medium and an internal memory. The non-volatile storage medium stores an operating system and a computer program. The internal memory provides an environment for the operation of the operating system and the computer program in the non-volatile storage medium. The network interface of the computer device is used to communicate with an external server through a network connection. When the computer program is executed by the processor, the functions or steps on the client side of a method for constructing a viscoelastic material detonation model are realized.

[0162] In one embodiment, the computer device can be used as one of the electronic devices. An electronic device is proposed, including a memory, a processor, and a computer program stored in the memory and executable on the processor. When the processor executes the computer program, the following steps are implemented: obtaining the formula data and density data of the target viscoelastic material, and obtaining at least two detonation velocity data of the target viscoelastic material when conducting a small partition experiment; determining the detonation reaction equation of the target viscoelastic material according to the formula data and the Binkley-Wilson method, and determining the local isentropic index of the target viscoelastic material based on the detonation reaction equation; constructing an initial detonation model according to the JWL equation; substituting the density data, the local isentropic index and each detonation velocity data into the initial detonation model to obtain a parameter group of the JWL equation; selecting target detonation velocity data from each detonation velocity data based on actual pressure data, and selecting target parameters from the parameter group; updating the initial detonation model based on the target detonation velocity data and the target parameters to obtain a detonation model of the target viscoelastic material.

[0163] The present application obtains the formula data and density data of the target viscoelastic material, and obtains at least two detonation velocity data of the target viscoelastic material when conducting a small partition experiment; determines the detonation reaction equation of the target viscoelastic material according to the formula data and the Binkley-Wilson method, and determines the local isentropic index of the target viscoelastic material based on the detonation reaction equation; constructs an initial detonation model according to the JWL equation; substitutes the density data, the local isentropic index and each detonation velocity data into the initial detonation model to obtain a parameter group of the JWL equation; selects target detonation velocity data from each detonation velocity data based on actual pressure data, and selects target parameters from the parameter group; updates the initial detonation model based on the target detonation velocity data and the target parameters to obtain a detonation model of the target viscoelastic material. By obtaining the measured data of the small partition test and combining it with numerical calculation (calculation of the JWL equation parameters), a detonation model of the viscoelastic material is constructed. Compared with the specific test (cylinder test) in the prior art, the small partition test is more economical, simpler and faster. Compared with the calculation process of the two-dimensional fluid dynamics calculation in the prior art, the numerical calculation (calculation of the JWL equation parameters) has fewer calculation steps and parameters involved in the calculation, and the calculation process is more convenient, which improves the efficiency of constructing the detonation model of the target viscoelastic material.

[0164] In one embodiment, a computer-readable storage medium is proposed, wherein the computer-readable storage medium stores a computer program, and when the computer program is executed by a processor, the following steps are implemented: obtaining formula data and density data of a target viscoelastic material, and obtaining at least two detonation velocity data of the target viscoelastic material when conducting a small partition experiment; determining a detonation reaction equation of the target viscoelastic material according to the formula data and the Binkley-Wilson method, and determining a local isentropic index of the target viscoelastic material based on the detonation reaction equation; constructing an initial detonation model according to a JWL equation; substituting the density data, the local isentropic index and each detonation velocity data into the initial detonation model to obtain a parameter group of the JWL equation; selecting target detonation velocity data from each detonation velocity data based on actual pressure data, and selecting target parameters from the parameter group; updating the initial detonation model based on the target detonation velocity data and the target parameters to obtain a detonation model of the target viscoelastic material.

[0165] The present application obtains the formula data and density data of the target viscoelastic material, and obtains at least two detonation velocity data of the target viscoelastic material when conducting a small partition experiment; determines the detonation reaction equation of the target viscoelastic material according to the formula data and the Binkley-Wilson method, and determines the local isentropic index of the target viscoelastic material based on the detonation reaction equation; constructs an initial detonation model according to the JWL equation; substitutes the density data, the local isentropic index and each detonation velocity data into the initial detonation model to obtain a parameter group of the JWL equation; selects target detonation velocity data from each detonation velocity data based on actual pressure data, and selects target parameters from the parameter group; updates the initial detonation model based on the target detonation velocity data and the target parameters to obtain a detonation model of the target viscoelastic material. By obtaining the measured data of the small partition test and combining it with numerical calculation (calculation of the JWL equation parameters), a detonation model of the viscoelastic material is constructed. Compared with the specific test (cylinder test) in the prior art, the small partition test is more economical, simpler and faster. Compared with the calculation process of the two-dimensional fluid dynamics calculation in the prior art, the numerical calculation (calculation of the JWL equation parameters) has fewer calculation steps and parameters involved in the calculation, and the calculation process is more convenient, which improves the efficiency of constructing the detonation model of the target viscoelastic material.

[0166] It should be noted that the above functions or steps that can be implemented by the computer-readable storage medium or computer device can refer to the relevant descriptions on the server side and the client side in the aforementioned method embodiment. To avoid repetition, they will not be described one by one here.

[0167] Those skilled in the art can understand that all or part of the processes in the above-mentioned embodiment methods can be completed by instructing the relevant hardware through a computer program, and the computer program can be stored in a non-volatile computer-readable storage medium. When the computer program is executed, it can include the processes of the embodiments of the above-mentioned methods. Among them, any reference to memory, storage, database or other media used in the embodiments provided in this application can include non-volatile and / or volatile memory. Non-volatile memory can include read-only memory (ROM), programmable ROM (PROM), electrically programmable ROM (EPROM), electrically erasable programmable ROM (EEPROM) or flash memory. Volatile memory can include random access memory (RAM) or external cache memory. As an illustration and not limitation, RAM is available in many forms, such as static RAM (SRAM), dynamic RAM (DRAM), synchronous DRAM (SDRAM), double data rate SDRAM (DDRSDRAM), enhanced SDRAM (ESDRAM), synchronous link (Synchlink) DRAM (SLDRAM), memory bus (Rambus) direct RAM (RDRAM), direct memory bus dynamic RAM (DRDRAM), and memory bus dynamic RAM (RDRAM).

[0168] Those skilled in the art can clearly understand that for the convenience and simplicity of description, only the division of the above-mentioned functional units and modules is used as an example. In actual applications, the above-mentioned functions can be distributed and completed by different functional units and modules as needed, that is, the internal structure of the device can be divided into different functional units or modules to complete all or part of the functions described above.

[0169] The embodiments described above are only used to illustrate the technical solutions of the present invention, rather than to limit the same. Although the present invention has been described in detail with reference to the aforementioned embodiments, those skilled in the art should understand that the technical solutions described in the aforementioned embodiments may still be modified, or some of the technical features may be replaced by equivalents. Such modifications or replacements do not deviate the essence of the corresponding technical solutions from the spirit and scope of the technical solutions of the embodiments of the present invention, and should all be included in the protection scope of the present invention.

Claims

1. A method for constructing a detonation model of a viscoelastic material, characterized in that: The method comprises: Obtaining formula data and density data of a target viscoelastic material, and obtaining at least two detonation velocity data of the target viscoelastic material when performing a small partition experiment; Determining a detonation reaction equation of the target viscoelastic material according to the formula data and the Binkley-Wilson method, and determining a local isentropic index of the target viscoelastic material based on the detonation reaction equation; Construct the initial detonation model according to the JWL equation; Substituting the density data, the local isentropic index and each detonation velocity data into the initial detonation model to obtain a parameter group of the JWL equation; Selecting target detonation velocity data from among the detonation velocity data based on the actual pressure data, and selecting a target parameter from among the parameter group; The initial detonation model is updated based on the target detonation velocity data and the target parameters to obtain a detonation model of the target viscoelastic material.

2. The method for constructing a viscoelastic material detonation model according to claim 1, characterized in that: The step of determining the local isentropic index of the target viscoelastic material based on the detonation reaction equation comprises: Determine various detonation products of the target viscoelastic material and the proportion data of each detonation product based on the detonation reaction equation; Determining a sub-local isentropic index of each detonation product according to preset local isentropic index data; The local isentropic index of the target viscoelastic material is calculated by using the sub-local isentropic index of each detonation product and the proportion data of each detonation product.

3. The method for constructing a viscoelastic material detonation model according to claim 2, characterized in that: The step of calculating the local isentropic index of the target viscoelastic material by using the sub-local isentropic index of each detonation product and the proportion data of each detonation product comprises: Calculate the product of the sub-local isentropic index and the proportion data of the same detonation product to obtain the sub-value of the detonation product; The sub-values ​​of the detonation products are accumulated to obtain the local isentropic index of the target viscoelastic material.

4. The method for constructing a viscoelastic material detonation model according to claim 1, characterized in that: The step of constructing the initial detonation model according to the JWL equation comprises: Determine the relationship equations of various physical quantities of the target viscoelastic material according to the JWL equation to obtain a relationship equation group; Each equation in the relationship equation group is combined to obtain the initial detonation model.

5. The method for constructing a viscoelastic material detonation model according to claim 1, characterized in that: The step of substituting the density data, the local isentropic index and each detonation velocity data into the initial detonation model to obtain the parameter group of the JWL equation includes: Substituting the density data, the local isentropic index and each detonation velocity data into the initial detonation model to obtain a γ-law state equation; Draw the isentropic expansion curve of the target viscoelastic material under different detonation velocity data according to the γ-law state equation; Fitting the isentropic expansion curve and the detonation velocity data corresponding to the isentropic expansion curve to obtain parameters of the JWL equation under different detonation velocity data; The parameters of the JWL equation under different detonation velocity data are summarized to form the parameter group.

6. The method for constructing a viscoelastic material detonation model according to claim 5, characterized in that: The step of selecting target detonation velocity data from various detonation velocity data based on actual pressure data comprises: A preselected JWL equation is established according to the parameters of the JWL equation under different detonation velocity data to obtain a preselected JWL equation group; Based on each preselected JWL equation in the preselected JWL equation group, a SPH-FEM coupling algorithm is used to simulate an internal explosion test to obtain a simulation data group; Acquiring actual pressure data of the target viscoelastic material in an actual explosion experiment, and matching the actual pressure data with each simulation data in the simulation data group; using simulation data matching the actual pressure data as target simulation data; The detonation velocity data corresponding to the target simulation data among the detonation velocity data is used as the target detonation velocity data.

7. The method for constructing a viscoelastic material detonation model according to claim 6, characterized in that: The step of selecting a target parameter in the parameter group comprises: The preselected JWL equation corresponding to the target simulation data in the preselected JWL equation group is used as a target JWL equation, and the parameters corresponding to the target JWL equation in the parameter group are used as target parameters.

8. A device for constructing a detonation model of a viscoelastic material, characterized in that: The device comprises: A data acquisition unit, used to acquire formula data and density data of a target viscoelastic material, and to acquire at least two detonation velocity data of the target viscoelastic material when performing a small partition experiment; A first calculation unit, configured to determine a detonation reaction equation of the target viscoelastic material according to the formula data and the Binkley-Wilson method, and determine a local isentropic index of the target viscoelastic material based on the detonation reaction equation; A model building unit, used for building an initial detonation model according to the JWL equation; A second calculation unit is used to substitute the density data, the local isentropic index and each detonation velocity data into the initial detonation model to obtain a parameter group of the JWL equation; a detonation velocity selection unit, for selecting target detonation velocity data from among the detonation velocity data based on the actual pressure data, and selecting a target parameter from among the parameter group; A model updating unit is used to update the initial detonation model based on the target detonation velocity data and the target parameters to obtain a detonation model of the target viscoelastic material.

9. An electronic 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 a viscoelastic material detonation model according to any one of claims 1 to 7 are implemented.

10. 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 a viscoelastic material detonation model as claimed in any one of claims 1 to 7 are implemented.

Images