Mapping cohesive fracture and fragmentation simulations to graphics processor units.

*(English)*Zbl 1352.65671Summary: A graphics processor units(GPU)-based computational framework is presented to deal with dynamic failure events simulated by means of cohesive zone elements. The work is divided into two parts. In the first part, we deal with pre-processing of the information and verify the effectiveness of dynamic insertion of cohesive elements in large meshes in parallel. To this effect, we employ a novel and simplified topological data structure specialized for meshes with triangles, designed to run efficiently and minimize memory occupancy on the GPU. In the second part, we present a parallel explicit dynamics code that implements an extrinsic cohesive zone formulation where the elements are inserted ‘on-the-fly’, when needed and where needed. The main challenge for implementing a GPU-based computational framework using an extrinsic cohesive zone formulation resides on being able to dynamically adapt the mesh, in a consistent way, by inserting cohesive elements on fractured facets. In order to handle that, we extend the conventional data structure used in the finite element method (based on element incidence) and store, for each element, references to the adjacent elements. This additional information suffices to consistently insert cohesive elements by duplicating nodes when needed. Currently, our data structure is specialized for triangular meshes, but an extension to tetrahedral meshes is feasible. The data structure is effective when used in conjunction with algorithms to traverse nodes and elements. Results from parallel simulations show an increase in performance when adopting strategies such as distributing different jobs among threads for the same element and launching many threads per element. To avoid concurrency on accessing shared entities, we employ graph coloring. In a pre-processing phase, each node of the dual graph (bulk elements of the mesh as graph nodes) is assigned a color different from the colors assigned to adjacent nodes. In that fashion, elements of the same color can be processed in parallel without concurrency. All the procedures needed for the insertion of cohesive elements along fracture facets and for computing nodal properties are performed by threads assigned to triangles, invoking one kernel per color. Computations on existing cohesive elements are also performed based on adjacent bulk elements. Experiments show that GPU speedup increases with the number of nodes and bulk elements.

PDF
BibTeX
XML
Cite

\textit{A. Alhadeff} et al., Int. J. Numer. Methods Eng. 103, No. 12, 859--893 (2015; Zbl 1352.65671)

Full Text:
DOI

##### References:

[1] | Seegyoung Seol, Efficient distributed mesh data structure for parallel automated adaptive analysis, Engineering with Computers 22 (3) pp 197– (2006) · Zbl 05192774 |

[2] | Kirk, Libmesh: a c++ library for parallel adaptive mesh refinement/coarsening simulations, Engineering with Computers 22 (3) pp 237– (2006) · Zbl 05192775 |

[3] | Klein, Physics-based modeling of brittle fracture: cohesive formulations and the application of meshfree methods, Theoretical and Applied Fracture Mechanics 37 (1-3) pp 99– (2001) |

[4] | Ortiz, Finite-deformation irreversible cohesive elements for three-dimensional crack-propagation analysis, International Journal for Numerical Methods in Engineering 44 (9) pp 1267– (1999) · Zbl 0932.74067 |

[5] | Cirak, A cohesive approach to thin-shell fracture and fragmentation, Computer Methods in Applied Mechanics and Engineering 194 (21-24) pp 2604– (2005) · Zbl 1082.74052 |

[6] | Zhang, Extrinsic cohesive modelling of dynamic fracture and microbranching instability in brittle materials, International Journal for Numerical Methods in Engineering 72 (8) pp 893– (2007) · Zbl 1194.74319 |

[7] | Belytschko, Dynamic crack propagation based on loss of hyperbolicity and a new discontinuous enrichment, International Journal for Numerical Methods in Engineering 58 (12) pp 1873– (2003) · Zbl 1032.74662 |

[8] | Zhang, Cohesive zone modeling of dynamic failure in homogeneous and functionally graded materials, International Journal of Plasticity 21 (6) pp 1195– (2005) · Zbl 1154.74391 |

[9] | Sharon, Local crack branching as a mechanism for instability in dynamic fracture, Physical Review Letters 74 (25) pp 5096– (1995) |

[10] | Sharon, Microbranching instability and the dynamic fracture of brittle materials, Physical Review B (Condensed Matter) 54 (10) pp 7128– (1996) |

[11] | Dooley, Parallel simulations of dynamic fracture using extrinsic cohesive elements, Journal of Scientific Computing 39 (1) pp 144– (2009) · Zbl 1203.74004 |

[12] | Espinha, Partops: compact topological framework for parallel fragmentation simulations, Engineering with Computers 25 (4) pp 345– (2009) · Zbl 05781547 |

[13] | Arias, Numerical modelling and experimental validation of dynamic fracture events along weak planes, Computer Methods in Applied Mechanics and Engineering 196 pp 3833– (2007) · Zbl 1173.74394 |

[14] | Molinari, The cohesive element approach to dynamic fragmentation: the question of energy convergence, International Journal for Numerical Methods in Engineering 69 pp 484– (2007) · Zbl 1194.74450 |

[15] | Paulino, Adaptive dynamic cohesive fracture simulation using nodal perturbation and edge-swap operators, International Journal for Numerical Methods in Engineering 84 pp 1303– (2010) · Zbl 1202.74148 |

[16] | Lawlor, Parfum: a parallel framework for unstructured meshes for scalable dynamic physics applications, Engineering with Computers 22 (3) pp 215– (2006) · Zbl 05192771 |

[17] | Radovitzky, A scalable 3D fracture and fragmentation algorithm based on a hybrid, discontinuous Galerkin, cohesive element method, Computer Methods in Applied Mechanics and Engineering 200 (1-4) pp 326– (2011) · Zbl 1225.74105 |

[18] | Bolz, Sparse matrix solvers on the GPU: conjugate gradients and multigrid, ACM Transactions on Graphics 22 pp 917– (2003) · Zbl 05457243 |

[19] | Wu, A hybrid condensed finite element model with GPU acceleration for interactive 3D soft tissue cutting: research articles, Computer Animation and Virtual Worlds 15 (3-4) pp 219– (2004) · Zbl 02179830 |

[20] | Krakiwsky, Acceleration of finite-difference time-domain (FDTD) using graphics processor units (GPU), Microwave Symposium Digest, 2004 IEEE MTT-S International 2 pp 1033– (2007) |

[21] | Tejada, Large steps in GPU-based deformable bodies simulation, Simulation Modelling Practice and Theory 13 (8) pp 703– (2005) |

[22] | Taylor, High-speed nonlinear finite element analysis for surgical simulation using graphics processing units, IEEE Transactions on Medical Imaging 27 (5) pp 650– (2008) |

[23] | Göddeke, Using GPUs to improve multigrid solver performance on a cluster, International Journal of Computational Science and Engineering (IJCSE) 4 (1) pp 36– (2008) |

[24] | Anderson, General purpose molecular dynamics simulations fully implemented on graphics processing units, Journal of Computational Physics 227 (10) pp 5342– (2008) · Zbl 1148.81301 |

[25] | Rodriguez-Navarro J Susin A Non structured meshes for cloth GPU simulation using FEM 2006 1 7 |

[26] | Göddeke D Strzodka R Turek S Accelerating double precision FEM simulations with GPUs Proceedings of ASIM 2005 - 18th Symposium on Simulation Technique SCS Publishing House e.V., ASIM 2005 139 144 |

[27] | Cecka, Assembly of finite element methods on graphics processors, International Journal for Numerical Methods in Engineering 85 (5) pp 640– (2011) · Zbl 1217.80146 |

[28] | Geveler, Towards a complete FEM-based simulation toolkit on GPUs: unstructured grid finite element geometric multigrid solvers with strong smoothers based on sparse approximate inverses, Computers & Fluids 80 pp 327– (2012) · Zbl 1284.76249 |

[29] | Komatitsch, Porting a high-order finite-element earthquake modeling application to Nvidia graphics cards using CUDA, Journal of Parallel and Distributed Computing 69 (5) pp 451– (2009) |

[30] | Liu Y Jiao S Wu W De S GPU accelerated fast FEM deformation simulation IEEE Asia Pacific Conference on Circuits and systems, 2008. APCCAS 2008 Macao 2008 606 609 |

[31] | Fan Z Qiu F Kaufman A Yoakum-Stover S GPU cluster for high performance computing Proceedings of the 2004 ACM/IEEE Conference on Supercomputing SC ’04 Pittsburgh, PA 2004 12 |

[32] | Godel, Scalability of higher-order discontinuous Galerkin FEM computations for solving electromagnetic wave propagation problems on GPU clusters, IEEE Transactions on Magnetics 46 (8) pp 3469– (2010) |

[33] | Kakay, Speedup of FEM micromagnetic simulations with graphical processing units, IEEE Transactions on Magnetics 46 (6) pp 2303– (2010) |

[34] | Ren, Power aware parallel 3-D finite element mesh refinement performance modeling and analysis with CUDA/MPI on GPU and multi-core architecture, IEEE Transactions on Magnetics 48 (2) pp 335– (2012) |

[35] | Markall, Finite element assembly strategies on multi-core and many-core architectures, International Journal for Numerical Methods in Fluids 71 (1) pp 80– (2013) |

[36] | Zegard, Toward GPU accelerated topology optimization on unstructured meshes, Structural and Multidisciplinary Optimization 48 pp 473– (2013) |

[37] | Kirk, Programming Massively Parallel Processors: A Hands-On Approach (2010) |

[38] | Cuda c programming guide 3.2 2010 |

[39] | Sanders, CUDA by Example: An Introduction to General-Purpose GPU Programming (2010) |

[40] | Celes W, A compact adjacency-based topological data structure for finite element mesh representation, International Journal for Numerical Methods in Engineering 64 (11) pp 1529– (2005) · Zbl 1122.74504 |

[41] | Paulino, A general topology-based framework for adaptive insertion of cohesive elements in finite element meshes, Engineering with Computers 24 (1) pp 59– (2008) · Zbl 05535658 |

[42] | Welsh, An upper bound for the chromatic number of a graph and its application to timetabling problems, The Computer Journal 10 (1) pp 85– (1967) · Zbl 0147.15206 |

[43] | Pandolfi, Solid modeling aspects of three-dimensional fragmentation, Engineering with Computers 14 (4) pp 287– (1998) · Zbl 0933.68141 |

[44] | Park, A unified potential-based cohesive model of mixed-mode fracture, Journal of the Mechanics and Physics of Solids 57 (6) pp 891– (2009) |

[45] | Camacho, Computational modelling of impact damage in brittle materials, International Journal of Solids and Structures 33 (20) pp 2899– (1996) · Zbl 0929.74101 |

This reference list is based on information provided by the publisher or from digital mathematics libraries. Its items are heuristically matched to zbMATH identifiers and may contain data conversion errors. It attempts to reflect the references listed in the original paper as accurately as possible without claiming the completeness or perfect precision of the matching.