PeriFast/Dynamics: A MATLAB Code for Explicit Fast Convolution-based Peridynamic Analysis of Deformation and Fracture

We present PeriFast/Dynamics, a compact and user-friendly MATLAB code for fast peridynamic (PD) simulations for deformation and fracture. PeriFast/Dynamics uses the fast convolution-based method (FCBM) for spatial discretization and an explicit time marching scheme to solve large-scale dynamic fracture problems. Different from existing PD solvers, PeriFast/Dynamics does not require neighbor search and storage, due to the use of the Fast-Fourier transform and its inverse to compute the integral operator. Run-times and memory allocation are independent of the number of neighbors inside the PD horizon, leading to faster computations and lower storage requirements. The governing equations and discretization method are briefly reviewed, the code structure explained, and individual modules described in detail. A demonstrative example on dynamic brittle fracture in 3D, with multiple crack branching events, is solved using three different constitutive models: a bond-based, an ordinary state-based, and a correspondence model. The small differences between results with the three different constitutive models are explained. Users are provided with a step-by-step description of the problem setup and execution of the code. PeriFast/Dynamics is a branch of the PeriFast suite of codes, and is available for download at the GitHub link provided in reference [1].