Parallel coupling material point finite difference finite element simulation method and device

Abstract

A parallel coupling material point finite difference finite element simulation method is applied to material explosion simulation and comprises the steps that background grids are generated, each background grid comprises a plurality of material points, material information of the material points is initialized, and the material information comprises at least one of mass information, momentum information, density information, stress strain information and deformation gradient information; initializing the rod unit; based on a finite difference finite element method, performing parallel simulation and communication among the background grids; and moving the grid data, including movement of the substance points and movement of the rod units. According to the method, the simulation precision can be improved when the problem of small deformation is solved, so that the simulation better fits the actual situation.

Description

Technical Field

[0001] The present invention relates to the field of large-scale supercomputing technology, and in particular to a parallel coupled material point finite difference finite element simulation method and device. Background Art

[0002] The material point method is very suitable for dealing with extreme events, such as hypervelocity collisions, geotechnical problems, explosions, crack evolution, incompressible fluids and other applications. For fluid-solid coupling processes with strong nonlinearity, the coupled material point finite difference method divides the problem domain into a finite difference area and a material point area, uses finite differences to simulate fluid propagation, and uses the material point method to simulate solid response, which not only takes advantage of the finite difference method for solving large-scale flow field problems, but also takes advantage of the material point method for solving large structural deformation problems. However, when faced with small deformation problems, such as the simulation of the steel reinforcement part of a building in an explosion application, the accuracy and efficiency of the material point method are lower than the finite element method. Therefore, a solution to small deformation problems is urgently needed. Summary of the invention

[0003] In order to solve the problems existing in the prior art, the embodiments of the present application provide a method, apparatus, computing device, computer storage medium and product including a computer program for parallel coupled material point finite difference finite element simulation, which can improve the simulation accuracy and make the simulation more in line with the actual situation.

[0004] In the first aspect, an embodiment of the present application provides a parallel coupled material point finite difference finite element simulation method, which is applied to material explosion simulation, and the method includes: generating background grids, each background grid including a number of material points, initializing material information of the material points, the material information including at least one of mass, momentum, density, stress strain, and deformation gradient information; initializing rod units; based on the finite difference finite element method, parallel simulation and communication between each background grid; moving grid data, including movement of material points and movement of rod units.

[0005] In some possible implementations, generating a background grid specifically includes: generating a background grid using an octree structure to construct the background grid, determining the final shape of the background grid by encrypting the octree, and initial parameters and encryption rules of the octree need to be determined based on the spatial characteristics and accuracy of the material explosion simulation scene.

[0006] In some possible implementations, initializing the rod unit includes: setting structural parameters for the rod unit, the structural parameters including the length, cross-sectional shape and size, elastic modulus, and Poisson's ratio information of the rod unit, establishing a mapping relationship between the point structure of the rod unit and the material points, wherein the point structure participates in the material point method grid point mapping and calculation, and the rod structure is updated based on the position information of the two end points of the rod unit and the set structural parameters.

[0007] In some possible implementations, based on the finite difference and finite element method, parallel simulation and communication between various background grids include: dividing the background grids into three categories: fluid, solid, and fluid-solid coupling; using finite differences to calculate fluid dynamics related equations for fluid grids; using material point method or finite element method to calculate solid mechanics related equations for solid grids; using immersed boundary method to deal with interactions for fluid-solid coupling grids; establishing a neighbor list mechanism; when data transmission involves cross-processes, searching the subscript of the target grid in the adjacent process based on the neighbor list, and sending the data through data packaging and transmission protocol.

[0008] In some possible implementations, the mobile grid data includes: for the movement of material points, by establishing a material point position query function and a process attribution judgment function, the newly calculated position of the material point is compared with the background grid range that the current process is responsible for. When the material point exceeds the current process range, all the material information of the material point is packaged and moved to the corresponding adjacent process in a preset data format; for the movement of the rod unit, based on the position change of the two end points of the rod unit, it is judged whether it needs to be moved. If it needs to be moved, in addition to moving the rod structure and point structure of the rod unit to the corresponding process at the same time during the movement process; the index mapping relationship of the rod unit in the new process is updated.

[0009] In the second aspect, an embodiment of the present application provides a parallel coupled material point finite difference finite element simulation device, which is deployed on a computing device, and the computing device is applied to material explosion simulation. The device includes: an initialization module, which is used to generate background grids, each background grid includes a number of material points, and initializes the material information of the material points, and the material information includes at least one of mass, momentum, density, stress strain, and deformation gradient information; the initialization module is also used to initialize the rod unit; a processing module is used to perform parallel simulation and communication between each background grid based on the finite difference finite element method; the processing module is also used to move grid data, including the movement of material points and the movement of rod units.

[0010] In some possible implementations, the initialization module is specifically used to: generate a background grid constructed using an octree structure, determine the final shape of the background grid by encrypting the octree, and the initial parameters and encryption rules of the octree need to be determined based on the spatial characteristics and accuracy of the material explosion simulation scene.

[0011] In some possible implementations, the initialization module is also used to: set structural parameters for the rod unit, the structural parameters include the length, cross-sectional shape and size, elastic modulus, and Poisson's ratio information of the rod unit, establish a mapping relationship between the point structure of the rod unit and the material point, wherein the point structure participates in the material point method grid point mapping and calculation, and the rod structure is updated based on the position information of the two end points of the rod unit and the set structural parameters.

[0012] In some possible implementations, the processing module is used to: divide the background grid into three categories: fluid, solid, and fluid-solid coupling; use finite difference to calculate fluid dynamics related equations for the fluid grid; use material point method or finite element method to calculate solid mechanics related equations for the solid grid; use immersed boundary method to handle interactions for the fluid-solid coupling grid; establish a neighbor list mechanism; when data transmission involves cross-processes, search for the subscript of the target grid in the adjacent process based on the neighbor list, and send the data through data packaging and transmission protocol.

[0013] In some possible implementations, the processing module is also used for: for the movement of material points, by establishing a material point position query function and a process attribution judgment function, comparing the newly calculated position of the material point with the background grid range that the current process is responsible for, when the material point exceeds the current process range, packaging all the material information of the material point in a preset data format and moving it to the corresponding adjacent process; for the movement of the rod unit, judging whether it needs to be moved based on the position change of the two end points of the rod unit, if it needs to be moved, in addition to moving the rod structure and point structure of the rod unit to the corresponding process at the same time during the movement; updating the index mapping relationship of the rod unit in the new process.

[0014] In a third aspect, an embodiment of the present application provides a computer-readable storage medium, comprising computer-readable instructions. When a computer reads and executes the computer-readable instructions, the computer executes the method described in any one of the first aspects.

[0015] In a fourth aspect, an embodiment of the present application provides a computing device, comprising a processor and a memory, wherein the memory stores computer program instructions, and when the computer program instructions are executed by the processor, the method described in any one of the first aspects is executed.

[0016] In a fifth aspect, an embodiment of the present application provides a product comprising a computer program. When the computer program product runs on a processor, the processor executes the method as described in any one of the first aspects. BRIEF DESCRIPTION OF THE DRAWINGS

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

[0018] Figure 1 It is a flow chart of a parallel coupled material point finite difference finite element simulation method provided in an embodiment of the present application;

[0019] Figure 2is a schematic diagram of a rod unit division provided in an embodiment of the present application;

[0020] Figure 3 is a schematic diagram of a rod unit data transmission provided by an embodiment of the present application;

[0021] Figure 4 This is a schematic diagram of the movement of a rod unit between processes provided by an embodiment of the present application;

[0022] Figure 5 It is a structural schematic diagram of a parallel coupled material point finite difference finite element simulation device provided in an embodiment of the present application;

[0023] Figure 6 It is a structural diagram of a computing device provided in an embodiment of the present application. DETAILED DESCRIPTION

[0024] In order to make the purpose, technical solution and advantages of the embodiments of the present invention clearer, the technical solution in the embodiments of the present invention will be clearly and completely described below in conjunction with the drawings 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.

[0025] The term "and / or" in this article is a description of the association relationship of associated objects, indicating that there can be three relationships. For example, A and / or B can represent: A exists alone, A and B exist at the same time, and B exists alone. The symbol " / " in this article indicates that the associated objects are in an or relationship, for example, A / B means A or B.

[0026] The terms "first" and "second" in the specification and claims herein are used to distinguish different objects rather than to describe a specific order of the objects. For example, a first response message and a second response message are used to distinguish different response messages rather than to describe a specific order of the response messages.

[0027] In the embodiments of the present application, words such as "exemplary" or "for example" are used to indicate examples, illustrations or descriptions. Any embodiment or design described as "exemplary" or "for example" in the embodiments of the present application should not be interpreted as being more preferred or more advantageous than other embodiments or designs. Specifically, the use of words such as "exemplary" or "for example" is intended to present related concepts in a specific way.

[0028] In the description of the embodiments of the present application, unless otherwise specified, "multiple" means two or more than two. For example, multiple processing units refer to two or more processing units, etc.; multiple elements refer to two or more elements, etc.

[0029] To facilitate the understanding of the embodiments of the present application, the following will further explain and illustrate specific embodiments in conjunction with the accompanying drawings. The embodiments do not constitute a limitation on the embodiments of the present invention.

[0030] The material point method (MPM) discretizes the material region through particles, and each particle carries all the material information, such as mass, position, momentum, etc. During the calculation process, the material points move on the background space grid, and the background grid is based on the Euler format and will not change with the movement of the material points, thus avoiding the problem of grid distortion. By using a regular Euler background grid to calculate the spatial derivatives and momentum equations, the material point method realizes the interaction and connection between the particles. In addition, the material point method also avoids the processing of convection terms, which makes it advantageous in algorithm stability and efficiency. The material point method has become an effective method for solving strong nonlinear problems such as hypervelocity collisions, explosions, and impact penetration. In hypervelocity collisions, the material point method can accurately simulate the deformation, fragmentation, and sputtering of objects during hypervelocity collisions. In explosion simulations, the material point method can capture complex physical phenomena such as shock waves and debris scattering generated by explosions. In impact penetration, the material point method can simulate the dynamic response of materials under impact loads, such as the process of projectiles penetrating a target plate. In dynamic fracture simulation, the material point method can capture the processes of crack initiation, propagation and penetration. In fluid-solid coupling, the material point method can handle the interaction between fluid and solid, such as the process of droplets hitting the solid surface. In multi-scale analysis, the material point method can also be combined with other numerical methods to perform multi-scale analysis to capture the behavior of materials at different scales. However, when faced with some small deformation problems such as the steel reinforcement of buildings in explosion applications, the accuracy and efficiency of the material point method are not satisfactory.

[0031] In view of this, an embodiment of the present application provides a parallel coupled material point finite difference finite element simulation method, which uses finite differences to simulate fluid propagation, uses material point method to simulate solid response, and uses finite elements to simulate small deformation problems, making the simulation closer to reality.

[0032] Specifically, the method can be divided into three parts. The first part is the initialization of the coupled material point finite difference finite element method parallel simulation. The second part is the coupled material point finite difference finite element method parallel simulation and communication. The third part is the movement of the objects (material points, rod elements) in the coupled material point finite difference finite element method parallel simulation.

[0033] In the first part, the initialization process first needs to build an octree and use the octree encryption to complete the generation of the background grid. Secondly, based on the initial octree, the background grid is divided, and each process obtains the background grid it manages. Each process initializes objects (material points and rod units) and materials through its own managed grid physical area; finally, the background grid is classified according to whether it contains material points or fluid grids.

[0034] In the second part, the background grid is divided into three categories: fluid, solid, and fluid-solid coupling. Finite difference and solid methods (material point method, finite element method) are used for fluid grids and solid grids respectively, and immersed boundary method is used for fluid-solid coupling grids. A neighbor list in 26 directions is established. If it is in the process, it is searched in this process. If it is not in this process, the subscript of the target grid in the adjacent process is searched, and the data is packaged and sent. A data transmission method is designed for the special storage structure of the rod unit.

[0035] In the third part, the newly calculated positions of the material points and rod units are used to determine whether they are still in the octree nodes that this process is responsible for. If not, all their attributes are packaged and moved to the corresponding adjacent processes through the lists of material point and rod unit queries.

[0036] For example, Figure 1 FIG. 1 is a flow chart of a parallel coupled material point finite difference finite element simulation method provided in an embodiment of the present application. Figure 1 As shown, the method may include the following steps:

[0037] S11: Generate background grids, each of which includes a number of material points, and initialize material information of the material points, which includes mass, momentum, density, stress strain, deformation gradient, etc.

[0038] In this embodiment, an octree is generated according to the needs of the physical background grid. Among them, the physical background grid is a discretized representation of the simulated physical scene in space. It carries physical information such as material points, rod units and different materials. Its characteristics and distribution will directly affect the accuracy of the simulation results. The octree is a data structure used for space division. The tree is divided into three directions, and the product of the three directions is the number of original trees. Through the octree, the root node and adjacent relationships are generated, and the background grid is encrypted. The space sub-regions are divided in parallel according to the background grid, and each sub-region grabs its own material points. Initialize the material information of each material point, including mass, momentum, density, stress strain, deformation gradient, etc. At the same time, the properties of the fluid grid are added to the background grid to classify the grid units. Grid units can be divided into three types, namely solid grid units, fluid grid units, and fluid-solid two-phase grid units.

[0039] S12: Initialize the rod unit.

[0040] In this embodiment, in the finite element method, the rod unit is a basic unit used to simulate the behavior of structural mechanics, used to simulate the slender components subjected to axial force, and to simplify the simulation of similar actual components in structural mechanics analysis. It consists of two nodes and straight rods connecting the nodes. The nodes are located at both ends of the rod and are connected to other units through the nodes to construct a complex structural model. The rod unit is assumed to only bear axial tension or pressure, and the effects of shear force, bending moment and torque are not considered. When an external force acts on the rod unit, the force is transmitted along the axis of the rod, causing the rod unit to produce tensile or compressive deformation. The displacement mode of the rod unit is usually described by a linear function, that is, it is assumed that the displacement in the rod unit changes linearly from one node to another. This simple displacement mode can not only meet the calculation accuracy requirements, but also greatly simplify the calculation process. Each rod unit has specific material properties, such as elastic modulus (reflecting the ability of the material to resist elastic deformation), cross-sectional area, etc. These properties determine the mechanical response of the rod unit when it is subjected to force.

[0041] Unlike material points, rod units are composed of rod structures and point structures. The point structure is similar to material points and participates in the mapping and calculation of material point method grid points. The rod structure updates other physical quantities based on the position information of the two end points of the rod unit. When initializing the rod unit, the structural parameters are set for the rod unit. The structural parameters include the length, cross-sectional shape and size, elastic modulus, Poisson's ratio information, etc. of the rod unit, and the mapping relationship between the point structure and the material point of the rod unit is established. Since the rod unit nodes will move across processes with the calculation, it is difficult to uniquely represent the mapping relationship by only representing the mapping of the rod unit rod structure to the point structure. Therefore, a global index can be added to the original rod unit data structure, and a mapping list can be used to establish a mapping from the global index of the point to the local index in the local rod unit node array. While ensuring that the mapping relationship of the rod unit structure is unique, the relationship between the original rod unit rod structure and the point structure remains unchanged.

[0042] For example, please refer to Figure 2 , Figure 2 Figure 2 shows a schematic diagram of the rod unit division. Figure 2 As shown in the figure, when the rod unit is initially divided, since the rod unit rod structure has no position information, only the rod unit nodes have position information. Therefore, during initialization, the nodes at both ends of the rod unit are used to determine whether the rod is in the current process. Similar to the material point method, when both nodes of the rod unit are in the grid governed by the current process, the rod unit structure is considered to be in the current process, and the nodes at both ends are real points; if only one end is in the current process, the rod unit is considered to cross the process, and the rod unit has a copy in each of the two processes and the nodes at both ends are a pair of virtual and real nodes.

[0043] S13: Based on the finite difference finite element method, parallel simulation and communication are performed between each grid.

[0044] In this embodiment, data is transmitted between each grid. Data transmission can exchange information related to calculations between different processes. For example, when calculating different types of grids such as fluids, solids, and fluid-solid coupling, grid data, rod unit data, etc. are transmitted between processes, which can ensure the accuracy and coordination of calculations. Each process is responsible for the finite difference method and material point method related calculations of the background grid in the octree node it manages, and the data transmission with the neighbor grid is completed through octree management. Among them, finite differences are used to calculate fluid dynamics related equations for fluid grids, material point method or finite element method is used to calculate solid mechanics related equations for solid grids, and immersed boundary method is used to process interactions for fluid-solid coupling grids. A neighbor list mechanism is established. When data transmission involves cross-process, the subscript of the target grid in the adjacent process is searched according to the neighbor list, and the data is sent through data packaging and transmission protocol.

[0045] The calculation of the rod element involves two parts. The calculation of the rod element node structure is similar to the material point and is calculated independently, while the calculation of the rod structure is related to the nodes at both ends of the rod element. To facilitate the parallel rod element calculation, the virtual and real attributes of the end nodes are added to the rod element across the node. The virtual nodes do not participate in the calculation related to the end nodes, but are synchronized before the calculation related to the rod structure and participate in the calculation related to the rod structure. The real point has only one and unique copy in all processes, while the virtual point needs to synchronize the data of the corresponding real point before updating the stress and strain of the rod element.

[0046] In addition to being consistent with the material point method octree parallelism, it is necessary to synchronize the contributions of each process to the grid points when calculating physical quantities. It is also necessary to synchronize the data of the corresponding real points before updating the stress and strain of the rod element to maintain the consistency of the calculation results. For example, Figure 3 This is a schematic diagram of a rod unit data transmission provided in an embodiment of the present application. Figure 3 As shown, the embodiment of the present application designs a virtual-real node synchronization function to communicate the cross-process rod unit data to the corresponding rod unit node of the corresponding process. The synchronized rod unit node knows which process to receive from, which process to send to, which local rod unit node data to synchronize, and then sends / receives its data.

[0047] S14: Moving mesh data, including the movement of material points and the movement of rod elements.

[0048] In this embodiment, the processes to which material points and rod units belong may change in physical space due to calculation results, so it is necessary to transfer their positions between processes. This process not only involves the transmission of data, but also includes operations such as updating their positions in the new process and maintaining related indexes and structural relationships.

[0049] Specifically, a list of material point queries is established to query which object the material points contained in the grid to be moved belong to and their serial numbers in the object. The moving grid traverses the neighbor process list to know which process to send to and from which process to receive, locally counts the number of material points sent / received, performs a handshake to allocate corresponding memory, and then sends / receives the grid and the material points it carries.

[0050] For example, please refer to Figure 4 , Figure 4 Schematic diagram of the movement of a rod unit between processes provided by an embodiment of the present application. Figure 4 As shown in the figure, when the rod unit moves across processes, since the rod structure of the rod unit is highly bound to the point structure, but the rod structure does not have position information, the rod unit is moved across processes based on the rod structure and the position information of the nodes at both ends of the rod. When moving, in addition to moving the rod structure and the point structure of the rod unit to the corresponding process at the same time, it is also necessary to maintain the mapping in the index mapping list and update the nodes at both ends of the rod to maintain the uniqueness of the mapping relationship of the rod unit structure.

[0051] Steps S13 and S14 are executed in a loop until the end.

[0052] The above is an introduction to the parallel coupled material point finite difference finite element simulation method provided by the embodiment of the present application. By generating background grids, each background grid contains a number of material points, and the material information such as the mass and momentum of the material points is initialized. At the same time, the rod unit is also initialized. Then, based on the finite difference finite element method, parallel simulation and communication are carried out between each background grid. In this process, the background grids are divided into three categories: fluid, solid, and fluid-solid coupling. The corresponding calculation methods are respectively adopted, and the cross-process data transmission is processed by establishing a neighbor list mechanism. Finally, the mesh data is moved. For material points and rod units, it is judged whether it exceeds the scope of the current process according to their newly calculated positions. If it exceeds, it is moved. The material points package and transfer the material information. The rod unit ensures that the rod structure and the point structure are transferred synchronously when moving and updates the index mapping relationship to ensure the accuracy and integrity of the entire simulation process. This method provides an effective technical means for material explosion simulation and improves the accuracy of the simulation.

[0053] It is understandable that the size of the sequence number of each step in the above-mentioned embodiments does not mean the order of execution, and the execution order of each process should be determined by its function and internal logic, and should not constitute any limitation on the implementation process of the embodiments of the present application. In addition, in some possible implementations, the steps in the above-mentioned embodiments can be selectively executed according to actual conditions, and can be partially executed or fully executed, which is not limited here. All or part of any feature of any embodiment of the present application can be freely and arbitrarily combined without contradiction. The combined technical solution is also within the scope of the present application.

[0054] Based on the method in the above embodiment, the embodiment of the present application also provides a parallel coupled material point finite difference finite element simulation device. For example, Figure 5 A parallel coupled material point finite difference finite element simulation device is shown, which is deployed on a computing device, and the computing device is used for material explosion simulation. Figure 5 As shown, the device 500 includes an initialization module 501 and a processing module 502 .

[0055] Among them, the initialization module 501 is used to generate background grids, each background grid includes a number of material points, initialize the material information of the material points, and the material information includes at least one of mass, momentum, density, stress strain, and deformation gradient information; the initialization module 501 is also used to initialize the rod unit. The processing module 502 is used to perform parallel simulation and communication between various background grids based on the finite difference finite element method; the processing module 502 is also used to move grid data, including the movement of material points and the movement of rod units.

[0056] In some possible embodiments, the initialization module 501 is specifically used to: generate a background grid using an octree structure to construct the background grid, determine the final shape of the background grid by encrypting the octree, and the initial parameters and encryption rules of the octree are determined according to the spatial characteristics and accuracy of the material explosion simulation scene.

[0057] In some possible embodiments, the initialization module 501 is also used to: set structural parameters for the rod unit, the structural parameters including the length, cross-sectional shape and size, elastic modulus, and Poisson's ratio information of the rod unit, establish a mapping relationship between the point structure of the rod unit and the material point, wherein the point structure participates in the material point method grid point mapping and calculation, and the rod structure is updated based on the position information of the two end points of the rod unit and the set structural parameters.

[0058] In some possible embodiments, the processing module 502 is used to: divide the background grid into three categories: fluid, solid, and fluid-solid coupling; use finite difference to calculate fluid dynamics related equations for the fluid grid; use material point method or finite element method to calculate solid mechanics related equations for the solid grid; use immersed boundary method to handle interaction for the fluid-solid coupling grid; establish a neighbor list mechanism, when data transmission involves cross-process, search for the subscript of the target grid in the adjacent process based on the neighbor list, and send the data through data packaging and transmission protocol.

[0059] In some possible embodiments, the processing module 502 is also used for: for the movement of material points, by establishing a material point position query function and a process attribution judgment function, comparing the newly calculated position of the material point with the background grid range that the current process is responsible for, when the material point exceeds the current process range, packaging all the material information of the material point in a preset data format and moving it to the corresponding adjacent process; for the movement of the rod unit, judging whether it needs to be moved based on the position change of the two end points of the rod unit, if it needs to be moved, in addition to moving the rod structure and point structure of the rod unit to the corresponding process at the same time during the movement; updating the index mapping relationship of the rod unit in the new process.

[0060] It should be understood that the above-mentioned device is used to execute the method in the above-mentioned embodiment. The implementation principle and technical effect of the corresponding program module in the device are similar to those described in the above-mentioned method. The working process of the device can refer to the corresponding process in the above-mentioned method, which will not be repeated here.

[0061] The present application also provides a computing device 600. Figure 6 As shown, computing device 600 includes: bus 602, processor 604, memory 606 and communication interface 608. Processor 604, memory 606 and communication interface 608 communicate through bus 602. Computing device 600 can be a server or a terminal device. It should be understood that the present application does not limit the number of processors and memories in computing device 600.

[0062] The bus 602 may be a peripheral component interconnect (PCI) bus or an extended industry standard architecture (EISA) bus. The bus may be divided into an address bus, a data bus, a control bus, etc. For ease of representation, Figure 6 The bus 604 is represented by only one line, but does not mean that there is only one bus or one type of bus. The bus 604 may include a path for transmitting information between various components of the computing device 600 (eg, the memory 606, the processor 604, and the communication interface 608).

[0063] The processor 604 may include any one or more of a central processing unit (CPU), a graphics processing unit (GPU), a microprocessor (MP), or a digital signal processor (DSP).

[0064] The memory 606 may include a volatile memory, such as a random access memory (RAM). The processor 604 may also include a non-volatile memory, such as a read-only memory (ROM), a flash memory, a hard disk drive (HDD), or a solid state drive (SSD).

[0065] The memory 606 stores executable program codes, and the processor 604 executes the executable program codes to respectively implement the functions of the initialization module 501 and the processing module 502, thereby implementing all or part of the steps of the method in the above embodiment. That is, the memory 606 stores instructions for executing all or part of the steps in the method in the above embodiment.

[0066] Alternatively, the memory 606 stores executable codes, and the processor 604 executes the executable codes to respectively implement the functions of the aforementioned parallel coupled material point finite difference finite element simulation device 500, thereby implementing all or part of the steps in the above-mentioned embodiment method. That is, the memory 606 stores instructions for executing all or part of the steps in the above-mentioned embodiment method.

[0067] The communication interface 608 uses a transceiver module such as, but not limited to, a network interface card or a transceiver to implement communication between the computing device 600 and other devices or communication networks.

[0068] Based on the methods in the above embodiments, an embodiment of the present application provides a computer-readable storage medium, which stores a computer program. When the computer program runs on a processor, the processor executes the methods in the above embodiments.

[0069] Based on the methods in the above embodiments, an embodiment of the present application provides a computer program product. When the computer program product runs on a processor, the processor executes the methods in the above embodiments.

[0070] It is understood that the processor in the embodiments of the present application may be a central processing unit (CPU), or other general-purpose processors, digital signal processors (DSP), application specific integrated circuits (ASIC), field programmable gate arrays (FPGA) or other programmable logic devices, transistor logic devices, hardware components or any combination thereof. The general-purpose processor may be a microprocessor or any conventional processor.

[0071] The method steps in the embodiments of the present application can be implemented by hardware or by a processor executing software instructions. The software instructions can be composed of corresponding software modules, and the software modules can be stored in random access memory (RAM), flash memory, read-only memory (ROM), programmable read-only memory (PROM), erasable programmable read-only memory (EPROM), electrically erasable programmable read-only memory (EEPROM), registers, hard disks, mobile hard disks, CD-ROMs, or any other form of storage medium known in the art. An exemplary storage medium is coupled to the processor so that the processor can read information from the storage medium and write information to the storage medium. Of course, the storage medium can also be a component of the processor. The processor and the storage medium can be located in an ASIC.

[0072] In the above embodiments, it can be implemented in whole or in part by software, hardware, firmware or any combination thereof. When implemented using software, it can be implemented in whole or in part in the form of a computer program product. The computer program product includes one or more computer instructions. When the computer program instructions are loaded and executed on a computer, the process or function described in the embodiment of the present application is generated in whole or in part. The computer may be a general-purpose computer, a special-purpose computer, a computer network, or other programmable device. The computer instructions may be stored in a computer-readable storage medium or transmitted through the computer-readable storage medium. The computer instructions may be transmitted from a website site, computer, server or data center to another website site, computer, server or data center by wired (e.g., coaxial cable, optical fiber, digital subscriber line (DSL)) or wireless (e.g., infrared, wireless, microwave, etc.) means. The computer-readable storage medium may be any available medium that a computer can access or a data storage device such as a server or data center that includes one or more available media integrated. The available medium may be a magnetic medium (e.g., a floppy disk, a hard disk, a tape), an optical medium (e.g., a DVD), or a semiconductor medium (e.g., a solid state disk (SSD)), etc.

[0073] It should be understood that the various numerical numbers involved in the embodiments of the present application are only used for the convenience of description and are not used to limit the scope of the embodiments of the present application.

Claims

1. A parallel coupled material point finite difference finite element simulation method, applied to material explosion simulation, characterized in that: The method comprises: Generate a background grid, each background grid includes a plurality of material points, initialize material information of the material points, the material information includes at least one of mass, momentum, density, stress strain, and deformation gradient information; Initialize the rod unit; Based on the finite difference finite element method, parallel simulation and communication are performed between the background grids; The moving mesh data includes the movement of the material points and the movement of the rod elements.

2. The method according to claim 1, characterized in that The generating of the background grid specifically comprises: The generated background grid is constructed using an octree structure, and the final form of the background grid is determined by an encryption operation on the octree. The initial parameters and encryption rules of the octree are determined according to the spatial characteristics and accuracy of the material explosion simulation scene.

3. The method according to claim 1, characterized in that The initialization rod unit comprises: Structural parameters are set for the rod unit, wherein the structural parameters include the length, cross-sectional shape and size, elastic modulus, and Poisson's ratio information of the rod unit, and a mapping relationship between the point structure of the rod unit and the material point is established, wherein the point structure participates in the material point method grid point mapping and calculation, and the rod structure is updated according to the position information of the two end points of the rod unit and the set structural parameters.

4. The method according to claim 1, characterized in that: The parallel simulation and communication between the background grids based on the finite difference finite element method include: The background grid is divided into three categories: fluid, solid, and fluid-solid coupling. Finite difference is used to calculate fluid dynamics related equations for fluid grids, material point method or finite element method is used to calculate solid mechanics related equations for solid grids, and immersed boundary method is used to deal with interactions for fluid-solid coupling grids. A neighbor list mechanism is established. When data transmission involves cross-processes, the subscript of the target grid in the adjacent process is searched according to the neighbor list, and the data is sent through data packaging and transmission protocol.

5. The method according to claim 1, characterized in that The mobile grid data includes: For the movement of the material point, by establishing a material point position query function and a process attribution judgment function, the newly calculated position of the material point is compared with the background grid range that the current process is responsible for. When the material point exceeds the range of the current process, all the material information of the material point is packaged and moved to the corresponding adjacent process according to a preset data format; Regarding the movement of the rod unit, whether it needs to be moved is determined based on the position changes of the two end points of the rod unit. If it needs to be moved, in addition to moving the rod structure and point structure of the rod unit to the corresponding process at the same time, the index mapping relationship of the rod unit in the new process is updated.

6. A parallel coupled material point finite difference finite element simulation device, deployed on a computing device, the computing device is applied to material explosion simulation, characterized in that: The device comprises: An initialization module, used to generate a background grid, each background grid includes a plurality of material points, and initialize material information of the material points, wherein the material information includes at least one of mass, momentum, density, stress strain, and deformation gradient information; The initialization module is also used to initialize the rod unit; A processing module, used for performing parallel simulation and communication between the background grids based on the finite difference finite element method; The processing module is further used to move the mesh data, including the movement of the material points and the movement of the rod units.

7. The device according to claim 1, characterized in that The initialization module is specifically used for: The generated background grid is constructed using an octree structure, and the final form of the background grid is determined by an encryption operation on the octree. The initial parameters and encryption rules of the octree are determined according to the spatial characteristics and accuracy of the material explosion simulation scene.

8. The device according to claim 1, characterized in that The initialization module is also used for: Structural parameters are set for the rod unit, wherein the structural parameters include the length, cross-sectional shape and size, elastic modulus, and Poisson's ratio information of the rod unit, and a mapping relationship between the point structure of the rod unit and the material point is established, wherein the point structure participates in the material point method grid point mapping and calculation, and the rod structure is updated according to the position information of the two end points of the rod unit and the set structural parameters.

9. The device according to claim 1, characterized in that The processing module is used for: The background grid is divided into three categories: fluid, solid, and fluid-solid coupling. Finite difference is used to calculate fluid dynamics related equations for fluid grids, material point method or finite element method is used to calculate solid mechanics related equations for solid grids, and immersed boundary method is used to deal with interactions for fluid-solid coupling grids. A neighbor list mechanism is established. When data transmission involves cross-processes, the subscript of the target grid in the adjacent process is searched according to the neighbor list, and the data is sent through data packaging and transmission protocol.

10. The device according to claim 1, characterized in that The processing module is also used for: For the movement of the material point, by establishing a material point position query function and a process attribution judgment function, the newly calculated position of the material point is compared with the background grid range that the current process is responsible for. When the material point exceeds the range of the current process, all the material information of the material point is packaged and moved to the corresponding adjacent process according to a preset data format; Regarding the movement of the rod unit, whether it needs to be moved is determined based on the position changes of the two end points of the rod unit. If it needs to be moved, in addition to moving the rod structure and point structure of the rod unit to the corresponding process at the same time, the index mapping relationship of the rod unit in the new process is updated.

Images