A parallel algorithm for unilateral contact problems

•This contact algorithm is designed for High-Performance Computing.•The multicode approach allows using different time integration schemes for each body.•The size of the linear system of equations remains fixed.•The mesh partitioning is performed once at the beginning of the simulation.•The contact...

Celý popis

Uloženo v:
Podrobná bibliografie
Vydáno v:Computers & structures Ročník 271; s. 106862
Hlavní autoři: Guillamet, G., Rivero, M., Zavala-Aké, M., Vázquez, M., Houzeaux, G., Oller, S.
Médium: Journal Article
Jazyk:angličtina
Vydáno: Elsevier Ltd 15.10.2022
Témata:
Contact mechanics
Dirichlet–Neumann boundary conditions
Domain decomposition approach
Finite element method
High performance computing
Finite element method
Dirichlet–Neumann boundary conditions
High performance computing
Domain decomposition approach
Contact mechanics
ISSN:0045-7949
On-line přístup:Získat plný text
Tagy: Přidat tag
Žádné tagy, Buďte první, kdo vytvoří štítek k tomuto záznamu!
Popis
Shrnutí:•This contact algorithm is designed for High-Performance Computing.•The multicode approach allows using different time integration schemes for each body.•The size of the linear system of equations remains fixed.•The mesh partitioning is performed once at the beginning of the simulation.•The contact algorithm is simple and flexible for large-scale problems. In this paper, we introduce a novel parallel contact algorithm designed to run efficiently in High Performance Computing based supercomputers. Particular emphasis is put on its computational implementation in a multiphysics finite element code. The algorithm is based on the method of partial Dirichlet–Neumann boundary conditions and is capable to solve numerically a nonlinear contact problem between rigid and deformable bodies in a whole parallel framework. Its distinctive characteristic is that the contact is tackled as a coupled problem, in which the contacting bodies are treated separately, in a staggered way. Then, the coupling is performed through the exchange of boundary conditions at the contact interface following a Gauss–Seidel strategy. To validate this algorithm we conducted several benchmark tests by comparing the proposed solution against theoretical and other numerical solutions. Finally, we evaluated the parallel performance of the proposed algorithm in a real impact test to show its capabilities for large-scale applications.
ISSN:0045-7949
DOI:10.1016/j.compstruc.2022.106862