Parallel implementation of implicit finite element model with cohesive zones and collision response using CUDA
Authors:
- 11 May 2018
- journal article
- research article
- Published by Wiley in International Journal for Numerical Methods in Engineering
- Vol. 115 (7) , 771-790
Abstract
Summary: The aim of this work is to efficiently implement the Park‐Paulino‐Roesler cohesive zone model with the objective of creating realistic high‐resolution simulations of material deformation, fracture, and postfracture behavior. Intrinsically, unstructured meshes can create more realistic fracture patterns in bulk material than structured meshes. Implicit methods, stable for much larger time steps, have greater potential to model both fracture and postfracture behavior without sacrificing speed of execution. Several technical contributions are presented, including (i) GPU‐accelerated implementation of the Park‐Paulino‐Roesler cohesive zone model, (ii) efficient creation of sparse matrix structure, and (iii) comparison of different unloading/reloading relations when using an implicit scheme. A potential‐based collision response scheme was implemented that allows one to model the interaction of fragmented material. Several test simulations are carried out to demonstrate the flexibility of the model and its ability to reproduce different materials under various loading conditions. Benchmarking results show that most of the computational time is spent by the third‐party solver library, meaning that other operations do not require optimization. The library is made available as open source.Keywords
Funding Information
- Natural Sciences and Engineering Research Council of Canada
- Research and Development Corporation of Newfoundland and Labrador
This publication has 23 references indexed in Scilit:
- Vega: Non‐Linear FEM Deformable Object SimulatorComputer Graphics Forum, 2012
- Computational implementation of the PPR potential-based cohesive model in ABAQUS: Educational perspectiveEngineering Fracture Mechanics, 2012
- Implicit FEM Solver on GPU for Interactive Deformation Simulation2012
- Large-scale 3D random polycrystals for the finite element method: Generation, meshing and remeshingComputer Methods in Applied Mechanics and Engineering, 2011
- Assembly of finite element methods on graphics processorsInternational Journal for Numerical Methods in Engineering, 2010
- MCCD: Multi-core collision detection between deformable models using front-based decompositionGraphical Models, 2010
- A unified potential-based cohesive model of mixed-mode fractureJournal of the Mechanics and Physics of Solids, 2009
- Conserving energy and momentum in nonlinear dynamics: A simple implicit time integration schemeComputers & Structures, 2007
- Energy conservation in Newmark based time integration algorithmsComputer Methods in Applied Mechanics and Engineering, 2006
- The Strength and Ductility of Ice Under TensionJournal of Offshore Mechanics and Arctic Engineering, 1988