Pressure state equation determination method and device, simulation equipment and computer storage medium

Abstract

The invention relates to a pressure state equation determination method and device, simulation equipment and a computer storage medium, and belongs to the technical field of explosion and impact dynamics, and the state equation determination method comprises the steps: taking real experiment data as the input of a numerical iteration model, and outputting to obtain initial simulation axial stress; adjusting the actually-measured axial stress based on the initial simulation axial stress to obtain an adjusted actually-measured axial stress, taking the actually-measured volume strain and the adjusted actually-measured axial stress as the input of a numerical iteration model, and outputting to obtain a transition simulation axial stress; and when the actually measured axial stress does not need to be adjusted, determining the current transition simulation axial stress as the final axial stress, and constructing a wide pressure range state equation of the experimental object by taking the final axial stress as the hydrostatic pressure. According to the method, the obtained simulated axial stress can serve as the real hydrostatic pressure of the experimental object, and then a wide-range pressure state equation of the experimental object can be accurately constructed.

Description

Technical Field

[0001] The present invention relates to the technical field of explosion and shock dynamics, and particularly to a method, device, simulation equipment and computer storage medium for determining a pressure equation of state. Background Art

[0002] In the high-precision numerical simulation of the dynamic damage and failure of concrete-like materials under explosion and shock loads, the commonly used model for concrete-like materials under explosion and shock loads is the fluid elastoplastic model. Under explosion and shock loads, the dynamic tensile damage and failure of concrete-like materials occur under low hydrostatic pressures, the dynamic shear failure occurs under relatively low or medium hydrostatic pressures, and the dynamic compressive damage and failure occur under high hydrostatic pressures. Therefore, the material model must be able to better describe the dynamic mechanical and damage behaviors of concrete-like materials under low, medium, and high pressures at the same time.

[0003] To quantitatively analyze various damage effects in the near-middle-far zones generated by explosion and shock loads and the stress wave propagation and attenuation laws, it is necessary to establish a pressure-density (volumetric strain) relationship curve that can completely describe the pressure range experienced by concrete-like materials during the entire load response process, that is, the equation of state for concrete-like materials in a wide pressure range. Limited by equipment and technology, the maximum hydrostatic pressure applied in the current hydrostatic compression test is below 1 GPa, which does not yet meet the high-pressure state of concrete-like materials in the near zone under explosion and shock loads at the order of 10 GPa. Therefore, the plate impact experiment is used to replace the hydrostatic compression test to obtain a higher hydrostatic pressure. However, in the plate impact experiment, if the concrete-like material cannot be regarded as a non-viscous compressible fluid, the existence of shear strength effects and strain rate effects will cause the measured axial stress to be much higher than the true hydrostatic pressure, resulting in a large error and making it impossible to accurately construct the equation of state for concrete-like materials in a wide range.

[0004] It can be seen that if we want to accurately construct the equation of state for concrete-like materials in a wide pressure range under explosion loads, it is necessary to use an appropriate method to eliminate the shear strength effects and strain rate effects in the plate impact experiment and correct the measured axial stress in the plate impact experiment to the true hydrostatic pressure required for constructing the equation of state. Summary of the Invention

[0005] In view of this, it is necessary to provide a method, device, simulation equipment and computer storage medium for determining a pressure equation of state to solve the problem that the true hydrostatic pressure cannot be obtained due to the inability to eliminate the shear strength effects and strain rate effects in the existing plate impact experiment, resulting in the inability to accurately construct the equation of state for concrete-like materials in a wide range.

[0006] To solve the above problems, in the first aspect, the present invention provides a method for determining a pressure equation of state, including: Taking the real experimental data as the input of the numerical iteration model, the initial simulated axial stress is output, and the real experimental data are the measured axial stress and the measured volume strain of the experimental object in the real plate impact experiment; Adjusting the measured axial stress based on the initial simulated axial stress to obtain the adjusted measured axial stress, and taking the measured volume strain and the adjusted measured axial stress as the input of the numerical iteration model, the transitional simulated axial stress is output; When the measured axial stress does not need to be adjusted, the current transitional simulated axial stress is determined as the final axial stress, and the final axial stress is used as the hydrostatic pressure to construct the state equation of the wide pressure range of the experimental object.

[0007] In a possible implementation manner, the real experimental data include multiple measured axial stresses and the measured volume strains corresponding to the measured axial stresses, and the measured axial stresses cover the wide pressure range borne by the experimental object under the explosion shock load.

[0008] In a possible implementation manner, taking the real experimental data as the input of the numerical iteration model and outputting the initial simulated axial stress includes: Taking the real test data as the input of the state equation representing the axial stress-strain volume relationship of the experimental object in the numerical iteration model, and adjusting the impact velocity of the flyer impacting the experimental object in the numerical iteration model to obtain the simulated axial stresses corresponding to different measured axial stresses.

[0009] In a possible implementation manner, the adjustment process of multiple measured axial stresses includes: Adjusting the measured axial stresses in sequence according to the magnitude order of the impact velocities corresponding to the measured axial stresses.

[0010] In a possible implementation manner, adjusting the measured axial stress based on the initial simulated axial stress to obtain the adjusted measured axial stress includes: Calculating the first difference between the measured axial stress and the initial simulated axial stress, and calculating the absolute value of the ratio of the first difference to the measured axial stress; When the absolute value of the ratio is greater than or equal to a preset threshold, taking the second difference between the measured axial stress and the first difference as the adjusted measured axial stress.

[0011] In a possible implementation manner, the condition that the measured axial stress does not need to be adjusted is that the absolute value of the ratio is less than the preset threshold.

[0012] In a possible implementation manner, taking the simulated axial stress that meets the preset condition as the hydrostatic pressure to construct the state equation of the wide pressure range of the experimental object includes: Construct the hydrostatic pressure - volume strain relationship of the experimental object by taking the simulated axial stress that meets the preset conditions as the hydrostatic pressure; Construct the state equation of the experimental object within a wide pressure range according to the hydrostatic pressure - volume strain relationship.

[0013] In a second aspect, the present invention further provides a state equation determination device, including: An initial simulated axial stress determination module, configured to take the real experimental data as the input of a numerical iteration model and output the initial simulated axial stress, where the real experimental data is the measured axial stress and the measured volume strain obtained by the experimental object in a real plate impact experiment; A transition simulated axial stress determination module, configured to adjust the measured axial stress based on the initial simulated axial stress to obtain an adjusted measured axial stress, take the measured volume strain and the adjusted measured axial stress as the input of the numerical iteration model, and output the transition simulated axial stress; A pressure state equation construction module, configured to, when the measured axial stress does not need to be adjusted, determine the current transition simulated axial stress as the final axial stress, and take the final axial stress as the hydrostatic pressure to construct the state equation of the experimental object within a wide pressure range.

[0014] In a third aspect, the present invention further provides a simulation device, including a memory and a processor, wherein, The memory is used to store programs; The processor is coupled to the memory and is configured to execute the programs stored in the memory to implement the steps in the state equation determination method according to any one of the above - mentioned embodiments.

[0015] In a fourth aspect, the present invention further provides a computer - readable storage medium, configured to store computer - readable programs or instructions, and when the programs or instructions are executed by a processor, the steps in the state equation determination method according to any one of the above - mentioned embodiments can be implemented.

[0016] The beneficial effects of the present invention are as follows: The method for determining the state equation provided by the present invention takes real test data as the input of the numerical iteration model, outputs the initial simulated axial stress, ensures the reliability of the simulated experimental data, adjusts the measured axial stress based on the initial simulated axial stress to obtain the adjusted measured axial stress, and takes the measured volume strain and the adjusted measured axial stress as the input of the numerical iteration model, and outputs the transitional simulated axial stress, which can overcome the shear strength effect and strain rate effect in the plate impact experiment, prevent the repeated simulation of the shear strength effect and strain rate effect in the simulated plate impact experiment, and further can use the obtained simulated axial stress as the true hydrostatic pressure of the experimental object, ensure the accuracy of the numerical simulation of the true hydrostatic pressure, and further accurately construct the pressure state equation of the experimental object in a wide range. BRIEF DESCRIPTION OF THE DRAWINGS

[0017] In order to more clearly illustrate the technical solutions in the embodiments of the present invention, the following will briefly introduce the drawings required for the description of the embodiments. Obviously, the following drawings are only some embodiments of the present invention. For those skilled in the art, without creative efforts, other drawings can be obtained based on these drawings.

[0018] Figure 1 It is a schematic flowchart of a method for determining a pressure state equation provided by an embodiment of the present invention; Figure 2 It is a numerical simulation structure diagram of a plate impact experiment provided by an embodiment of the present invention; Figure 3 It is a schematic flowchart of iterative correction of simulated axial pressure provided by an embodiment of the present invention; Figure 4 It is a data diagram of three simulated plate impact experiments provided by an embodiment of the present invention; Figure 5 It is a schematic flowchart of a method for determining a state equation in a wide pressure range provided by an embodiment of the present invention; Figure 6 It is a verification comparison diagram provided by an embodiment of the present invention; Figure 7 It is a schematic structural diagram of a state equation determination device provided by an embodiment of the present invention; Figure 8 It is a schematic structural diagram of a simulation device provided by an embodiment of the present invention. DETAILED DESCRIPTION OF THE EMBODIMENTS

[0019] The following will specifically describe the preferred embodiments of the present invention in conjunction with the drawings. Among them, the drawings form a part of this application and are used together with the embodiments of the present invention to explain the principle of the present invention, and are not used to limit the scope of the present invention.

[0020] In the embodiments of the present invention, the descriptions such as "first" and "second" are only for descriptive purposes, and cannot be construed as indicating or implying their relative importance or implicitly specifying the quantity of the indicated technical features. Therefore, the technical features defined with "first" and "second" may explicitly or implicitly include at least one such feature.

[0021] Reference to "embodiment" in this text means that a particular feature, structure, or characteristic described in connection with the embodiment can be included in at least one embodiment of the present invention. The phrase appears in various places in the specification and does not necessarily refer to the same embodiment, nor is it an independent or alternative embodiment mutually exclusive with other embodiments. Those skilled in the art will explicitly and implicitly understand that the embodiments described herein can be combined with other embodiments.

[0022] A specific embodiment of the present invention, as Figure 1 shown, discloses a method for determining a pressure state equation, including: S101, taking real experimental data as the input of a numerical iteration model, and outputting an initial simulated axial stress, where the real experimental data is the measured axial stress and the measured volume strain of the experimental object in a real plate impact experiment.

[0023] In the embodiments of the present invention, the experimental object is a concrete-like material, such as concrete, rock, cement mortar, etc. Specifically, taking limestone as an example, a real plate impact experiment is carried out on limestone, and the real experimental data is the relevant data of the experimental object obtained in the real experiment, including the measured axial stress and the measured volume strain.

[0024] In some possible embodiments of the present invention, the experimental data includes a plurality of measured axial stresses and the measured volume strains corresponding to each measured axial stress, and the measured axial stresses cover a wide pressure range that the experimental object bears under explosive shock loads.

[0025] In the embodiments of the present invention, in order to ensure the accuracy of the experimental data and to ensure that the real experimental data can cover the wide pressure range from megapascals (MPa) to gigapascals (GPa) that the concrete-like material bears when subjected to explosive shock loads, it is necessary to conduct real plate impact experiments multiple times to obtain a plurality of measured axial stresses and the measured volume strains corresponding to each measured axial stress. As shown in Table 1, a corresponding relationship table of the measured axial stress and the measured volume strain can be obtained: Table 1: Corresponding relationship table of measured axial stress and measured volume strain

[0026] In the embodiments of the present invention, to implement the numerical simulation of the plate impact experiment, as Figure 2As shown, a numerical model for numerical iteration can be established using the Smoothed Particle Galerkin (SPG) algorithm, where the target is the fluid elastoplastic model of the experimental object. The test data in Table 1 above is used as the input of the tabular equation of state in the fluid elastoplastic model to simulate the plate impact experiment. Among them, the tabular equation of state is used to represent the axial stress-strain volume relationship that the experimental object bears in the plate impact experiment. By adjusting the impact velocity when the flyer impacts the target, different simulated axial stresses are obtained. The tabular equation of state refers to the equation of state that shows the corresponding relationship between input and output in a tabular manner. This equation of state does not give or does not need to give a specific function expression, and can give the single-point output value corresponding to the equation of state for a single-point input value.

[0027] S102. Adjust the measured axial stress based on the initial simulated axial stress to obtain the adjusted measured axial stress. Use the measured volume strain and the adjusted measured axial stress as the input of the numerical iteration model, and output to obtain the transitional simulated axial stress.

[0028] In the embodiment of the present invention, due to the coupled influence of the shear strength effect and the strain rate effect in the simulated plate impact experiment, the simulated axial stress output by the plate impact experiment is greater than the measured axial stress. It is necessary to correct the measured axial stress that is the input of the tabular equation of state, and use the corrected adjusted measured axial stress and the measured volume strain as the input of the new tabular equation of state, and iteratively correct the measured axial stress until the simulated axial stress in the simulated plate impact experiment meets the preset conditions. Among them, the correction method for the measured axial stress will be described in detail later in the present invention.

[0029] S103. When the measured axial stress does not need to be adjusted, determine the current transitional simulated axial stress as the final axial stress, and use the final axial stress as the hydrostatic pressure to construct the wide pressure range equation of state of the experimental object.

[0030] In the embodiment of the present invention, when the measured axial stress does not need to be adjusted, because the shear strength effect and the strain rate effect are eliminated in the numerical simulation process of the plate impact experiment, the simulated axial stress can be directly used as the hydrostatic pressure of the experimental object. Therefore, the accurate hydrostatic pressure of the experimental object is obtained, and the wide range pressure equation of the experimental object can be constructed based on this hydrostatic pressure.

[0031] The method for determining the equation of state provided by the present invention takes real experimental data as the input of a numerical iteration model, outputs the initial simulated axial stress, ensures the reliability of the simulated experimental data, adjusts the measured axial stress based on the initial simulated axial stress to obtain the adjusted measured axial stress, and takes the measured volumetric strain and the adjusted measured axial stress as the input of the numerical iteration model, and outputs the transitional simulated axial stress, which can overcome the shear strength effect and strain rate effect in the plate impact experiment, prevent the repeated simulation of the shear strength effect and strain rate effect in the simulated plate impact experiment, and then the obtained simulated axial stress can be used as the true hydrostatic pressure of the experimental object, ensuring the accuracy of the numerical simulation of the true hydrostatic pressure, and further accurately constructing the pressure equation of state for a wide range of experimental objects.

[0032] In some possible embodiments of the present invention, the adjustment process of multiple measured axial stresses includes: Adjust the measured axial stresses in sequence according to the magnitude order of the impact velocities corresponding to the measured axial stresses.

[0033] In the embodiments of the present invention, since there are multiple measured axial stresses, to ensure the accuracy of the final hydrostatic pressure - volumetric strain relationship curve, when correcting the measured axial stresses, it is necessary to correct them in sequence. Only after the correction of the previous measured axial stress is completed, can the next measured axial stress be corrected. The specific correction sequence is the magnitude order of the flyer impact velocities corresponding to each measured axial stress. As shown in Table 1, the correction should start from the bottom - most measured axial stress in the table and proceed upwards in sequence until all the measured axial stresses are corrected.

[0034] The embodiments of the present invention correct each measured axial stress in sequence, ensure the accuracy of the relationship between each measured axial stress and the volumetric strain, and finally can determine a more accurate hydrostatic pressure - volumetric strain relationship curve.

[0035] In some possible embodiments of the present invention, as Figure 3 shown, the adjustment of the measured axial stress based on the initial simulated axial stress to obtain the adjusted measured axial stress includes: S301, calculate the first difference between the measured axial stress and the initial simulated axial stress, and calculate the absolute value of the ratio of the first difference to the measured axial stress; S302, when the absolute value of the ratio is greater than or equal to a preset threshold, use the second difference between the measured axial stress and the first difference as the adjusted measured axial stress.

[0036] In the embodiments of the present invention, first, the measured axial stress and the measured strain volume (where \(i = 1, 2,\cdots, n\), representing the \(i\)-th input data point in the tabular equation of state) is used as the first input to the tabular equation of state. For the sake of illustration, taking a specific embodiment as an example, taking the numerical simulation of the plate impact experiment with a flyer impact velocity of 1280 m / s in Table 1 as an example, the measured axial stress of the first input is i , the volumetric strain is , and the axial stress output from the simulated plate impact experiment is , . Due to the coupled influence of the shear strength effect and the strain rate, must be greater than . At this time, it is necessary to calculate , and then according to correct to obtain . replaces the measured axial stress with a flyer impact velocity of 1280 m / s in the measured axial stress, and then and are used as the second input to the tabular equation of state to conduct the second simulated plate impact experiment with a flyer impact velocity of 1280 m / s. The simulated axial stress output from the second simulated plate impact experiment is . Generally, at this time, still has a large difference from , and the difference between the two is . Similarly, use to correct to obtain , and then use and as the third input to the tabular equation of state, conduct the simulated plate impact experiment again, and judge whether holds. When holds, is the true hydrostatic pressure with a flyer impact velocity of 1280 m / s, where \(j\) is the \(j\)-th simulated plate impact experiment. is the difference between the axial stress and the hydrostatic pressure caused by the coupled influence of the shear strength effect and the strain rate effect in the plate impact experiment with a flyer impact velocity of 1280 m / s. As Figure 4 shows, it is a data graph of three simulated plate impact experiments. Based on this principle, the measured axial stress can be corrected multiple times.

[0037] In the embodiment of the present invention, the measured axial stress is corrected by an iterative strategy, which can effectively remove the difference between the axial stress and the hydrostatic pressure caused by the coupled influence of the shear strength effect and the strain rate effect, and thus obtain an accurate hydrostatic pressure.

[0038] In some possible embodiments of the present invention, such as Figure 5 shown, a wide pressure range equation of state for the experimental object is constructed by taking the simulated axial stress that meets the preset conditions as the hydrostatic pressure, including: S501, constructing the hydrostatic pressure - volume strain relationship of the experimental object by taking the simulated axial stress that meets the preset conditions as the hydrostatic pressure; S502, constructing the wide pressure range equation of state of the experimental object according to the hydrostatic pressure - volume strain relationship.

[0039] In the embodiments of the present invention, after obtaining the hydrostatic pressure data based on the iterative strategy in the foregoing embodiments, the hydrostatic pressure - volume strain relationship is constructed. Since in the case where the measured axial stress is much higher than the concrete shear strength due to the high flyer impact velocity, the material can be regarded as a non - viscous compressible fluid. At this time, the shear strength effect and strain rate effect can be ignored, so there is no need for numerical iteration. Therefore, for the simulated axial stress in the high - pressure state, it can be not corrected by iteration. The medium - low pressure equation of state of concrete after removing the shear strength effect and strain rate effect and the high - pressure equation of state considering concrete - like materials as non - viscous compressible fluids are combined to establish the wide pressure range equation of state of concrete - like materials.

[0040] Furthermore, based on the established wide pressure range equation of state of concrete - like materials, combined with the strength surface parameters calibrated by triaxial compression test data and the strain rate effect parameters calibrated by Hopkinson bar test data, a high - precision limestone fluid elastoplastic model that can better describe the dynamic mechanical and damage behavior of concrete under low, medium, and high pressures is established. Using the existing experimental data of the penetration depth of limestone targets at different penetration velocities, a numerical model is established according to the specific penetration experiment conditions. By comparing the predicted values of numerical simulation and the measured values of the experiment, the reliability and accuracy of the wide pressure range equation of state of limestone are verified. As Figure 6 shown, the numerical simulation results based on the established wide pressure range equation of state of limestone are in good agreement with the measured values, while the numerical simulation results based on the original equation of state without removing the shear strength effect and strain rate effect deviate greatly from the measured values, thus verifying that the wide pressure range equation of state based on limestone can achieve high - precision numerical simulation of the dynamic damage of concrete - like materials under explosive shock loads.

[0041] In order to better implement the equation of state determination method in the embodiments of the present invention, correspondingly, on the basis of the equation of state determination method, such as Figure 7 shown, the embodiments of the present invention also provide an equation of state determination device. The equation of state determination device 700 includes: The initial simulated axial stress determination module 701 is configured to use the real experimental data as the input of the numerical iteration model and output the initial simulated axial stress. The real experimental data is the measured axial stress and the measured volume strain obtained by the experimental object in the real plate impact experiment. The transitional simulated axial stress determination module 702 is configured to adjust the measured axial stress based on the initial simulated axial stress to obtain the adjusted measured axial stress, and use the measured volume strain and the adjusted measured axial stress as the input of the numerical iteration model to output the transitional simulated axial stress. The pressure equation of state construction module 703 is configured to, when the measured axial stress does not need to be adjusted, determine the current transitional simulated axial stress as the final axial stress, and use the final axial stress as the hydrostatic pressure to construct the wide pressure range equation of state of the experimental object. The equation of state determination device 700 provided in the above embodiments can implement the technical solutions described in the above equation of state determination method embodiments. The specific implementation principles of the above modules or units can be referred to the corresponding content in the above equation of state determination method embodiments, and will not be elaborated here.

[0042] As Figure 8 shown, the present invention also correspondingly provides a simulation device 800. The simulation device 800 includes a processor 801, a memory 802, and a display 803. Figure 8 Only some components of the simulation device 800 are shown, but it should be understood that it is not required to implement all the shown components, and more or fewer components can be alternatively implemented.

[0043] In some embodiments, the processor 801 may be a central processing unit (CPU), a microprocessor, or other data processing chips, and is configured to run the program code stored in the memory 802 or process data, such as the equation of state determination method in the present invention.

[0044] In some embodiments, the processor 801 may be a single server or a server group. The server group may be centralized or distributed. In some embodiments, the processor 801 may be local or remote. In some embodiments, the processor 801 may be implemented on a cloud platform. In some embodiments, the cloud platform may include a private cloud, a public cloud, a hybrid cloud, a community cloud, a distributed cloud, an internal cloud, a multi-cloud, etc., or any combination of the above.

[0045] The memory 802 can be an internal storage unit of the analog device 800 in some embodiments, such as the hard disk or memory of the analog device 800. The memory 802 can also be an external storage device of the analog device 800 in other embodiments, such as a plug-in hard disk equipped on the analog device 800, a Smart Media Card (SMC), a Secure Digital (SD) card, a Flash Card, etc.

[0046] Furthermore, the memory 802 can also include both the internal storage unit and the external storage device of the analog device 800. The memory 802 is used to store the application software and various types of data for installing the analog device 800.

[0047] The display 803 can be an LED display, a liquid crystal display, a touch liquid crystal display, and an OLED (Organic Light-Emitting Diode) toucher, etc. in some embodiments. The display 803 is used to display the information of the analog device 800 and to display a visual user interface. The components 801 - 803 of the analog device 800 communicate with each other through a system bus.

[0048] In some embodiments, when the processor 801 executes the state equation determination program in the memory 802, the following steps can be implemented: Using the real experimental data as the input of the numerical iteration model, an initial simulated axial stress is output, where the real experimental data is the measured axial stress and the measured volume strain of the experimental object in a real plate impact experiment; Adjusting the measured axial stress based on the initial simulated axial stress to obtain an adjusted measured axial stress, and using the measured volume strain and the adjusted measured axial stress as the input of the numerical iteration model, and outputting a transitional simulated axial stress; When the measured axial stress does not need to be adjusted, determining the current transitional simulated axial stress as the final axial stress, and using the final axial stress as the hydrostatic pressure to construct the state equation of the wide pressure range of the experimental object.

[0049] It should be understood that when the processor 801 executes the state equation determination program in the memory 802, in addition to the above functions, other functions can also be implemented. For specific details, please refer to the description of the corresponding method embodiments above.

[0050] Furthermore, embodiments of the present invention do not specifically limit the type of the mentioned analog device 800. The analog device 800 may be a portable analog device such as a mobile phone, a tablet computer, a personal digital assistant (PDA), a wearable device, a laptop, etc. Exemplary embodiments of the portable analog device include, but are not limited to, portable analog devices running IOS, android, microsoft, or other operating systems. The above-mentioned portable analog device may also be other portable analog devices, such as a laptop with a touch-sensitive surface (e.g., a touch panel). It should also be understood that in some other embodiments of the present invention, the analog device 800 may not be a portable analog device, but a desktop computer with a touch-sensitive surface (e.g., a touch panel).

[0051] Correspondingly, embodiments of the present application also provide a computer-readable storage medium. The computer-readable storage medium is used to store computer-readable programs or instructions. When the programs or instructions are executed by a processor, the steps or functions in the state equation determination method provided by the above-mentioned method embodiments can be implemented.

[0052] Those skilled in the art can understand that all or part of the processes of implementing the methods in the above embodiments can be completed by instructing relevant hardware through a computer program. The program can be stored in a computer-readable storage medium. Among them, the computer-readable storage medium is a disk, an optical disk, a read-only memory, or a random access memory, etc.

[0053] The above is only a preferred specific embodiment of the present invention, but the protection scope of the present invention is not limited thereto. Any changes or substitutions that can be easily thought of by those skilled in the art within the technical scope disclosed by the present invention should be covered by the protection scope of the present invention.

Claims

1. A method for determining a pressure state equation, characterized in that: include: Using real experimental data as input of the numerical iteration model, and outputting an initial simulated axial stress, wherein the real experimental data is the measured axial stress and the measured volume strain of the experimental object in a real flat plate impact experiment; The measured axial stress is adjusted based on the initial simulated axial stress to obtain the adjusted measured axial stress, the measured volume strain and the adjusted measured axial stress are used as inputs of the numerical iteration model, and the transition simulated axial stress is obtained as output; When the measured axial stress does not need to be adjusted, the current transitional simulated axial stress is determined as the final axial stress, and the final axial stress is used as the hydrostatic pressure to construct a wide pressure range state equation of the experimental object.

2. The method for determining the pressure state equation according to claim 1, characterized in that: The actual experimental data include a plurality of measured axial stresses and measured volumetric strains corresponding to the measured axial stresses, and the measured axial stresses cover a wide pressure range that the experimental object is subjected to under explosion impact loads.

3. The method for determining the pressure state equation according to claim 2, characterized in that: The actual experimental data is used as the input of the numerical iteration model, and the initial simulated axial stress is obtained as the output, including: The actual test data is used as the input of the state equation for representing the axial stress-strain volume relationship of the experimental object in the numerical iteration model, and the impact velocity of the flying piece hitting the test object in the numerical iteration model is adjusted to obtain the simulated axial stress corresponding to different measured axial stresses.

4. The method for determining the pressure state equation according to claim 3, characterized in that: The adjustment process for multiple measured axial stresses includes: The measured axial stress is adjusted in sequence according to the magnitude order of the impact velocity corresponding to the measured axial stress.

5. The method for determining the pressure state equation according to claim 1, characterized in that: The step of adjusting the measured axial stress based on the initial simulated axial stress to obtain the adjusted measured axial stress includes: Calculating a first difference between the measured axial stress and the initial simulated axial stress, and calculating an absolute value of a ratio of the first difference to the measured axial stress; When the absolute value of the ratio is greater than or equal to a preset threshold, a second difference between the measured axial stress and the first difference is used as the adjusted measured axial stress.

6. The method for determining the pressure state equation according to claim 5, characterized in that: The condition that the measured axial stress does not need to be adjusted is that the absolute value of the ratio is less than the preset threshold.

7. The method for determining the pressure state equation according to item 1 is characterized in that: The method of constructing a wide pressure range state equation of the experimental object by using the simulated axial stress that meets the preset conditions as the hydrostatic pressure includes: Using the simulated axial stress that meets the preset conditions as the hydrostatic pressure to construct a hydrostatic pressure-volume strain relationship of the experimental object; The state equation of the experimental object over a wide pressure range is constructed based on the hydrostatic pressure-volume strain relationship.

8. A device for determining a pressure state equation, characterized in that: include: An initial simulated axial stress determination module is used to use real experimental data as input of a numerical iteration model and output an initial simulated axial stress, wherein the real experimental data is the measured axial stress and measured volume strain obtained by the experimental object in a real flat plate impact experiment; A transition simulation axial stress determination module is used to adjust the measured axial stress based on the initial simulation axial stress to obtain the adjusted measured axial stress, and use the measured volume strain and the adjusted measured axial stress as inputs of the numerical iteration model to obtain the transition simulation axial stress as output; The pressure state equation construction module is used to determine the current transition simulation axial stress as the final axial stress when the measured axial stress does not need to be adjusted, and to construct the wide pressure range state equation of the experimental object using the final axial stress as the hydrostatic pressure.

9. A simulation device, characterized in that: comprising a memory and a processor, wherein: The memory is used to store programs; The processor is coupled to the memory and is used to execute the program stored in the memory to implement the steps in the method for determining the pressure state equation as described in any one of claims 1 to 7 above.

10. A computer-readable storage medium, characterized in that: Used to store computer-readable programs or instructions, which, when executed by a processor, can implement the steps in the method for determining the pressure state equation as described in any one of claims 1 to 7.

Images