Dual-layer hybrid parallel optimal transportation meshfreep (pOTM) method for solving large deformation of material

Abstract

The invention relates to a dual-layer hybrid parallel optimal transportation meshfreep method (pOTM) for solving large deformation of a material. The pOTM is a large-scale parallel method for a dual-layer hybrid parallel architecture based on distributed multithreading parallelization and shared memory multithreading parallelization; a set of Shadow point / node calculating data division scheme is developed on an MPI (Message Passing Interface) layer, and OTM simulation is subjected to distributed computing decomposition; on a multithreading layer, a multithreading parallelization computing scheme with a balanced dynamic load is developed. According to the pOTM scheme, a heterogeneous supercomputing cluster is effectively utilized to perform linear acceleration or even superlinear acceleration, so that the computing efficiency is greatly improved.

Description

technical field [0001] The invention relates to an efficient and stable numerical method for solving continuum mechanics problems involving extreme deformation, dynamic crack growth, high-speed impact and geometric distortion, material fission, metal material forming, multi-phase transformation and material dynamic response related to loading history, In particular, the meshless solution method for general continuum mechanics problems is a two-layer mixed parallel optimal transport meshless computing method that combines distributed multi-process parallelization and shared memory multi-thread parallelization. field of computational mechanics. Background technique [0002] Traditional numerical methods based on Lagrangian grids and Eulerian grids, such as finite element method, finite volume method, and finite difference method, etc., limit their performance in extremely large deformation and dynamic crack propagation due to their inherent grid characteristics. , high-speed ...

AI Technical Summary

Problems solved by technology

However, the development of mesh-free methods is also facing difficulties that need to be solved urgently: ① most of the mesh-free interpolation functions do not satisfy the Kronecker property, making it very difficult to apply displacement boundary conditions; ② in the absence of mesh Under the condition of the gridless Galerkin method, it will be more difficult to carry out the equivalent integration in the "weak" form. A large number of research literatures show that several possible integration schemes (node ​​integration scheme, stress point integration scheme , background grid integration schemes) are faced with difficulties in dealing with tensile stress instability and numerical integration errors to varying degrees; ③lack of rigorous mathematical analysis of convergence and error theory; ④due to the use of high-order interpolation in meshless methods function, which poses a great challenge to high-performance computing

Images