High-precision digital penetration simulation method, device and equipment and storage medium

Abstract

The invention discloses a high-precision digital penetration simulation method. The method comprises the following steps: constructing a penetration body finite element model and a target finite element model; obtaining penetration conditions of the penetration body, wherein the penetration conditions comprise an attack angle, a penetration rate and a penetration angle; obtaining a contact relation and a boundary condition of the penetration body and the target; on the basis of the conditions, simulation of the penetration process is carried out; comparing a simulation result of the penetration process with an actual test result; and if the deviation between the simulation result of the penetration process and the actual test result is within the threshold value, outputting a penetration body overload curve, and if the deviation is not within the threshold value, correcting the penetration body finite element model and the target finite element model, and repeating the steps. The invention further discloses high-precision digital penetration simulation equipment and a storage medium. According to the high-precision digital penetration simulation method provided by the invention, a digital penetration simulation result can be ensured to be consistent with a real object test result, and the precision of the digital penetration simulation method is improved.

Description

Technical Field

[0001] The present application relates to the field of digital simulation technology, and more specifically, to a high-precision digital penetration simulation method, device and storage medium. Background Art

[0002] Digital penetration simulation technology is widely used in the product development and research of equipment systems due to its many advantages such as flexibility, repeatability and economy.

[0003] However, due to the complexity of the penetration control system and the target, the similarity between the digital simulation model and the physical model determines the accuracy and precision of the penetration simulation results. Low-precision or even erroneous simulation results not only fail to play a guiding role in engineering activities such as product development and production testing, but may also provide wrong guidance, leading to the emergence of "problematic" products and "problematic" tests, seriously affecting the normal performance of the penetration control system and even the entire equipment system.

[0004] The existing technology is unable to locate the key factors that affect the digital penetration simulation results, and it is difficult to ensure that the digital penetration simulation results are consistent with the actual test results, which reduces the accuracy of the digital penetration simulation method. Summary of the invention

[0005] In response to at least one defect or improvement need of the prior art, the present invention provides a high-precision digital penetration simulation method, device and storage medium, which can at least solve one of the above-mentioned problems in the prior art.

[0006] To achieve the above object, according to a first aspect of the present invention, a high-precision digital penetration simulation method is provided, the method comprising:

[0007] S1. Construct the finite element model of the penetrator and the finite element model of the target;

[0008] S2. Obtaining the penetration conditions of the penetrator, wherein the penetration conditions include the angle of attack, the penetration speed and the penetration angle; obtaining the contact relationship and boundary conditions between the penetrator and the target;

[0009] S3. Based on the penetrator finite element model, the target finite element model, the penetration condition, and the contact relationship and boundary conditions between the penetrator and the target, perform a penetration process simulation;

[0010] S4. Compare the results of the penetration process simulation with the actual test results. If the deviation between the results of the penetration process simulation and the actual test results is within a preset threshold, output the penetrator overload curve; if the deviation between the results of the penetration process simulation and the actual test results is not within the threshold, correct the penetrator finite element model and the target finite element model, and repeat steps S2 to S4.

[0011] Furthermore, in the above-mentioned high-precision digital penetration simulation method, the finite element model of the penetrator adopts a method of mapping the penetrator head block, and the penetrator is divided into hexahedral units.

[0012] Furthermore, in the above-mentioned high-precision digital penetration simulation method, the threaded connection between the penetrating body components is simulated by using a bound contact or common node method.

[0013] Furthermore, in the above-mentioned high-precision digital penetration simulation method, the construction of the finite element model of the penetrator is based on a constitutive relationship model, and the constitutive relationship model is

[0014]

[0015] Among them, σ y is the yield stress, A, B, n, C and m are the constitutive model parameters of metal materials, obtained through experiments, ε p is the effective plastic strain, is the reference strain rate, is the actual strain rate, T is the actual temperature, T m is the melting point of the material and T0 is the reference temperature.

[0016] Furthermore, in the above-mentioned high-precision digital penetration simulation method, in the target finite element model, the ratio of the penetrator unit size to the target unit size is between 1 and 1.2.

[0017] Furthermore, in the above-mentioned high-precision digital penetration simulation method, when the target is a non-homogeneous target with uneven internal structure, if the target unit size is smaller than the penetrator head unit size, a point-to-surface erosion contact keyword between the target and the penetrator is set in the simulation calculation file, so that the target is formed as a slave surface and the penetrator is formed as a master surface in the target finite element model.

[0018] Furthermore, in the above-mentioned high-precision digital penetration simulation method, when the target is a non-homogeneous target with an uneven internal structure, if the target unit size is larger than the penetrator head unit size, the number of grids in the target finite element model that are in contact with the penetrator is 2 to 5 times the number of grids that are not in contact with the penetrator.

[0019] Furthermore, in the above-mentioned high-precision digital penetration simulation method, the construction of the finite element model of the penetrator is based on the material damage model of the penetrator, and the material damage model of the penetrator includes a damage evolution equation and a yield surface equation, wherein the damage evolution equation is

[0020]

[0021] Among them, Δε pis the equivalent plastic strain within a calculation cycle, Δμ p is the plastic volume strain within a calculation cycle, is the equivalent plastic strain of crushing under normal pressure, is the plastic volume strain of crushing under normal pressure;

[0022] The yield surface equation is

[0023]

[0024] Among them, σ * is the actual equivalent stress, P * is the hydrostatic pressure; is the strain rate; D is the damage degree; A' and B' are the strength parameters of the material.

[0025] According to a second aspect of the present invention, there is also provided a high-precision digital penetration simulation device, which includes at least one processing unit and at least one storage unit, wherein the storage unit stores a computer program, and when the computer program is executed by the processing unit, the processing unit executes the steps of any one of the above-mentioned methods.

[0026] According to the third aspect of the present invention, there is also provided a storage medium storing a computer program executable by a high-precision digital penetration simulation device, wherein when the computer program runs on the high-precision digital penetration simulation device, the high-precision digital penetration simulation device executes the steps of any one of the above-mentioned methods.

[0027] In general, the above technical solutions conceived by the present invention can achieve the following beneficial effects compared with the prior art:

[0028] The present invention provides a high-precision digital penetration simulation method, which optimizes the digital simulation model of the penetrator from the aspects of penetrator finite element model construction, material constitutive selection, material dynamic mechanical parameter setting, etc., and can obtain high-precision digital penetration simulation results by solving the digital simulation model, which can ensure that the digital penetration simulation results are consistent with the actual test results, thereby improving the accuracy of the digital penetration simulation method. BRIEF DESCRIPTION OF THE DRAWINGS

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

[0030] Figure 1A flowchart of a high-precision digital penetration simulation method provided in an embodiment of the present application;

[0031] Figure 2 A schematic diagram of a digital penetration simulation calculation process provided in an embodiment of the present application;

[0032] Figure 3 Schematic diagram of a finite element network model of a penetration control system formed using a common node method provided in an embodiment of the present application

[0033] Figure 4 A schematic diagram of a finite element mesh of the unit sizes of the penetrator and target provided in an embodiment of the present application;

[0034] Figure 5 A schematic diagram of a finite element mesh showing the sparse and dense distribution of target units provided in an embodiment of the present application. DETAILED DESCRIPTION

[0035] In order to make the purpose, technical solutions and advantages of the present invention more clearly understood, the present invention is further described in detail below in conjunction with the accompanying drawings and embodiments. It should be understood that the specific embodiments described herein are only used to explain the present invention and are not intended to limit the present invention. In addition, the technical features involved in the various embodiments of the present invention described below can be combined with each other as long as they do not conflict with each other.

[0036] The terms "first", "second", "third", etc. in the specification and claims of this application and the above-mentioned drawings are used to distinguish different objects, rather than to describe a specific order. In addition, the terms "including" and "having" and any variations thereof are intended to cover non-exclusive inclusions. For example, a process, method, system, product or device that includes a series of steps or units is not limited to the listed steps or units, but optionally includes steps or units that are not listed, or optionally includes other steps or units inherent to these processes, methods, products or devices.

[0037] Figure 1 is a flow chart of a high-precision digital penetration simulation method provided by an embodiment of the present application, such as Figure 1 As shown, including:

[0038] S1. Construct the finite element model of the penetrator and the finite element model of the target.

[0039] Specifically, digital penetration simulation belongs to the category of transient impact dynamics simulation. The finite element explicit solution algorithm can be used to calculate the large deformation dynamic response of nonlinear structural materials. The Lagrangian description increment method and the central difference method are preferably used for time integration in the algorithm.

[0040] In each time increment step, a series of calculations are performed on each Lagrangian subgrid. The node acceleration is calculated through boundary conditions, internal forces, and external forces. The node velocity is obtained after integration, and the node displacement is obtained after integration again. The volume and strain rate are calculated by the mass conservation equation combined with the initial conditions. The pressure, internal energy, and deviatoric stress are calculated through the material model and energy equation. The node force is calculated through the momentum conservation equation to complete a calculation cycle. Figure 2 shown.

[0041] The existing digital penetration simulation results are obtained by solving the digital penetration simulation model, and the digital penetration simulation model contains the simplification degree of the penetration control system geometric model, the finite element mesh model of the penetration control system and the target, material properties, loads and boundary conditions, and the interaction between the components inside the penetration control system and the target. Among them, for the simplification of the penetration control system geometric model, since the structural composition of different penetration control systems is different, the simplification degree of the geometric model needs to be determined according to the specific problem, which is not a common factor affecting the digital penetration simulation results, and therefore is not a factor affecting the accuracy of digital penetration simulation.

[0042] The key influencing factors in the digital penetration simulation model are the finite element mesh model of the penetration control system and the target, the material properties, and the interaction between the components inside the penetration control system and between the penetration control system and the target. These key influencing factors and other information reflecting the characteristics of the physical model constitute a complete digital penetration simulation model, which is a program that can be recognized by the computer, and the digital penetration simulation results are obtained after solving.

[0043] Therefore, based on the material state equation, constitutive relationship model, material damage model, boundary constraints and contact connection relationship of the penetrator, the finite element model of the penetrator is obtained; based on the material composition of the target, the finite element model of the target is obtained. Among them, the method of constructing the finite element model is a prior art means in this field, and this application will not be repeated here.

[0044] S2. Obtain the penetration conditions of the penetrator, which include the angle of attack, penetration speed and penetration angle; obtain the contact relationship and boundary conditions between the penetrator and the target.

[0045] Specifically, for loads and boundary conditions, the penetration speed, impact angle, attack angle and target fixing method of the penetrator in the digital penetration simulation are easy to keep consistent with the physical test. The penetration conditions of the penetrator are obtained, including the attack angle, penetration speed and penetration angle; the contact relationship and boundary conditions between the penetrator and the target are obtained.

[0046] S3. Based on the finite element model of the penetrator, the finite element model of the target, the penetration conditions, the contact relationship between the penetrator and the target, and the boundary conditions, the penetration process is simulated.

[0047] Specifically, the penetration process simulation can be performed using any simulation software such as PANDA-Impact, LS-DYNA or Explicit Dynamics. The penetration process simulation can be performed and the simulation results can be obtained by inputting the finite element model of the penetrator, the finite element model of the target, the penetration conditions, the contact relationship between the penetrator and the target, and the boundary conditions into the simulation software.

[0048] S4. Compare the results of the penetration process simulation with the actual test results. If the deviation between the results of the penetration process simulation and the actual test results is within the threshold, output the penetrator overload curve; if the deviation between the results of the penetration process simulation and the actual test results is not within the threshold, correct the penetrator finite element model and the target finite element model, and repeat steps S2 to S4.

[0049] Specifically, the threshold value represents the threshold value at which the deviation between the penetration simulation result and the actual test result is within an acceptable range, and can be set according to engineering needs, and this application does not make specific requirements.

[0050] If the deviation between the results of the penetration process simulation and the actual test results is within the threshold, it means that the deviation between the penetration simulation results and the actual test results is within an acceptable range. At this time, the settings of the penetration simulation can be used for the actual simulation process. At this time, the penetrator overload curve is output. The penetrator overload curve is the penetrator acceleration change curve. This penetrator overload curve can be used in subsequent simulations.

[0051] If the deviation between the results of the penetration process simulation and the actual test results is not within the threshold, the penetrator finite element model and the target finite element model are corrected according to the penetration process simulation results, and the corrected penetrator finite element model and the target finite element model are re-input into the penetration process simulation software, and steps S1 to S4 are repeated to judge the penetration process simulation results again until the deviation between the results of the penetration process simulation and the actual test results is within the threshold.

[0052] A high-precision digital penetration simulation method provided in an embodiment of the present application optimizes the digital simulation model of the penetration control system from many aspects, such as the construction of a finite element model of the penetration control system, the selection of material constitutive properties, and the setting of material dynamic mechanical parameters. By solving the digital simulation model, high-precision digital penetration simulation results can be obtained, which can ensure that the digital penetration simulation results are consistent with the actual test results, thereby improving the accuracy of the digital penetration simulation method.

[0053] Optionally, in the high-precision digital penetration simulation method provided in the embodiment of the present application, the finite element model of the penetrator adopts a method of mapping the head part of the penetrator to divide the entire penetrator into hexahedral units.

[0054] Specifically, the head of the penetration control system is generally truncated oval, with a smaller blunt or rounded head. If the penetrator is swept along the circumference, the penetrator head will inevitably generate some pentahedral units, and the size of the generated pentahedral units along the circumference will be very small. If the penetrator head unit contacts the much larger unit in the target, unit penetration or unit distortion will occur during the solution, causing the program to terminate with an error. For the structural features of the penetrator head, the method of mapping the penetrator head partial blocks can be used to ensure that the penetrator is completely divided into hexahedral units and the size of the penetrator head unit can be controlled within a reasonable range.

[0055] Optionally, in the high-precision digital penetration simulation method provided in the embodiment of the present application, the threaded connection between the penetration body components is simulated using a bound contact or common node method.

[0056] Specifically, in general, the penetration body and the control body base, and the control body base and the control body are connected by threads. It is very difficult to establish the geometric entity of the threaded pair, which not only affects the layout of the unit size of each component in the overall simulation model, but also increases the number of units in the overall simulation model.

[0057] Figure 3 Figure 1 is a schematic diagram of the finite element network model of the penetration control system formed using the common node method. Figure 3 As shown, for the threaded connection between components, bound contact or common nodes can be used for processing. Binding contact does not require node coordination of the threaded connection part, but because its contact search is unidirectional, the master and slave surfaces in the contact must be strictly distinguished; common nodes do not need to be set up for contact in the K file, but the node coordination of the threaded connection part must be guaranteed in the pre-processing. Among them, the K file is a commonly used file in the penetration simulation software, which can be used to set the point-surface erosion contact keywords between the target and the penetrator, and this application will not go into details here. The position where the interaction occurs between non-threaded connectors in the penetration control system is defined by automatic contact, and the force load is transmitted through the nodes on each component during the penetration process.

[0058] Optionally, the high-precision digital penetration simulation method provided in the embodiment of the present application, the constitutive relationship model is a metal material constitutive model based on temperature and strain rate effects, and its expression is as follows:

[0059]

[0060] Among them, σy is the yield stress, A, B, n, C and m are the constitutive model parameters of metal materials, obtained through experiments, ε p is the effective plastic strain, is the reference strain rate, is the actual strain rate, T is the actual temperature, T m is the melting point of the material and T0 is the reference temperature.

[0061] Specifically, the material constitutive model refers to a mathematical model that reflects the specific macroscopic properties of the material. The material constitutive model and parameter settings are key factors affecting the accuracy of digital simulation. When the penetrator and the target meet, the material is in a high temperature and high strain rate state, and the influence of temperature and strain rate on the mechanical properties of the material is very significant. Therefore, in order to improve the calculation accuracy, it is not appropriate to use a material constitutive model without temperature and strain rate settings when performing digital penetration simulation. It is necessary to select a material constitutive model that takes into account the influence of temperature and strain rate.

[0062] Johnson-Cook is a constitutive model for metal materials that takes into account the effects of temperature and strain rate.

[0063]

[0064] Optionally, in the high-precision digital penetration simulation method and target finite element model provided in an embodiment of the present application, when the target is a homogeneous target with a uniform internal structure, the ratio of the penetrator unit size to the target unit size is between 1 and 1.2, for example, 1, 1.1, 1.15 or other values ​​within this range.

[0065] Specifically, in order to ensure the stability of the solution, the target finite element mesh division needs to comply with the mesh unit size where the target may contact the penetrator and the mesh unit size of the penetrator head. In the case of a homogeneous target with a uniform internal structure, the ratio of the penetrator unit size to the target unit size is between 1 and 1.2.

[0066] When the target is a homogeneous target with a uniform internal structure, that is, in the case of a homogeneous target, for example, the target material is a steel plate, the ratio of the penetrator unit size to the target unit size is between 1 and 1.2. When setting the contact keywords of the subsequent K file, the main and secondary surface positions of the penetrator and the target can be correctly defined.

[0067] Optionally, in the high-precision digital penetration simulation method and target finite element model provided in the embodiment of the present application, when the target is a non-homogeneous target with uneven internal structure, if the target unit size is smaller than the penetrator head unit size, the point-to-surface erosion contact keywords between the target and the penetrator are set in the K file to form the target as the slave surface and the penetrator as the master surface to meet the contact definition criterion requirements.

[0068] Specifically, in the case where the target is a non-homogeneous target with an uneven internal structure, in one embodiment, the target is made of reinforced concrete. According to the connection method between the steel beam unit and the concrete entity unit, if the target unit size is smaller than the penetrator head unit size, the steel beam unit size is the same as the concrete unit size, that is, the steel beam unit size is also smaller than the penetrator head unit size, then the connection between the steel bar and the concrete adopts a common node method, and it is no longer necessary to define the keywords for the contact between the steel bar and the concrete in the K file. In the K file, the point-to-surface erosion contact keywords between the steel bar and the penetrator are set to form a steel bar node group as a slave surface to meet the contact definition criteria. In the surface-to-surface erosion contact keywords between the concrete and the penetrator are set, the target is used as a slave surface and the penetrator is used as a master surface to meet the contact definition criteria.

[0069] Optionally, in the high-precision digital penetration simulation method and target finite element model provided in the embodiment of the present application, when the target is a non-homogeneous target with an uneven internal structure, if the target unit size is larger than the penetrator head unit size, the number of grids in the contact part of the penetrator is 2 to 5 times the number of grids in the non-contact part of the penetrator, for example, 2 times, 3 times, 4 times, 5 times or other multiples within this range.

[0070] Specifically, in one embodiment, the target is made of reinforced concrete. If the target unit size is larger than the unit size of the penetrator head, and the steel beam unit size is also smaller than the unit size of the penetrator head, the steel bar and concrete cannot be connected in a common node manner, because the common node connection must ensure that the steel beam unit size is equal to the concrete unit size. At this time, the steel beam unit size is larger than the penetrator head unit size, and in the point-surface erosion contact keyword between the steel bar and the penetrator, the secondary surface must be the node group that constitutes the steel bar. At this time, the secondary surface unit is larger than the main surface unit, which violates the contact definition criterion requirements, and contact failure will occur in the solution, and the penetrator and the steel bar will penetrate each other. At this time, under the condition of "steel beam unit size < penetrator head unit size < concrete unit size", the point-surface binding contact is defined for the steel bar and the concrete, and the surface-surface erosion contact between the concrete and the penetrator and the point-surface penetration contact between the steel bar and the penetrator play a role at the same time, establishing the interaction relationship between it and the penetrator, concrete, and steel bar.

[0071] While controlling the size of the target unit in contact with the penetrator, the density of the unit is considered. The mesh in contact with the penetrator is finer, and the mesh in non-contact with the penetrator is sparser, that is, the number of meshes in contact with the penetrator is 2 to 5 times that of meshes in non-contact with the penetrator. The schematic diagram of the finite element mesh reflecting the size of the penetrator and target unit and the schematic diagram of the finite element mesh of the target unit density distribution are shown in Figure 1. Figure 4 , Figure 5 shown.

[0072] Optionally, in the high-precision digital penetration simulation method provided in the embodiment of the present application, the material damage model of the penetrator includes a damage evolution equation and a yield surface equation, wherein the damage evolution equation is

[0073]

[0074] Among them, Δε p is the equivalent plastic strain within a calculation cycle, Δμ p is the plastic volume strain within a calculation cycle, is the equivalent plastic strain of crushing under normal pressure, is the plastic volume strain of crushing under normal pressure;

[0075] The yield surface equation is

[0076]

[0077] Among them, σ * is the actual equivalent stress, P * is the hydrostatic pressure; is the strain rate; D is the damage degree; A' and B' are the strength parameters of the material.

[0078] Specifically, the concrete material can preferably adopt the JOHNSON_HOLMQUIST_CONCRETE constitutive model, which has a yield surface equation and a damage evolution equation.

[0079] In one embodiment, the yield surface equation is

[0080]

[0081] Among them, σ * is the actual equivalent stress, P * is the hydrostatic pressure; is the strain rate; D is the damage degree; A' and B' are the strength parameters of the material.

[0082] The damage evolution equation is:

[0083]

[0084] In one embodiment, the target is made of concrete, and its constitutive mechanical parameters can be set as shown in the following Table 1. However, the constitutive mechanical parameters of the present application are not limited thereto, and other parameters that can meet the simulation requirements are possible.

[0085] Table 1 Constitutive mechanical parameters of concrete materials

[0086]

[0087]

[0088] The present application also provides a computer-readable storage medium on which a computer program is stored, and when the program is executed by a processor, the steps of the above method are implemented. The computer-readable storage medium may include, but is not limited to, any type of disk, including a floppy disk, an optical disk, a DVD, a CD-ROM, a micro drive, and a magneto-optical disk, a ROM, a RAM, an EPROM, an EEPROM, a DRAM, a VRAM, a flash memory device, a magnetic card or an optical card, a nanosystem (including a molecular memory IC), or any type of medium or device suitable for storing instructions and / or data.

[0089] It should be noted that, for the aforementioned method embodiments, for the sake of simplicity, they are all expressed as a series of action combinations, but those skilled in the art should be aware that the present application is not limited by the described order of actions, because according to the present application, certain steps can be performed in other orders or simultaneously. Secondly, those skilled in the art should also be aware that the embodiments described in the specification are all preferred embodiments, and the actions and modules involved are not necessarily required by the present application.

[0090] In the above embodiments, the description of each embodiment has its own emphasis. For parts that are not described in detail in a certain embodiment, reference can be made to the relevant descriptions of other embodiments.

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

[0092] The units described as separate components may or may not be physically separated, and the components shown as units may or may not be physical units, that is, they may be located in one place or distributed on multiple network units. Some or all of the units may be selected according to actual needs to achieve the purpose of the solution of this embodiment.

[0093] In addition, each functional unit in each embodiment of the present application may be integrated into one processing unit, or each unit may exist physically separately, or two or more units may be integrated into one unit. The above-mentioned integrated unit may be implemented in the form of hardware or in the form of software functional units.

[0094] If the integrated unit is implemented in the form of a software functional unit and sold or used as an independent product, it can be stored in a computer-readable memory. Based on this understanding, the technical solution of the present application is essentially or the part that contributes to the prior art or all or part of the technical solution can be embodied in the form of a software product, and the computer software product is stored in a memory, including a number of instructions to enable a computer device (which can be a personal computer, server or network device, etc.) to perform all or part of the steps of the method described in each embodiment of the present application. The aforementioned memory includes: U disk, read-only memory (ROM), random access memory (RAM), mobile hard disk, disk or optical disk and other media that can store program code.

[0095] A person of ordinary skill in the art can understand that all or part of the steps in the various methods of the above embodiments can be completed by instructing related hardware through a program, and the program can be stored in a computer-readable memory, which can include: a flash drive, a read-only memory (ROM), a random access memory (RAM), a magnetic disk or an optical disk, etc.

[0096] The above is only an exemplary embodiment of the present disclosure, and the scope of the present disclosure cannot be limited thereto. That is, any equivalent changes and modifications made according to the teachings of the present disclosure are still within the scope of the present disclosure. After considering the specification and practicing the disclosure here, those skilled in the art will easily think of the implementation scheme of the present disclosure. This application is intended to cover any modification, use or adaptation of the present disclosure, which follows the general principles of the present disclosure and includes common knowledge or customary technical means in the technical field not recorded in the present disclosure. The description and examples are regarded as exemplary only, and the scope and spirit of the present disclosure are defined by the claims.

[0097] The technical features of the above embodiments may be combined arbitrarily. To make the description concise, not all possible combinations of the technical features in the above embodiments are described. However, as long as there is no contradiction in the combination of these technical features, they should be considered to be within the scope of this specification.

[0098] It will be easily understood by those skilled in the art that the above description is only a preferred embodiment of the present invention and is not intended to limit the present invention. Any modifications, equivalent substitutions and improvements made within the spirit and principles of the present invention should be included in the protection scope of the present invention.

Claims

1. A high-precision digital penetration simulation method, characterized in that: include: S1. Construct the finite element model of the penetrator and the finite element model of the target; S2. Obtaining the penetration conditions of the penetrator, wherein the penetration conditions include the angle of attack, the penetration speed and the penetration angle; obtaining the contact relationship and boundary conditions between the penetrator and the target; S3. Based on the penetrator finite element model, the target finite element model, the penetration condition, and the contact relationship and boundary conditions between the penetrator and the target, perform a penetration process simulation; S4. Compare the result of the penetration process simulation with the actual test result. If the deviation between the result of the penetration process simulation and the actual test result is within a preset threshold, output the overload curve of the penetrator; If the deviation between the result of the penetration process simulation and the actual test result is not within the threshold, the penetrator finite element model and the target finite element model are corrected, and steps S2 to S4 are repeated.

2. The high-precision digital penetration simulation method according to claim 1, characterized in that: The penetrator finite element model adopts a method of mapping the penetrator head block, and the penetrator is divided into hexahedral units.

3. The high-precision digital penetration simulation method according to claim 1 or 2, characterized in that: The threaded connection between the penetrating body components is simulated by using a binding contact or a common node method.

4. The high-precision digital penetration simulation method according to claim 1, characterized in that: The construction of the finite element model of the penetrator is based on a constitutive relationship model, which is Among them, σ y is the yield stress, A, B, n, C and m are the constitutive model parameters of metal materials, obtained through experiments, ε p is the effective plastic strain, is the reference strain rate, is the actual strain rate, T is the actual temperature, T m is the melting point of the material and T0 is the reference temperature.

5. The high-precision digital penetration simulation method according to claim 1, characterized in that: In the target finite element model, the ratio of the penetrator unit size to the target unit size is between 1 and 1.

2.

6. The high-precision digital penetration simulation method according to claim 1, characterized in that: In the case where the target is a non-homogeneous target with an uneven internal structure, if the target unit size is smaller than the penetrator head unit size, a point-surface erosion contact keyword between the target and the penetrator is set in the simulation calculation file so that the target is formed as a slave surface and the penetrator is formed as a master surface in the target finite element model.

7. The high-precision digital penetration simulation method according to claim 1, characterized in that: In the case where the target is a non-homogeneous target with an uneven internal structure, if the target unit size is larger than the penetrator head unit size, the number of meshes in the penetrator contact portion in the target finite element model is 2 to 5 times the number of meshes in the non-penetrator contact portion.

8. The high-precision digital penetration simulation method according to claim 1, characterized in that: The construction of the penetrator finite element model is based on the material damage model of the penetrator, and the material damage model of the penetrator includes a damage evolution equation and a yield surface equation, wherein the damage evolution equation is Among them, Δε p is the equivalent plastic strain within a calculation cycle, Δμ p is the plastic volume strain within a calculation cycle, is the equivalent plastic strain of crushing under normal pressure, is the plastic volume strain of crushing under normal pressure; The yield surface equation is Among them, σ * is the actual equivalent stress, P * is the hydrostatic pressure; is the strain rate; D is the damage degree; A' and B' are the strength parameters of the material.

9. A high-precision digital penetration simulation device, characterized in that: The method comprises at least one processing unit and at least one storage unit, wherein the storage unit stores a computer program, and when the computer program is executed by the processing unit, the processing unit executes the steps of the method according to any one of claims 1 to 8.

10. A storage medium, characterized in that: It stores a computer program executable by a high-precision digital penetration simulation device. When the computer program runs on the high-precision digital penetration simulation device, the high-precision digital penetration simulation device executes the steps of the method described in any one of claims 1 to 8.

Images